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INTRODUCTION 


The  Institute  for  Naval  Oceanography  (INO)  sponsored  Workshop  on  Data 
Assimilation  in  Coupled  Air-Sea  Nowcast/Forecast  Systems  was  held  Febmary  27 
and  28,  1992,  to  promote  discussions  and  develop  a  strategy  for  research  and 
development  of  data  assimilation  techniques  in  coupled  atmospheric/ocean 
nowcast/forecast  modeling  systems  for  the  Mediterranean  Sea.  The  motivation 
for  the  workshop,  in  terms  of  INO’s  program,  is  the  changing  priority  of  Navy 
program  sponsors  to  an  emphasis  in  coastal  and  semi-enclosed  seas  research  and 
development  efforts.  Within  this  area,  the  Mediterranean  Sea  has  been 
specifically  identified  by  the  Navy  as  priority. 

The  purpose  of  developing  a  strategy  for  air/sea  coupling  in  the 
Mediterranean  is  to  identify  opportunities  for  collaborative  efforts  in  basic 
research  with  regard  to  data  assimilation  and  for  eventual  development  of 
coupled  nowcast/forecast  systems  in  coastal  regions  and  marginal  seas,  beginning 
with  the  Mediterranean  Sea.  It  is  recognized  that  there  has  been  considerable 
efforts  in  gaining  a  fundamental  understanding  of  the  circulation  of  this  region 
through  previous  and  ongoing  observational  programs,  numerical  modeling 
efforts,  and  theoretical  work,  and  that  there  is  continuing,  strong  scientific 
interest  in  these  topics.  There  has  been  progress  in  applying  data  assimilation 
methods  to  ocean  and  atmospheric  models  independent  of  each  other.  These 
uncoupled  models  can  reasonably  operate  within  presently  available  computing 
resource  environments.  Further,  some  understanding  of  how  to  couple  air-sea 
models  exists,  and  there  is  recognition  at  this  time  that  an  approach  could  be 
developed  to  integrate  this  knowledge  into  a  design  for  a  coupled  air/sea 
nowcast/forecast  system. 

OBJECTIVES  OF  THE  WORKSHOP 

1 .  Conduct  open  forum  discussions  regarding  the  fundamental  scientific  issues 
in  data  assimilation  in  ocean  circulation  models,  in  atmospheric  models, 
and  in  coupled  ocean/atmosphere  nowcast/forecast  systems. 

2.  Develop  a  strategy  of  collaborative  research  for  approaches  in  addressing 
data  assimilation  issues  in  coupled  systems  (data  types  are  those 
“commonly”  observed  or  derived  from  routine  observations), 

3.  Design  a  developmental  plan  to  support  research  efforts  such  that  a 
research/development  coupled  nowcast/forecast  system  for  the 
Mediterranean  could  reside  fully  within  ECMOP. 
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4.  Produce  a  manuscript  of  workshop  results/conclusions. 

Hence,  the  participants  were  charged  to  propose  a  strategy  for 
coUaborative  air-sea  interaction  research,  and  to  outline  a  developmental  plan  for 
supporting  the  proposed  research.  A  creative  approach  to  addressing  the 
objectives  was  motivated  through  discussions  on  relevant  research  issues,  that 
articulated  the  technical  support  required,  for  the  design  of  a  coupled 
nowcast/forecast  system  in  the  Mediterranean  Sea.  To  that  end,  there  were 
presentations  on  the  known  existing  components  that  would  comprise  the  coupled 
system.  These  were  followed  by  working  group  and  plenary  discussions. 

WORKING  GROUPS 

Questions  for  the  working  groups  to  address  were  posed  prior  to  the 
workshop.  These  were  not  considered  to  be  all  inclusive  but  to  serve  as  a 
template  to  aid  in  meeting  the  workshop  objectives.  The  questions  were; 

1.  How  can  one  determine  the  resolution  needed  to  adequately  address 
baroclinic  effects  with  sigma  coordinates  given  Mediterranean  topography?  How 
can  one  determine  the  maximum  bottom  slope  that  can  be  used  in  a  numerical 
model  and  still  maintain  adequate  accuracy  for  the  baroclinic  terms  at  a  given 
resolution  and  typical  Mediterranean  buoyancy  distribution? 

2.  What  are  the  physical  process  scales  that  must  be  addressed  (i.e., 
cannot  be  well  represented  by  present  tractable  turbulence  closure  schemes)? 

3.  What  software  functionality  currently  used  by  this  community  should 
be  included  in  future  versions  of  ECMOP? 

4.  When  and  where  are  free-surface  effects  important  in  the 
Mediterranean? 

5.  What  available  data  and  data  assimilation  schemes  are  of  potential 
value  in  a  coupled  air-sea  system?  Which  of  these  are  applicable  for  the 
Mediterranean?  How  do  we  validate  data  assimilation  when  the  model  validation 
uses  the  same  data  types  for  both? 

6.  Are  present  mixing  and/or  instant  convective  adjustment  schemes 
adequate  for  addressing  Mediterranean  water  mass  transformations?  Are  there 
available  schemes  that  are  applicable  to  the  two-phase  wavy  air-sea  boundary 
layer?  (NOTE:  Observations  of  weakly  stratified  "chimneys"  in  the  North 
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Atlantic  suggest  that  such  schemes  might  not  be  good  for  the  North  Atlantic.  Do 
such  phenomena  occur  as  a  result  of  rapid  surface  cooling  through  sensible  heat 
exchange  or  evaporation  in  the  Mediterranean?  If  so,  how  important  are  they? 
Must  we  develop  new  methods  to  adequately  represent  them?) 

7.  Should  we  use  vertical  boundary-fitted  (sigma)  coordinates  in  order 
to  accurately  represent  bottom  boundary  layer  effects?  (An  alternative  is  to  use  a 
separate  bottom  boundary  layer  model  to  determine  vertical  velocity  and 
momentum,  heat,  and  salinity  fluxes  from  the  bottom  boundary  layer.) 

8.  What  are  the  important  technical  issues  that  need  to  be  addressed  to 
successfully  model  the  interaction  between  the  atmosphere  and  ocean? 

9.  Should  we  apply  lateral  boundary-fitted  coordinates  to  the 
Mediterranean?  If  so,  what  kind  of  grid(s)  is  desirable? 

The  group  discussions,  (objective  1  of  the  workshop),  lead  to  the 
development  of  a  Strategy  for  Collaborative  Research  and  a  design  for  the 
Developmental  Plan,  fulfilling  workshop  objectives  2  and  3.  This  report  fulfills 
objective  4. 

CONCLUSIONS  AND  RECOMMENDATIONS 

After  the  presentations  (see  Appendix  A),  two  working  groups  were 
formed.  Each  group  addressed  the  questions  posed  prior  to  the  workshop  plus 
any  additions  they  felt  were  needed.  The  original  charter  was  expanded  to 
include  fundamental  scientific  and  technical  issues  to  provide  a  more 
comprehensive  approach.  This  allowed  the  groups  to  also  address  data-related 
issues  explicitly.  Although  not  all  topics  were  covered  as  thoroughly  as  desired, 
the  following  conclusions  and  recommendations  were  made. 

•  Subgrid  scale  process  modeling  was  discussed  in  detail.  “Bulk” 
methods,  such  as  mixed  layers,  seem  more  compatible  with  layer 
models,  while  diffusive-based  methods  fit  more  naturally  into  level 
models.  It  was  felt  that  the  parameterizations  would  need  to 
incorporate  the  findings  of  non-hydrostatic  studies,  making  them 
subphysics  parameterizations  as  well.  In  the  Mediterranean,  it  is 
most  important  to  let  the  physical  dynamics  dictate  the 
parameterizations  needed  to  allow  for  the  episodic  nature  of  the  deep 
convection.  That  is,  there  is  a  need  to  choose  relevant  assumptions. 
Towards  this  end  it  was  felt  that  the  importance  of  pre-conditioning 
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events  for  the  episodic  deep  convection  events  be  given  visibility  in 
the  program. 

The  physics  of  the  deepening  processes,  such  as  those  that 
form  the  deep  water  and  the  Levintine  intermediate  water, 
need  to  be  clearly  identifiable  in  the  subgrid  scale 
parameterization  scheme  used.  Any  such  parameterized 
physics  must  be  consistent  throughout  all  interacting  models  being 
used,  such  as  a  coupled  general  circulation  model  (GCM)  and  a 
regional  model. 

Twin  experiments  should  be  mn  before  beginning  the  actual  data 
assimilation,  or  even  before  beginning  to  do  the  data  set  development 
and  collection,  for  several  reasons:  (a)  data  assimilation  schemes 
need  to  be  thoroughly  tested,  (b)  “value  added”  determination  is 
needed  for  various  data  types,  and  (c)  concerns  were  voiced 
regarding  the  availability  of  high-quality  satellite  data.  Although  the 
Earth  Resources  Satellite  (ERS-1)  is  providing  data,  the  errors  are 
still  too  large,  and  the  data  is  not  widely  available.  Although  existing 
atmospheric  models  have  water  vapor  correction  algorithms,  they 
are  not  at  sufficient  resolution  for  ocean  models.  It  was  suggested 
that  a  preliminary  study  could  be  made  using  the  Geodetic  Earth 
Orbiting  Satellite  (GEOS AT)  data  from  1987-1988,  under  the 
known  caveat  that  the  available  atmospheric  model  data  to 
complement  such  a  study  is  not  near  as  good  as  current  model  output 
due  to  coarser  resolution  and  “bad”  wind  data. 

Other  twin  experiments  can  be  run. 

(1)  Previous  wind  sensitivity  studies  could  be  reproduced, 
accounting  for  the  fact  that  the  ocean  needs  wind  stress  rather 
than  the  wind  field.  Research  is  needed  to  determine  the  wave 
effects,  because  the  stress  in  a  wave  model  is  not  the  exact 
same  stress  for  the  wind  driven  ocean  circulation.  A  wave 
model,  such  as  the  WAM,  would  need  to  be  used  to  specify  the 
percentage  of  stress  that  is  successfully  transferred  into  the 
ocean.  WAM  is  operational  in  the  Mediterranean  at  coarse 
resolution. 
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(2)  Tide,  expendable  bathythermograph  (XBT),  drifter,  and 
sea-surface  temperature  (SST)  twin  experiments  need  to  be 
run.  The  relationship  between  atmospheric  and  ocean  models 
in  terms  of  consistency,  is  unknown  in  both  performance  and 
effect.  For  example,  the  atmospheric  models  need  sea-surface 
skin  temperature,  while  ocean  models  need  finite  depth 
sea-surface  temperature,  and  the  relationship  between  the  two 
is  not  well  understood. 

(3)  There  is  also  a  need  for  three-dimensional  (3-D) 
assimilation,  such  as  can  be  generated  with  an  objective 
analysis  coupled  with  a  feature  model  (like  the  Optimal 
TTiermal  Interpolation  System  [OTIS]),  but  much  more  work 
is  needed  to  develop  a  product  that  meets  the  ocean  modeling 
needs,  and  no  thought  has  been  given  to  any  possible  coupled 
air-sea  interaction  needs. 

•  The  question  of  fluxes  between  the  atmospheric  and  ocean  models 
was  also  addressed.  The  different  resolution  needs  of  the  air  and  sea 
probably  will  not  affect  the  atmospheric  models,  but  probably  will 
affect  the  coupled  ocean  model,  because  it  needs  higher  accuracy. 
But  there  are  additional  existing  problems.  Flux  parameterizations 
approach  zero  as  the  wind  speed  approaches  zero.  And  in  the 
summer,  the  Mediterranean  can  reside  in  a  low  wind  speed  regime 
for  days  at  a  time.  During  such  episodes,  the  models  would  be 
completely  uncoupled.  Hence,  it  was  suggested  that  we  may  need  an 
intermediate  step  between  the  two  coupled  models  that  first  identifies 
which  wind  speed  regime  is  active  and  then  switches  to  an 
appropriate  scheme.  It  was  pointed  out  that,  at  the  mesoscales,  a 
1/4°C  change  in  temperature  can  make  the  difference  between  a 
stable  and  an  unstable  atmospheric  event;  hence,  the  resolution- 
feedback  interaction  may  significantly  affect  the  atmospheric 
forecast. 

•  A  free  surface  is  a  necessity  because,  although  the  tidal  variations 
overall  are  of  small  amplitude,  there  is  a  strong  tidal  effect  at 
Gibraltar  (see  location  of  tide  guage  stations  in  Figure  1).  Further, 
storm  surge  winds  can  play  a  large  dynamic  roll  in  the  region.  As  a 
side  effect  it  also  provides  a  more  elegant  way  to  assimilate  altimeter 
data. 
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•  Fluxes  between  the  atmosphere  and  ocean  seem  a  more  physically 
compelling  medium  of  exchange  than  SST,  and  the  choice  needs  to  be 
investigated  indepth.  It  was  recommended  that  the  regional  output 
from  the  Navy  Operational  Regional  Atmospheric  Prediction 
System  (NORAPS),  fluxes  included,  begin  to  be  saved  as  soon  as 
possible. 

•  Coupled  air-sea  models  must  compute  surface  fluxes  in  a  “consistent” 
manner.  A  caveat  was  added  stating  that  the  coupling  studies  should 
start  as  one-way  communication  only,  reproducing  current  research 
in  this  area. 

•  Some  current  atmospheric  models  average  over  three  hour  periods 
trying  to  incorporate  a  diurnal  signature,  while  others  try  to  average 
out  any  diurnal  variability.  Such  differences  need  to  be  considered 
when  choosing  the  atmospheric  model  to  use. 

•  Air-sea  interaction  consistency  would  dictate  a  finite  depth  SST,  but 
sometimes  the  atmospheric  models  may  want  a  true  skin 
temperature. 

•  A  major  area  of  research  entails  determining  the  interactive  air-sea 
processes  that  dictate  the  feedback  coupling  mechanisms  that  occur. 
Determining  the  necessity,  if  any,  of  additional  communication 
besides  SST  is  a  large  body  of  work.  The  ocean  modelers  currently 
believe  they  need  accurate  fluxes  and  wind  stress  (both  magnitude 
and  direction  to  obtain  lines  of  zero  curl)  to  begin  a  first  order 
attack.  Again,  it  is  unknown  if  additional  information  is  needed. 
This  is  the  major  fundamental  issue  around  which  a  strategy  for 
collaboration  needs  to  be  constructed. 

•  An  interesting  interaction  occurs  between  the  atmosphere  and  the 
ocean  during  storms  where  the  ocean  sea-surface  height  (SSH)  can 
vary  over  a  large  area,  effectively  varying  Zq. 

•  Sub-basin  to  mesoscale  interaction  is  more  important  in  the 
Mediterranean  than  in  the  open  ocean,  and  the  data  sets  must  reflect 
this. 
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A  resolution  study  for  the  Mediterranean  was  also  suggested  that 
would  determine  a  sufficient  sampling  strategy  for  the  Ekman  layer 
and  the  other  present  vertical  structures. 

Ocean  models  will  need  good  correlation  functions  for  the 
Mediterranean  for  many  reasons,  including  sparse  data  and  the  need 
for  3-D  fields. 

Ocean  models  need  better  local  wind  forcing  to  resolve  coastal 
upwelling  (downwelling).  It  is  unknown  what  will  be  a  sufficient 
resolution  for  the  very  shallow  coastal  surf  zone  forcing. 

Island  effects  are  crucial  in  the  Mediterranean  where  high  mountains 
can  strongly  affect  the  wind  field. 

Steep  topography  in  the  Mediterranean  can  be  solved  by 
“calibrating”  the  initial  background  density  field,  and  should  not 
introduce  any  new  numerical  problems. 

A  5  km  horizontal  resolution  combined  with  0(30)  vertical  levels 
would  be  desirable,  should  future  computer  hardware  make  this 
feasible. 

It  is  recommended  that  a  study  be  made  on  the  numerical 
characteristics,  including  truncation  error,  associated  with  different 
grids  including  rectilinear,  orthogonal  curvilinear,  and 
nonorthogonal  grids. 

Understand  and  learn  from  the  existing  experience  in  atmospheric 
modeling,  but  at  the  same  time,  avoid  archaic  practices.  This 
recommendation  can  be  taken  at  many  different  levels,  and  INO 
recommends  that  one  scientist  be  on  the  project  who  has  sufficient 
atmospheric  modeling  expertise  to  guide  the  remaining  ocean 
modelers  to  the  relevant  literature.  This  person  does  not  need  to  be 
an  INO  scientist;  they  may  be  a  consultant. 

In  some  regions  of  the  Mediterranean,  such  as  in  well  defined  basins, 
reduced  models  have  been  getting  good  results  (e.g..  Harvard 
quasi-geostrophic  model).  Although  this  is  a  longer  term  issue, 
coupling  to  simplified  models  for  computational  efficiency  may  be 
considered. 


•  The  need  for  better  planetary  boundary  layers  was  considered.  It 
was  indicated  that  there  is  a  need  to  move  from  first  order  closure 
through  one-and-a-half  order  to  second  order  closure  schemes. 

•  Modeling  mixed-layer  and  convective  processes  well  is  a  high 
priority,  but  this  is  not  a  compelling  reason  to  move  to 
“non-hydrostatic”  models.  Current  results  in  the  Mediterranean, 
with  the  existing  Princeton  model  turbulence  schemes,  already 
reproduce  many  aspects  of  convective  cells  (chimneys);  it  is  unclear 
whether  the  evaporation/precipitation  process  models  or  the  flux 
parameterization  values  are  responsible  for  the  current  discrepancies 
between  models  and  data.  Hence,  more  measurements  are  needed. 

•  Concluding  that  a  Kalman  filtering  assimilation  scheme  is  inherently 
limited  by  its  linear  nature,  it  was  suggested  that  both  adjoint  and 
optimal  interpolation  (01)  schemes  be  pursued  concurrently.  The  01 
work  would  involve  an  extensive  regional  effort  in  the 
Mediterranean  to  address  the  complex,  multiscale,  inhomogeneous, 
anisotropic  nature  of  some  of  the  dynamical  processes  in  the  region. 
In  the  larger  group,  it  was  recommended  to  always  let  the  model 
define  the  correlations,  but  Allan  Robinson  stressed  the  need  to  also 
consider  the  data  characteristics.  While  the  developed  01  scheme 
would  be  limited  by  the  research  results  for  the  region,  it 
would  nonetheless  work  fairly  well  for  many  flow  regimes  and  be 
relatively  efficient.  The  second  adjoint  approach  would  focus  on  an 
entirely  separate  problem:  how  to  simplify  the  method  to  make  it 
more  efficient  and,  hence,  tractable. 

•  ECMOP  needs  more  testing  and  use  by  the  ocean  and  atmospheric 
communities  before  recommendations  can  be  made  about  choosing  it 
as  the  Mediterranean  platform  for  future  development  work.  Due  to 
existing  limited  knowledge,  it  was  not  possible  to  dictate  the  software 
facilities  needed  in  ECMOP.  On  the  contrary,  it  was  suggested  that 
ECMOP  take  on  this  project  as  a  central  application  effort  instead, 
which  would  allow  the  software  support  to  be  put  together  on  an  “as 
needed”  basis. 

ISSUES  NOT  ADDRESSED 

In  workshops  many  more  issues  are  briefly  mentioned  than  can  be 
addressed  by  the  group  at  large.  These  are  included  here  for  reference. 
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Model  Initialization  Schemes 


Very  selective  set  of  modes. 

•  Significance  of  large  scales  on  mesoscale  events. 

•  Significance  of  the  short  atmospheric  synoptic  time  scales  that  may 
be  responsible  for  longer  time  scale  ocean  phenomena.  What 
coupling  is  necessary  in  air-sea  interaction  to  generate  El  Nino  for 
example. 

•  Success  and  failure  of  bogus  feature  model  data  in; 

Atmospheres. 

Oceans. 

•  How  better  physics  parameterizations  interact  with  finer  grid 
resolution. 

•  Do  atmospheric  models  contain  enough  correct  information  to  be 
useful  to  the  ocean  and  vice  versa? 

•  Are  any  existing  two-way  coupling  schemes  sufficient? 

•  Is  the  coupling  problem  “just  computational”  or  is  more  theoretical 
work  necessary? 

•  How  much  resolution  is  needed  in  the  different  regions  of  the 
Mediterranean? 

•  How  important  is  model  clim.atology  in  an  area  where  the  circulation 
is  strongly  dependent  on  so  many  interacting  rings? 

•  Significance/impact  of  higher  order  advection  diffusion  differencing. 

•  Boundary  layer  characteristics  need  more  theoretical  descriptions. 

•  Need  to  move  away  from  “it  looks  good”  analysis. 

•  Topographic  vortex  stretching  and  influence  on  mesoscale  events. 
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•  What  is  dominant  time  scale  of  eastern  Mediterranean  for  general 
circulation? 

•  What  are  major  features  of  the  POEM  data  set  that  we  can  utilize? 

•  What  is  the  dominant  forcing  in  each  region  for  each  season? 

•  What  is  the  roll  of  ambient  potential  vorticity  (PV)  in  controlling 
permanent  and  semi-permanent  features? 

•  Need  to  develop  evaluation  tests  that  exercise  different  aspects  of  the 
coupled  systems. 

•  Grid  choices 

Masking  versus  boundary-fitted  grid. 

Multiple  domain  decomposition;  coastal  ocean  not  high 
priority,  but  sub-basin  model  coupling  may  be, 

•  Navy  priority  in  Mediterranean:  Ambient  wave  noise,  surface 
velocity,  and  high  frequencies. 


The  workshop  charter  called  for  the  development  of  a  strategic  plan  that 
outlines  the  significant  research  issues  that  must  be  addressed.  The  following  list 
contains  the  agreed  upon  needed  research. 


Develop  consistent  (air-sea)  parameterization  of: 


•  Deepening  convective  processes,  including  preconditioning,  that 
adequately  models  the  formation  of  the  Levintine  Basin  intermediate 
water  and  the  Ionian  deep  water. 


Quantify  the  separate  atmospheric/ocean  modeling  needs  for  air-sea 
fluxes  including  the  SST  definitions. 


Improve  the  “matching”  of  the  atmospheric  and  ocean  boundary 
layers  with  the  air-sea-wave  layer  physics.  Determine  which  air- 
sea-wave  layer  interaction  processes  dictate  the  coupling  mechanisms 
that  occur  between  the  respective  boundary  layers.  Examine  the 
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utility,  for  coupling,  of  the  available  data  “types.”  Couple  the  0(1) 
SSH  “zo”.  Determine  the  resolution  needed  for  the  respective 
coupled  models.  Improve  our  modeling  capabilities  for  planetary 
boundary  layers. 

•  Collaborate  with  both  experimentalists  and  modelers  to  collect 
appropriate  air-sea  interaction  data  in  the  region,  including 
evaporation  and  precipitation  data  and  better  (denser)  wind  data. 

•  Develop  an  01  assimilation  scheme  for  the  Mediterranean  that 
addresses  the  complex  covariance  structure  of  region  and  the 
complex  physics  that  is  multiscale,  anisotropic,  and  inhomogeneous. 
Then  produce  and  evaluate  the  model  covariance  structure  and 
compare  with  the  analytic  theory  and  with  data. 

•  Examine  ways  to  make  adjoint  methods  computationally  feasible. 

•  Determine  the  effect  of  multiple  scale  interactions,  in  both  space  and 
time,  on  the  dynamics  of  the  region  and  investigate  their  impact  on 
the  numerical  approximations.  The  dominant  scales,  as  well  as  the 
important  intermittent  scales,  need  to  be  identified. 

DEVELOPMENTAL  PLAN 


The  workshop  charter  did  not  stop  at  a  strategic  plan,  but  also  called  for  an 
outline  of  a  developmental  plan.  This  section  outlines  the  proposed  plan. 

•  Survey  existing  atmospheric  modeling 

Identify  archaic  parts. 

Identify  parts  that  need  to  be  consistent  for  (air-sea)  coupling. 
Identify  work  that  should  be  carried  over  to  ocean  models. 

•  Explore  anticipated  hardware  architectures  to  estimate  maximum 
amount  of  horizontal  and  vertical  resolution  that  will  be  available  at 
the  end  of  the  project. 

•  Compare  studies  on  the  Mediterranean  Basin. 
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•  Study  findings  of  non-hydrostatic  research  and,  where  possible, 
incorporate  into  our  “subphysics”  parameterizations. 

•  Conduct  preliminary  studies  on  basin  scales. 

•  Conduct  a  resolution  study  to  determine  sampling  strategies  for; 

Ekman  layer. 

Deep  convecting  chimneys. 

Data  sets. 

•  Reproduce  previous  wind  sensitivity  studies. 

•  Conduct  twin  experiments. 

Test  assimilation  schemes,  tide,  XBT,  SST,  “OTIS”  type 

objective  analysis  (OA)  data. 

Determine  value  added. 

•  Conduct  a  preliminary  study  of  assimilation  using  GEOSAT  1987- 
1988. 

•  Study  numerical  characteristics  (truncation  error,  resolution  needs, 
related  analytic  problem  solved,  cell  Reynolds  instabilities, 
convergence  order,  etc.)  associated  with  various  gridding  strategies, 
(rectilinear,  orthogonal  curvilinear,  etc.). 

a  versus  z  levels. 

Topographic  effects:  numerical  versus  analytical. 

Island  representations  and  their  effects. 

Boundary  fitting  effects. 
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APPENDIX  A 


Individual  Reports 

Jim  Goerss  gave  a  briefing  on  the  existing  Navy  Operational  Global 
Atmospheric  Prediction  System  (NOGAPS).  He  stressed  the  importance  of  data 
quality  in  sparse  data  coverage  regions.  During  his  talk,  discussion  about  the 
relative  advantages  and  disadvantages  of  periodic  model  reinitialization  compared 
to  continuous  data  assimilation  was  lively  (only  one  large  forecast  center  uses  a 
continuous  assimilation  scheme).  The  atmospheric  forecasters  present  suggested 
that  reinitialization,  every  six  hours  for  example,  allowed  better  quality  control 
because  it  rejects  fewer  good  data  points,  and  provides  more  control  of  model 
statistics.  His  write-up  is  on  pages  A-4  to  A-6. 

Richard  Hodur  introduced  the  Navy  Operational  Regional  Atmospheric 
Prediction  System  (NORAPS),  and  the  Coupled  Ocean/Atmosphere  Mesoscale 
Prediction  System  (COAMPS).  He  stressed  that  the  regional  boundary  conditions 
(taken  from  the  global  model,  interpolated  in  space  using  bilinear  interpolation 
every  twelve  hours,  and  smoothed  linearly  over  five  grid  points)  can  seriously 
impact  the  regional  small-scale  noise  over  monthly  time  scales.  Allan  Robinson 
pointed  out  that,  for  assimilation  into  a  local  model,  finer  resolution  is  a  bonus. 
Subsequent  discussion  included  the  observation  that  the  finer  scales  produced 
better  resolved  topography,  and  these  topographic  effects  exacerbate  small-scale 
noise  (the  same  kind  of  small-scale  noise  introduced  by  the  coupled  boundary 
conditions).  His  write-up  is  on  pages  A-7  to  A-8. 

Pat  Phoebus  spoke  on  Westerly  Wind  Bursts,  which  are  zonal,  organized 
winds  that  can  extend  over  20°  in  length.  Her  work  is  in  the  equatorial  Pacific. 
The  hypothesis  was  that  air-sea  coupling  may  work  as  follows:  westerly 
atmospheric  wind  bursts  may  generate  Kelvin  waves  in  the  underlying  ocean 
leading  to  upper  ocean  warming,  possibly  El  Nino.  She  showed  several  examples 
of  westerly  wind  bursts  interacting  with  tropical  storms.  Paola  Rizzoli  asked  if 
counter-  and  contrapositive  examples  to  the  stated  hypothesis  had  been  searched 
for,  but  they  had  not  been.  Pat  suggested  that  atmospheric,  synoptic-scale  wind 
forcing  may  be  necessary  for  ocean  models  to  reproduce  El  Nino.  She  suggested 
that  wind  forcing  on  these  scales  can  be  generated  using  an  atmospheric  model  to 
produce  denser  wind  coverage  in  both  space  and  time  than  is  currently  available 
from  data  products.  Her  write-up  is  on  pages  A-9  to  A- 18. 
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Tom  Rosmond  introduced  an  ongoing  project,  sponsored  by  Lawrence 
Livermore  National  Laboratory,  to  compare  atmospheric  models  forced  with  sea 
surface  temperature  (SST)  data  that  is  updated  every  model  time  step.  The 
models  use  prescribed  initial  conditions.  The  aim  is  to  produce  homogeneous 
data  sets.  He  stressed  the  need  to  “practice  walking  before  we  learn  to  run.”  He 
posed  the  question  “Are  atmospheric  models  good  enough”  (i.e.,  do  they  contain 
sufficient  correct  information)?  He  also  reminded  us  that  currently  all  two-way 
coupled  models  still  use  a  flux  correction  method,  which  Allan  Robinson  agreed 
was  “abhorrent.”  His  write-up  is  on  pages  A- 19  to  A-22, 

Le  Ly  spoke  about  an  ocean  model  coupled  to  a  wave  layer  model  where 
ocean  waves,  the  ocean,  and  the  atmosphere  interact.  The  model  was  developed 
to  use  an  “E-e”  turbulence  scheme.  Based  on  a  scaled  analysis  of  this 
parameterized  “E-e”  turbulent  energy  scheme  he  solves  a  reduced  e  equation 
analytically.  The  major  issue  here  is  matching  the  parameterization  scheme  to  the 
data. 


Allan  Robinson  gave  an  excellent  overview  talk  on  the  Physical 
Oceanography  in  the  Eastern  Mediterranean  (POEM)  project.  Their  results  to 
date  show  that  the  general  circulation  of  the  eastern  Mediterranean  is  a 
multiple-scale  problem  in  both  space  and  time.  The  basin-scale  circulation 
appears  to  be  a  series  of  interacting  mesoscale  events  of  both  quasi-permanent  and 
transient  rings.  They  appear  to  be  the  “building  blocks”  of  the  general 
circulation. 

George  Mellor  presented  runs  of  the  Princeton  primitive  equation  (PE) 
model  in  the  Mediterranean.  He  estimated  that  perhaps  25  to  30  vertical  levels 
will  be  needed.  He  stressed  that  further  convective  studies  to  model  the  episode 
chimneys  are  needed. 

George  Hebum  also  presented  some  model  runs.  He  has  multiple  year 
simulations  of  the  Mediterranean.  He  pointed  out  that  strong  vertical  sheer  in  too 
shallow  a  region,  such  as  occurs  in  straits,  is  highly  baroclinically  unstable. 

Paola  Rizzoli  presented  results  from  her  large-scale  forcing  study.  She 
showed  that  the  impact  of  the  externally  forced  mesoscale  features  on  the  general 
circulation  can  be  modeled  without  sufficient  resolution  to  resolve  the  mesoscale 
features,  provided  you  have  good  external  forcing  in  your  model.  She  also 
discussed  an  adjoint  assimilation  method,  stressing  that  it  was  an  assimilation 
technique  that  was  consistent  with  both  the  model  and  the  data. 


A-2 


Ziv  Sirkes  discussed  his  ideas  for  an  extension  of  the  steady  state  Kalran 
Filter  developed  by  Ichiro  Fukumori,  Jerome  Benveniste,  Carl  Wunsch,  and 
Dale  B.  Haidvogel  that  would  make  the  Kalman  filter  more  cost  effective.  These 
include  separating  the  model  time  step  from  the  assimilation  time  step,  and  using 
the  doubling  algorithm  to  speed  up  calculations. 

Jim  Corbin  reviewed  the  status  of  the  Experimental  Center  for  Mesoscale 
Ocean  Prediction  (ECMOP),  stressing  that  functional  development  support  is 
available. 

David  Dietrich  presented  results  in  the  Gulf  of  Mexico  and  pointed  out  that 
some  of  the  terms  in  the  equations  contribute  dominant  errors  and  should  be 
discretized  more  carefully.  He  also  suggested  the  vertical  resolution  can  seriously 
affect  the  vorticity  evolution.  His  summary  of  the  Plenary  group  discussion  is  on 
pages  A-153  to  A-154. 

Joe  Thompson  has  provided  us  with  an  overview  of  block  structured  grids 
and  general  curvilinear  coordinate  systems.  He  is  the  Director  of  the 
Engineering  Research  Center  at  Mississippi  State  University,  and  is  anxious  to 
assist  individual  modelers  with  grid  problems.  His  write-up  is  on  pages  A-89  to 
A-155. 
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ATMOSPHERIC  DATA  ASSIMILATION 


Jim  Goerss 
NRL  Monterey 


The  major  elements  of  atmospheric  data  assimilation  systems  run  by  the 
world's  forecast  centers  are  quality  control  of  data,  analysis,  initialization,  and  a 
prediction  model.  Data  quality  control  is  vital  to  the  success  of  a  data 
assimilation  system  and  is  normally  performed  prior  to  analysis  and  within  the 
analysis  itself  (where  observations  can  be  readily  compared  with  other 
observations  as  well  as  with  the  analysis  background  or  first-guess).  The  analysis 
combines  information  contained  in  the  observational  data  and  in  the  background 
fields  (a  forecast  made  by  the  prediction  model  valid  at  analysis  time)  to  produce 
a  set  of  initial  fields  for  the  prediction  model.  The  initialization  makes 
adjustments  to  these  fields  so  that  the  prediction  model  can  begin  its  integration 
with  very  little  noise.  In  order  to  have  a  successful  data  assimilation  system,  the 
analysis  must  be  such  that  it  effectively  includes  the  information  contained  in  the 
observations  to  produce  fields  which  are  changed  very  little  by  the  initialization 
step.  Furthermore,  the  prediction  model  must  be  capable  of  explaining  most  of 
the  variance  of  the  changes  in  the  atmosphere  since  the  previous  analysis  so  that 
the  current  analysis  is  only  required  to  make  relatively  small  changes.  For  global 
data  assimilation  systems,  atmospheric  forecast  models  have  reached  the  state 
where  the  error  variance  in  the  analysis  background  fields  is  comparable  to  that 
for  the  observational  data.  However,  the  global  assimilation  problem  is  the 
easiest  one  for  the  atmosphere.  As  the  domain  of  an  atmospheric  data 
assimilation  system  gets  smaller,  the  assimilation  problem  becomes  more 
difficult.  The  global  domain  is  relatively  rich  in  observational  data,  but  small 
mesoscale  domains  can  be  data  poor.  While  atmospheric  forecast  models  are 
quite  skillful  at  predicting  the  features  which  dominate  the  atmosphere  on  the 
global  scale,  a  mesoscale  domain  can  be  dominated  by  small-scale  features  which 
are  not  well  depicted  by  today's  models  (fronts,  convective  complexes,  etc.). 
Furthermore,  as  the  size  of  the  domain  shrinks  so  does  the  length  of  "useful" 
forecasts  by  the  prediction  model.  While  global  models  are  capable  of  making 
skillful  forecasts  of  the  atmosphere  7-10  days  in  the  future,  the  length  of  skillful 
forecasts  for  successful  mesoscale  models  will  be  measured  in  hours  not  days. 

Operational  Global  Data  Assimilation 

The  Navy  Operational  Global  Atmospheric  Prediction  System  (NOGAPS) 
is  the  Navy's  global  data  assimilation  system.  It  consists  of  pre-processing  quality 
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control  of  observational  data  (with  special  emphasis  on  radiosonde  data),  a 
multivariate  optimum  interpolation  (MVOI)  analysis,  non-linear  normal  mode 
initialization,  and  a  spectral  forecast  model.  The  MVOI  analysis  is  performed  at 
16  pressure  levels  from  1000  mb  to  10  mb  upon  the  1.5  degree  resolution 
Gaussian  grid  of  the  forecast  model.  In  order  to  maximize  the  use  of 
observational  data  at  each  analysis  grid  point,  the  globe  is  divided  into  over 
800  volumes.  The  volumes  have  been  chosen  so  that,  given  normal  receipt  of 
observational  data,  about  300  observations  are  available  for  each  volume.  The 
global  analysis  is  produced  by  assembling  the  individual  volume  analyses.  The 
MVOI  analysis  is  a  statistical  analysis  scheme  in  which  the  weights  applied  to  the 
observational  data  are  such  that  analysis  error  is  minimized.  The  weights  are  a 
function  of  the  error  properties  of  the  observational  data  and  the  analysis 
background,  a  6-hour  forecast  from  the  spectral  forecast  model.  In  data  rich 
areas  (like  Northern  Hemisphere  continental  areas)  the  background  errors  and 
observation  errors  are  comparable  in  size,  while  in  data  poor  areas  (like  the 
Southern  Hemisphere)  the  observation  errors  are  smaller  and  more  weight  is  put 
upon  the  observational  data  by  the  analysis.  The  T79  spectral  forecast  model 
(79  waves,  triangular  truncation)  has  18  levels  from  the  surface  to  approximately 
8  mb  and  parameterizations  of  all  major  diabatic  processes  (radiation,  convection, 
shallow  convection,  clouds,  etc.).  At  the  end  of  each  watch  at  FNOC  the  data 
assimilation  cycle  is  continued  by  the  performance  of  a  re-analysis  for  the 
synoptic  time  (OOZ  or  12Z),  a  6-hour  model  forecast,  an  off-time  analysis 
(06Z  or  18Z),  and  a  second  6-hour  model  forecast.  Because  of  the  temperature 
soundings  derived  from  4  polar  orbiting  satellites,  the  global  data  coverage  every 
6  hours  is  good.  Unfortunately,  the  observation  errors  for  these  data  are  double 
those  for  conventional  soundings  (radiosondes).  The  conventional  data  are  still 
the  best  quality  data  received  and  are  primarily  confined  to  the  continental  areas 
of  the  globe.  Therefore,  there  is  a  need  to  supplement  the  database  with  synthetic 
observations  in  certain  areas  of  the  world  (the  Southern  Hemisphere,  the  tropics, 
and  the  North  Pacific  and  Atlantic).  Synthetic  surface  pressure  observations  have 
long  been  produced  for  the  Southern  Hemisphere  and  disseminated  by  the 
Australian  Bureau  of  Meterology.  These  observations  are  assimilated  into 
NOGAPS  along  with  synthetic  "ship"  reports  which  can  be  generated  by  the 
watch  at  FNOC  in  the  vicinity  of  extra-tropical  cyclones  in  the  North  Pacific  and 
Atlantic.  The  impact  of  the  latter  is  to  improve  NOGAPS  analysis  and  short-term 
forecasts  for  these  cyclones.  Because  of  the  forecast  skill  of  the  NOGAPS 
spectral  model  in  the  tropics,  we  have  been  able  to  successfully  assimilate 
synthetic  tropical  wind  soundings  in  the  vicinity  of  tropical  cyclones.  These 
soundings  are  automatically  generated  using  information  contained  in  the 
warnings  issued  by  the  Joint  Typhoon  Warning  Center  (JTWC),  Guam  and  by  the 
National  Hurricane  Center,  Miami.  With  the  assimilation  of  these  soundings  it 
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was  found  that  for  the  1991  typhoon  season  NOGAPS  produced  tropical  cyclone 
track  forecasts  that  were  significantly  better  than  those  produced  by  any  of  the 
track  forecast  aids  used  by  JTWC. 
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MESOSCALE  MODELING  AT  NRL  MONTEREY 


Richard  M.  Hodur 


I.  INTRODUCTION 

Many  types  of  mesoscale  weather  events  can  adversely  affect  military 
operations.  These  events  include  flow  over  and  around  steep  terrain,  land  and  sea 
breezes  and  organized  convection,  such  as  that  found  in  tropical  cyclones  and 
squall  lines.  Even  with  the  recent,  rapid  advances  in  computer  hardware,  global 
models  are  still  unable  to  adequately  resolve  many  of  these  types  of  phenomena. 
This  can  only  be  accomplished  using  mesoscale  models  witli  sufficient  horizontal 
and  vertical  resolutions. 

Two  models  have  been  developed  to  examine  these  forecasting  problems. 
First,  the  Navy  Operational  Regional  Atmospheric  Prediction  System  (NORAPS) 
is  used  operationally  to  produce  3 -dimensional  atmospheric  simulations  over 
limited-areas  with  a  grid  resolution  about  3  times  higher  than  the  Navy  global 
model  (NOGAPS;  Navy  Operational  Global  Atmospheric  Prediction  System). 
Second,  the  Coupled  Ocean/Atmosphere  Mesoscale  Prediction  System  (COAMPS) 
features  a  nonhydrostatic  atmospheric  model  and  a  hydrostatic  ocean  model. 
CO  AMPS  is  used  for  3 -dimensional  atmospheric  and  oceanic  simulations  of 
mesoscale  phenomena  such  as  arctic  leads,  convective  storms,  and  air-sea  fluxes. 

II.  NORAPS 

NORAPS  has  been  used  operationally  at  the  Fleet  Numerical  Oceanography 
Center  since  1982  to  provide  12  to  48  hour  forecasts  for  limited  areas  of  the 
world.  NORAPS  uses  a  globally  relocatable  grid,  variable  resolutions  in  the 
horizontal  and  vertical  and  can  be  used  for  providing  mid-latitude  forecasts  and 
tropical  forecasts,  including  the  path  of  tropical  cyclones.  NORAPS  features  a 
12  h  data  assimilation  cycle,  a  multivariate  optimum  interpolation  analysis, 
nonlinear  vertical  mode  initialization,  a  state-of-the-art  mesoscale  prediction 
model  and  post-processing  software.  NORAPS  is  currently  used  operationally 
over  four  areas  of  the  world:  Mediterranean,  western  Pacific,  western  Adantic, 
and  the  Persian  Gulf,  with  resolutions  ranging  from  40  to  100  km. 

Several  forecasts  have  been  made  recently  that  demonstrate  the  usefulness 
of  the  high  resolution  of  NORAPS  for  mesoscale  events.  Using  80  km  resolution, 
NORAPS  was  able  to  correctly  forecast  a  mistral  of  approximately  50  knots 
while  NOGAPS,  using  a  resolution  of  about  150  km  forecast  maximum  wind 


speeds  of  only  20-25  knots.  In  another  case  of  55  knot  winds  near  Sicily,  a 
40  km  version  of  NORAPS  forecast  maximum  winds  of  45  knots,  an  80  km 
version  of  NORAPS  forecast  maximum  winds  of  35  knots,  and  a  version  of 
NOGAPS  using  approximately  125  km  resolution  forecast  maximum  winds  of 
about  30  knots.  These  forecasts  indicate  that  NORAPS  can  serve  as  a  model 
which  is  capable  of  predicting  mesoscale  events,  particularly  in  areas  such  as  the 
Mediterranean,  where  small-scale  forcing  from  terrain  and  detailed  coastlines  are 
important. 

m.  COAMPS 

COAMPS  is  being  developed  as  a  community  R&D  model  for  mesoscale 
studies.  COAMPS  features  a  nonhydrostatic,  compressible  atmospheric  model 
coupled  to  a  hydrostatic,  incompressible  ocean  model.  The  atmospheric  model 
uses  1-1/2  order  turbulence  closure  and  explicit  moist  physics,  while  the  ocean 
model  uses  first-order  turbulence  closure  and  a  free  surface.  COAMPS  is 
appropriate  for  studying  mesoscale  phenomena  such  as  convection,  tropical 
cyclones,  arctic  leads,  and  coastal  and  mountain  flows  with  horizontal  resolutions 
ranging  from  20  km  down  to  1  m. 

COAMPS  has  been  used  to  simulate  hurrirane  Gilbert  (September  1988) 
which  attained  a  minimum  central  prt;ssure  of  approximately  888  mb  with 
maximum  wind  speeds  of  approximately  160  kts.  Using  only  the  atmospheric 
model,  with  idealized  initial  conditions,  COAMPS  fo'^ecast  the  storm  to  intensify 
from  1000  mb  to  890  mb  in  66  hours  with  peak  winds  of  163  kts.  When  the 
ocean  model  is  used  (i.e.,  the  sea  surface  temperature  is  allowed  to  vary  in  time), 
the  model  storm  intensifies  at  a  slower  rate  due  to  the  lower  sea  surface 
temperature  caused  by  the  upwelling  created  by  the  storm  circulation.  Presently, 
the  coupled  experiment  has  only  been  run  to  48  h,  and  the  forecast  intensity  is 
about  10  mb  higher  than  the  fixed  sea  surface  temperature  case  at  the  same  time. 

IV.  CONCLUSION 

NRL-Monterey  is  aggressively  pursuing  the  development  of  mesoscale 
modeling  for  optimum  use  by  the  Navy.  NORAPS  has  been  developed  and  is 
being  continually  improved  to  serve  as  the  central  site  mesoscale  model  for  12  to 
48  h  guidance  in  coastal  areas.  COAMPS  is  being  developed  as  a  community 
research  model  to  study  mesoscale  phenomena  that  cannot  be  resolved  by 
NORAPS  for  the  goals  of  better  understanding  these  phenomena,  as  well  as  for 
improving  physical  parameterizations  for  the  larger-scale  models. 
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COUPLED  ATMOSPHERE/OCEAN  MODELING  IN  THE  TROPICAL  PACIFIC 

Patricia  A.  Phoebus 
Naval  Research  Laboratory,  Monterey 


1.  Introduction 

Keen  (1982)  suggested  that  ocean  Kelvin  waves  moving  from  west  to  east  in  the  equatorial  Pacific 
waveguide  may  be  at  least  partially  responsible  for  the  appearance  of  extremely  warm  water  in  the 
central  and  eastern  Pacific,  a  phenomena  that  is  commonly  referred  to  as  El  Nino.  These  waves 
may  be  triggered  by  short  episodes  of  westerly  wind  anomalies  near  the  equator  in  the  western 
part  of  the  basin.  These  westerly  wind  anomalies  are  particularly  strong  when  there  is  a  tropical 
cyclone  couplet  present  in  the  western  Pacific.  Because  this  area  is  devoid  of  radiosonde  stations 
(Fig.  1),  previous  studies  to  detect  westerly  wind  events  have  concentrated  on  analyzing  data  from 
island  stations  (Harrison  and  Giese,  1991)  and  from  shipping  lane  traffic  (Morrissey,  1990).  For  the 
same  reason,  attempts  to  model  the  ocean’s  response  to  this  anomalous  forcing  have  relied  upon 
specified  wind  patches  rather  than  actual  observed  wind  anomalies  (McPhaden  et  al.,  1987). 

This  study  proposes  that  a  four-dimensional  atmospheric  data  assimilation  system  can  provide  the 
best  estimate  of  the  structure  of  the  atmosphere  in  data  void  areas.  This  statement  is  particularly 
valid  in  areas  such  as  the  western  equatorial  Pacific,  where  there  are  ample  observations  available 
surrounding  the  area  of  interest.  A  data  cissimilation  system  can  use  these  observations  to  produce 
three-dimensional  wind  and  mass  fields  that  are  consistent  with  known  dynamical  and  physical 
properties  of  the  atmosphere.  While  modeling  the  tropical  atmosphere  presents  a  great  challenge, 
the  Navy  Operational  Global  Atmospheric  Prediction  System  (NOGAPS)  has  performed  admirably 
in  the  western  Pacific  since  mid- 1990  (Goerss  et  aJ.,  1991).  Thus,  the  NOGAPS  low-level  wind  fields 
may  provide  additional  information  about  the  temporal  and  spatial  scales  of  westerly  wind  events 
in  the  equatorial  Pacific. 

NOGAPS  contains  four  components  -  data  quality  control,  data  analysis,  model  initialization,  and 
model  foreciist  (see  Goerss,  this  report)  -  and  is  run  with  a  6-hr  update  cycle.  That  is,  the  6- 
hr  model  forecast  is  used  as  the  first  guess  for  the  next  analysis,  which  in  turn  uses  the  new'est 
observations  to  make  incremental  change?  to  the  first-guess  fields.  Conventional  observations  are 
available  from  radiosondes,  pilot  balloons,  aircraft,  ships,  and  buoys.  Remotely  sensed  observations 
include  temperature  soundings  from  the  polar-orbiting  satellites,  cloud-tracked  wind  vectors  from 
the  geostationary  satellites,  and  SSM/I  wind  speed  observations  from  the  DMSP  polar  orbiters  (a 
data  source  uniquely  used  by  the  Navy).  NOGAPS  also  makes  use  of  synthetically  generated  ob¬ 
servations  to  better  depict  the  strength  of  oceanic  extratropical  and  tropical  cyclones  (Phoebus  and 
Goerss.  1992).  NOGAPS  is  provided  a  real-time  sea  surface  temperature  field  from  an  operational 
ocean  data  assimilation  system  that  includes  a  three-dimensional  ocean  thermal  analysis  coupled 
to  a  mixed-layer  forecast  model  (Clancy  et  al.,  1991)  that  is  forced  with  NOGAPS  wind  stresses 
and  surface  fluxes. 

Since  the  tropical  ocean  responds  very  rapidly  to  large-scale  changes  in  the  surface  winds,  the  Naval 
Research  Laboratory  (NRL)  is  also  investigating  the  impact  of  modeled  westerly  wind  events  on  an 
ocean  circulation  model  (John  Kindle.  NRL,  Mississippi).  While  the  model  is  not  operational  at  this 
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time,  it  is  being  run  in  research  mode  with  strictly  a  one-way  coupling.  That  is.  there  is  no  feedback 
from  the  ocean  circulation  model  to  any  other  component  of  the  operational  ocean/atmosphere  data 
assimilation  systems.  Most  of  the  e.xperiments  are  being  conducted  with  a  1.5  layer  global  model 
that  includes  highly  parameterized  thermodynamics.  The  model  resolution  is  0.5°  latitude  by  0.7° 
longitude,  and  it  is  forced  with  the  12-hourly  NOGAPS  surface  wind  stresses.  By  studying  the 
ocean  model's  response  to  the  westerly  wind  forcing,  we  hope  to  distinguish  which  wind  events  are 
significant  and  what  impact  ocean  Kelvin  waves  have  on  sea  level  changes  in  the  eastern  equatorial 
Pacific. 

2.  Results 

Figure  2  shows  a  times  series  of  the  NOGAPS  10-m  zonal  wind  stress  component  compared  to  the 
zonal  wind  stress  computed  from  moored  buoy  observations  at  165E,  2N.  Westerly  wind  events  are 
noted  by  a  strong  positive  peak  in  wind  speed  (light  wind  speeds  are  indicative  of  either  variable  or 
more  meridional  wind  directions  and  are  not  of  interest  here).  Overall,  there  is  very  good  agreement 
between  the  model  and  the  observations,  even  though  data  from  this  buoy  was  only  sporadically 
available  to  NOGAPS.  The  strongest  westerlies  were  observed  in  mid-to-late  November  1990.  and 
were  associated  with  developing  storms  Page  and  Owen,  both  of  which  became  super-typhoons 
(Fig.  3).  But,  even  though  the  westerlies  are  very  strong  at  2N,  there  is  a  much  less  organized  flow 
along  the  equator,  so  the  ocean  response  might  be  somewhat  less  than  expected. 

The  peak  in  December  1990,  while  weaker  at  2N  than  the  previous  case,  was  actuaJly  stronger 
along  the  equator  (Fig.  4).  This  westerly  event  was  associated  with  coupled  northern  (Russ)  and 
southern  hemisphere  tropica  storms  (Fig.  5),  much  like  those  described  by  Keen  (1982).  However, 
where  Keen  relied  on  the  relative  positions  and  timing  of  the  storms  to  define  a  couplet,  by  using 
the  model  output,  we  can  clearly  see  that  the  circulations  defining  these  two  storms  axe  dynamically 
coupled.  Another  couplet  (Kelvin/Sharon)  developed  in  early  March  1991.  Although  strong  zonal 
winds  were  present  along  the  equator  (Fig.  6),  there  was  little  evidence  of  the  event  at  165E,  2N 
(Fig.  2).  These  situations  clearly  illustrate  the  potential  hazards  involved  with  using  isolated  ship 
and  island  observations  to  draw  conclusions  about  the  strengths  and  locations  of  westerly  wind 
events. 

Without  showing  figures  for  each  of  these  events,  I  will  note  that  the  westerly  peak  in  mid-March 
1991  was  associated  with  Typhoon  Tim,  which  developed  very  near  the  equator.  In  early  May 
1991.  another  couplet  formed  (Walt/Lisa)  around  155E,  concentrating  the  westerly  winds  between 
them.  .4  monsoon  depression  formed  along  the  monsoon  trough  in  June  1991,  while  July  1991  had 
a  series  of  tropical  storms  in  the  northern  hemisphere  (Amy,  Brendan,  Caitlin).  The  westerlies  in 
early  September  1991  were  associated  with  Typhoon  Ivy,  and  in  late  September  the  early  signatures 
of  Orchid  and  Pat  were  evident  near  the  equator,  with  westerlies  along  the  southern  boundaries 
of  these  circulations.  The  same  type  of  pattern  was  evident  throughout  the  remainder  of  1991, 
with  major  westerly  wind  events  appearing  in  conjunction  with  tropical  storms  developing  near  the 
equator  in  one  or  both  hemispheres. 

Finally.  Figure  7  shows  the  observed  sea  level  variation  at  Baltra  (90W.  0.5S),  one  of  the  Galapagos 
Islands.  Given  that  it  takes  around  two  months  for  an  ocean  Kelvin  wave  to  transverse  the  equa- 
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torial  Pacific,  the  rise  in  sea  level  in  February  1991  is  likely  due  to  an  ocean  Kelvin  wave  generated 
by  the  strong  westerlies  along  the  equator  in  December  1990.  Similarly,  the  rises  beginning  in  late 
April  on  into  May  would  be  expected  if  the  early  March  (Sharon/Kelvin)  westerly  event  triggered 
a  Kelvin  wave,  followed  shortly  thereafter  by  another  wave  generated  by  the  westerlies  associated 
with  Typhoon  Tim  in  mid-March.  Similar  results  have  been  seen  in  the  ocean  circulation  model, 
which  lends  further  support  to  this  theory. 

3.  Summary 

In  almost  every  case  of  westerly  winds  along  the  equator,  there  has  been  a  tropical  storm  developing, 
which  lends  further  support  to  Keen's  original  observation,  although  not  every  case  involves  a 
couplet.  Furthermore,  in  most  cases,  the  westerlies  appear  in  the  very  early  stages  of  the  storm's 
development,  usually  before  it  reaches  tropical  storm  strength.  Thus,  there  is  some  question  as 
to  whether  the  westerlies  are  causing  the  storm’s  development,  or  vice  versa.  These  types  of 
cause/effect  questions  are  at  the  heart  of  the  El  Nino/Southern  Oscillation  question,  and  there 
is  still  much  to  be  learned  about  what  role  is  played  by  each  of  the  pieces.  While  observed  and 
modeled  changes  in  sea  level  due  to  individual  ocean  Kelvin  waves  are  not  large  enough  to  fully 
explain  the  El  Nino  phenomena,  it  is  reasonable  that  an  appropriately  timed  series  of  such  events 
could  result  in  larger  overall  changes  in  the  eastern  Pacific.  Global  data  assimilation  systems  and 
coupled  ocean/atmosphere  models  may  eventually  be  able  to  provide  some  of  the  answers. 
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Fig.  2.  Time  series  of  NOGAPS  10-m  zonal  wind  stress  (solid)  compared  to  buoy  observed 
zonal  wind  stress  (dashed)  at  165?  E,  2?N.  Units  are  Pascals,  and  a  5-day  filter  has  been 
applied  to  the  12-hourly  data. 
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Fig.  3.  NOGAPS  WOO  mb  streamline/wind  analysis,  OOZ  21  November 
1990.  The  early  circulations  associated  with  Supertyphoon  Page  (1423  E) 
and  Supertyphoon  Owen  (16TE)  are  clearly  developed  by  the  model  before 
any  synthetic  observations  are  inserted. 
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Fig.  6.  NOGAPS  1000  mb  wind  analy.'fis,  OOZ  5  March  1991.  Wind.^  along 
the  equator  between  Tropical  Storm  Kelvin  and  tropical  depression  Sharon 
are  almost  due  westerly  at  25  kt. 


A-17 


NOGAPS  DEVELOPEMENT  IN  SUPPORT  OF  COUPLED 
ATMOSPHERE/OCEAN  MODELING 


Tom  Rosmond 
NRL-Monterey 


1.  INTRODUCTION 

High  priority  Navy  ocean  model  research  has  stimulated  the  development 
effort  of  the  Navy  Operational  Global  Atmospheric  Prediction  System  (NOGAPS) 
in  support  of  this  effort.  The  most  significant  areas  of  NOGAPS  research  are 

(1)  reduction  of  model  systematic  error,  particularly  of  the  surface  fluxes,  and 

(2)  development  of  a  fully  coupled  NOGAPS/TOPS  prediction  system.  The 
Thermodynamic  Ocean  Prediction  System  (TOPS)  is  the  Navy’s  ocean  mixed  layer 
model  run  operationally  by  Fleet  Numerical  Oceanography  Center  (FNOC). 

Systematic  or  bias  error  reduction  is  an  ongoing  priority  for  NOGAPS 
development  time  because  of  its  importance  not  only  for  the  coupled  model 
problem,  but  also  because  it  is  a  major  obstacle  to  successful  medium  range 
forecasts  and  climate  related  model  research.  A  useful  tool  for  diagnosing  bias  is 
extended  integration  with  prescribed  boundary  conditions.  This  report  describes  a 
current  project  of  this  type,  NOGAPS  participation  in  a  DOE  sponsored 
Atmosphere  Model  Intercomparison  Program  (AMIP). 

A  fully  coupled  NOGAPS/TOPS  global  forecast  system  has  been  under 
development  at  NRL  Monterey  for  some  time.  A  practical  problem  with  any 
coupled  atmosphere/ocean  model  is  the  mismatch  of  time  and  space  scales  between 
the  atmosphere  and  ocean.  The  Navy's  interests  in  ocean  modeling  ultimately 
require  hi^  resolution  both  in  the  vertical  and  horizontal  to  resolve  oceanic  eddys 
and  fronts.  The  computational  requirements  of  this  are  formidable,  even  exceeding 
the  capabilties  of  class  VTI  computers.  In  this  report  a  coupled  NOGAPS/TOPS 
system  with  TOPS  covering  only  the  Mediterranean  and  Black  Seas  is  described. 
This  system  is  a  prototype  for  other  more  elaborate  coupled  atmosphere/ocean 
systems  that  could  satisfy  Navy  needs  without  prohibitive  computational  cost. 

2.  NOGAPS  AMIP  PARTICIPATION 

Recent  national  interest  in  global  climate  change  has  stimulated  extra  effort 
in  the  global  modeling  commmunity  to  improve  the  models.  The  most  glaring 
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deficiency  in  all  global  general  jirculation  models  (GCM’s)  is  the  systematic  error 
or  bias,  which  makes  application  of  these  models  to  the  problems  of  climate 
change  or  even  extended  range  forecasting  questionable.  Nowhere  is  bias  error 
more  critical  than  in  interactive  atmosphere/ocean  model  systems,  where  the 
prediction  or  simulation  of  realistic  sea  surface  temperature  (SST)  is  central  to  the 
entire  concept  of  coupled  model  research.  To  date  efforts  to  run  fully  interactive 
coupled  systems  without  any  constraints  on  the  predicted  SST  show  unacceptable 
'climate  drift’.  Clearly  we  must  improve  the  GCM’s  before  practical  coupled 
atmosphere/ocean  models  that  satisfy  Navy  requirements  are  possible. 

The  DOE  AMIP  experiment  is  an  ambitious  effort  to  address  these  common 
deficiencies  of  all  GCM's.  The  Lawrence  Livermore  National  Laboratory  (LLNL) 
has  provided  a  10  year  history  (1979-1988)  of  SST  and  sea  ice  boundary 
conditions  to  all  interested  participants.  Over  30  global  modeling  groups  from 
around  the  world  are  running  their  models  over  this  period.  The  list  is  a  who’s  who 
of  these  modeling  groups.  NOGAPS  is  a  natural  participant  because  the  Navy  has 
much  to  gain  from  the  knowledge  gained  not  only  from  NOGAPS  results,  but  also 
the  results  from  all  the  other  models  as  well.  LLNL  will  archive  the  models 
history  files,  including  special  diagnostics  used  to  compare  the  models  systematic 
errors. 


NRL  Monterey  has  completed  the  first  5  years  of  the  AMIP  integration 
(1979-1983).  Note  that  this  period  includes  the  1982-83  El  Nino  episode. 
Preliminary  examination  of  NOGAPS  results  show  dramatically  different 
interannual  patterns  of  surface  fluxes,  cloud  cover,  and  other  important 
atmospheric  parameters  over  the  entire  Pacific  basin.  After  completing  the  full 
10  years  of  the  control  integration,  selective  rerunning  of  interesting  periods  is 
planned.  Of  particular  interest  will  be  running  the  coupled  NOGAPS/TOPS  to  see 
how  well  the  observed  SST  patterns  can  be  reproduced  with  the  predictive  SST. 
Careful  analysis  and  comparison  of  surface  energy  budgets  between  the  coupled 
and  fixed  SST  cases  will  give  insight  into  the  reasons  for  NOGAPS  systematic 
errors.  The  NOGAPS  AMIP  data  base  will  be  documented  and  made  available  to 
all  interested  parties.  It  will  provide  a  unique  data  set  of  long  term,  homogeneous 
atmospheric  fields  which  have  relatively  high  resolution  in  both  time  and  space. 
Ocean  modelers  in  particular  should  benefit  from  this  data  set,  both  because  their 
own  modeling  needs  are  met  and  also  because  they  can  provide  valuable  feedback 
to  NOGAPS  model  development. 

Anything  we  can  do  to  improve  NOGAPS  will  benefit  the  Navy  in 
many  important  ways.  Reduced  NOGAPS  systematic  error  will  have  benefits 
from  better  short  range  tactical  weather  forecast  to  better  climate 
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change  simulations,  all  of  which  are  in  the  Navy’s  operational  and  research 
interests. 

3.  COUPLED  NOGAPS/REGIONAL  TOPS 

The  full  global  NOGAPS/TOPS  is  a  useful  research  tool  for  studying  many 
aspects  of  atmosphere/ocean  interaction.  However,  the  global  system  has 
horizontal  resolutions  limited  by  the  computational  costs  of  running  an 
atmospheric  GCM.  Typically  this  resolution  is  100  to  200  kilometers  both  in 
NOGAPS  and  TOPS.  For  modeling  the  direct  response  of  the  ocean  mixed  layer 
to  atmospheric  forcing,  this  resolution  is  adequate,  but  in  areas  where  ocean  fronts 
and  eddies  are  present  it  is  far  too  coarse.  Typical  horizontal  scales  for  oceanic 
fronts  and  eddies  are  lO’s  of  kilometers,  to  TOPS  must  have  horizontal  grid 
spacing  fine  enough  to  resolve  these  features. 

The  Mediterranean  Sea  provides  an  excellent  test  bed  for  a  regional 
NOGAPS/TOPS  model.  Complex  small  scale  circulation  patterns  dominate  the 
basin,  and  intense  atmospheric  forcing  is  typical  in  many  of  the  important  observed 
features.  The  NOGAPS/TOPS  model  was  modified  to  include  only  a  regional 
TOPS  which  covered  the  Mediterranean  and  Black  Seas.  The  horizontal  grid  of 
the  regional  TOPS  was  about  20  kilometers,  1/10  the  NOGAPS  resolution.  Every 
^e  step  NOGAPS  surface  fluxes  are  interpolated  to  the  TOPS  fine  grid  and  the 
fine  grid  TOPS  predicted  SST  is  averaged  over  the  grid  area  corresponding  to  the 
coarse  NOGAPS  grid.  The  model  has  been  integrated  for  20  days  for  a 
January  case  with  no  apparent  problems.  Predicted  TOPS  SSTs  show  realistic 
patterns  when  driven  by  the  NOGAPS  surface  fluxes.  This  is  in  spite  of  the  fact 
Aat  the  initial  SST  distribution  was  derived  from  a  coarse  mesh  climatological 
SST.  There  are  two  serious  deficiencies  still  to  be  overcome  for  a  successful 
coupled  regional  NOGAPS/TOPS,  however.  First,  no  thermohaline  circulation 
currents  were  available  to  add  to  the  TOPS  generated  Ekman  currents.  These  will 
have  to  eventually  come  from  a  3-dimensional  circulation  model  for  the 
Mediterranean  and  Black  Seas.  Second,  the  thermodynamic  fluxes  produced  by 
NOGAPS  are  not  sensitive  to  the  high  resolution  SST  from  TOPS,  since  to 
OGAPS  the  details  of  the  TOPS  SST  are  sub-grid  scale.  A  possible  solution  is  as 
follows:  (1)  interpolate  the  NOGAPS  fluxes  to  the  TOPS  fine  mesh.  (2)  Assume 
dwt  air/sea  temperature  difference  variations  over  the  TOPS  grid  are  dominated  by 
SST  variations,  i.e.,  the  NOGAPS  predicted  surface  air  temperature  is  constant 
over  a  NOGAPS  grid  square.  (3)  Define  a  linear  parametric  model  for  the  variation 
0  a  reference  surface  flux  as  a  function  of  the  variation  from  a  reference  air/sea 
difference.  In  this  case  the  references  are  the  fluxes  and  temperature 
difference  at  the  NOGAPS  coarse  mesh  grid  point.  (4)  Use  the  constraint  that  the 
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area  integrared  surface  fluxes  over  the  TOPS  fine  mesh  points  must  equal  the  area 
integrated  fluxes  for  the  NOGAPS  grid  point  that  corresponds  to  the  TOPS  points 
to  determine  coefficients  and  parameters  of  the  surface  flux  variation  model.  Other 
constraints  may  also  be  possible  and  desirable,  but  the  approach  will  be  similar. 
The  benefit  of  such  an  adjustment  model  will  be  energetically  consistent  high 
resolution  surface  fluxes  that  will  provide  realistic  high  resolution  air/sea  feedback 
between  the  atmosphere  and  ocean  mixed  layer. 

The  principal  appeal  of  the  coupled  regional  NOGAPS/TOPS  is  the  ability 
to  realistically  model  high  vertical  and  horizontal  resolution  in  the  mixed  layer 
without  excessive  computational  cost.  The  Mediterranean  and  Black  Sea  system  is 
especially  attractive  because  of  its  closed  nature,  only  the  Strait  of  Gibralter 
requiring  inflow/outflow  information  which  is  relatively  well  known.  This  is  in 
contrast  to  open  ocean  regional  domains  where  lateral  boundary  information  is 
poorly  measured  and  understood.  Circulation  models  of  the  Mediterranean 
combined  with  a  regional  NOGAPS/TOPS  can  be  an  important  tool  for  studying 
general  problems  of  atmosphere/ocean  interaction  which  will  have  general 
application  to  other  coupled  model  applications. 
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1  Introduction 


The  study  of  air-sea-wave  interaction  is  one  of  the  most  important  problems  of  physics  of 
the  atmosphere  and  ocean.  This  problem  is  very  important  not  only  in  the  theoretical  study 
of  the  mechanism  of  energy  transfer  in  an  air-sea-wave  system,  but  also  in  parameterization 
of  the  atmospheric  and  oceanic  boundary- layer  processes  in  ocean  circulation  and  the  at¬ 
mospheric  climate  modeling.  The  transport  of  momentum,  heat,  humidity  and  salt  always 
occurs  at  the  air-sea  interface  and  is  affected  by  the  ocean  surface  waves.  The  character  of  the 
transport  is  regulated  by  the  atmospheric  and  oceanic  turbulent  boundary  layers  (.AOBL). 
The  influence  of  the  ocean  surface  waves  on  the  atmospheric  and  oceanic  surface  layers  has 
been  well  noted  by  Geernaert  et  al.  (1988),  Dubov  (1974),  Kitaigorodski  (1973),  Lai  and 
Shemdim  (1971),  Roll  (1965),  and  others.  This  process  is  very  complex  and  still  poorly 
understood. 

The  difficult  problem  in  modeling  air-sea  interaction  is  how  to  describe  the  ocean  surface 
wave  layer  (order  of  wave  height).  We  don’t  have  sufficient  data  to  specify  influences  of  ocean 
surface  waves  on  AOBL  structure,  and  don’t  know  yet  how  to  write  hydrodynamic  equations 
for  this  wave  layer  (WL).  The  WL  effects  are  expressed  through  a  roughness  length,  zq,  in 
traditional  physical  oceanography.  The  WL  should  not  be  considered  as  a  plane  surface. 
Discontinuity  of  various  physical  characteristics  are  expected  across  the  air-sea  interface. 

A  coupled  model  is  developed  to  study  surface  wave  effects  and  variation  of  the  atmo¬ 
spheric  stratification  on  structure  of  AOBL.  This  is  a  part  of  developing  an  air-sea-wave 
interaction  model.  Both  E-i  [i  is  length  scale,  E  is  turbulence  kinetic  energy)  and  E-s 
turbulence  schemes  are  used  in  this  coupled  model. 
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2  The  Model 


The  model  consists  of  a  set  of  equations  for:  (1)  momentum.  (2)  turbulent  kinetic  energy 
■TKEi.  i31  length  scale  [E-(.  turbulence  scheme).  (4)  energy-dissipation  li:).  (o)  turbulent 
e.xchange  coefficient  (TEC)  expressed  in  terms  of  TKE.  £  ( E-i  turbulence  scheme),  and  c  ( E- 
z  turbulence  scheme),  and  (6)  stratifications  in  the  atmosphere  and  ocean.  These  equations 
are  written  in  the  same  form  for  both  the  atmospheric  and  oceanic  boundary  layers  (  Ly.  1986. 
1990.  1991 ).  Vertical  profiles  of  air  temperature  and  seawater  density  are  approximated  with 
functions  to  represent  observations.  Stratification  in  the  atmosphere  is  parameterized 
by  the  nondimensional  surface  heat  flux.  The  model  includes  equations  for  VVL  describing 
the  transport  of  mean  energy  and  TKE  fluxes  across  the  air-sea  interface.  The  equation  for 
turbulent  kinetic  energy  flux  can  be  written  (Ly,  1990) 


20a 


+  (/0(P,^)  U.=  C,^,u;. 


(1) 


The  equation  for  the  mean  energy  fluxes  are  written  (Ly,  1992): 

/  \  ,  I  *7 

(KaPa-^Ua)  Ua  +  |^.=  CiPaul„ 


{KaPa^Va)  Uo,  +  {K,p,^V,)  |^,= 


oy 


(2) 


(3) 


where  C\  and  a  companion  Cj  are  nondimensional  parameters  representing  the  surface  wind 
wave  effect.  Subscript  a  refers  to  air  and  s  to  seawater.  The  roughness  lengths  at  the  upper 
and  lower  boundaries  of  the  wave  layer  are  roo  and  zq,.,  the  air  and  seawater  densities  are 
and  p,,  and  the  velocity  components  of  both  fluids  are  u  and  v. 
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3  An  E—:  Turbulence  Scheme  And  Ocean  Breaking  Wave  Effects 


The  theory  of  length  scale  has  had  great  success  in  modeling  geophysical  flows  over  fiat 
terrain.  It  is  not  always  obvious  how  to  apply  a  length  scale  to  flows  over  complex  terrain 
(ocean  surface  waves)  and  in  the  middle  and  upper  (lower)  parts  of  AOBL.  The  E-e  closure 
method  makes  use  of  an  £-  equation  instead  of  a  length  scale  equation.  A  few  applications 
of  the  E-£  method  have  been  made  for  the  atmospheric  boundary  layer  (Lee  and  Kao.  1979: 
Detering  and  Etling,  1985:  Stubley  and  Rooney,  1986),  the  ocean  boundary  layer  (Marchuck 
et  al..  1977:  Svensson,  1979),  and  both  AOBL  (Ly,  1991).  The  advantage  of  the  E-~  method 
is  that  it  is  a  new  tool  for  studying  turbulence  structure,  especially  since  e  can  be  measured 
directly  in  the  ocean  (Monin  and  Ozmidov,  1985).  It  is  easier  to  include  the  breaking  wave 
effect  into  a  model.  The  disadvantage  of  this  method  is  that  the  e-equation  introduces  a  new 
set  of  constants.  The  E-s  turbulence  scheme  includes  TKE,  e,  and  Kolmogorov  equations. 


The  TKE  for  the  atmosphere  and  ocean  is  determined  from  the  budget  equation 
dEi  ®3.-  g  dSi,  ,  d  ^^dEi  Ej 


Soi  dzi 


'dzi  '  dZi 


where  subscript  i  =  a  for  the  atmosphere  and  i  =  s  for  the  oce2ui;  Soa  and  5a  are  the  mean 
and  local  values  of  potential  temperature,  9a  or  of  seawater  density;  5o»  and  5,  are  the  mean 
and  local  values  of  density,  /?,. 

Energy-dissipation,  £j,  for  both  the  atmosphere  and  ocean  is  determined  from  the  follow¬ 
ing  equation. 


dei 

~di 


0U  5o.  dz,  ^  dzi  dzi  £, 


(5) 


The  Kolmogorov  equation  in  which  the  TEC  is  expressed  in  terms  of  the  TKE  and 


energy-dissipation  is 


K.  =  Ofc.  £,V£. 


(6) 
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where  a^,  =  0.046  determined  from  experimental  data  in  the  surface  layer  for  both  the 
atmospheric  and  oceanic  equations.  The  constants  q,.  i,  and  others  constants,  and  more 
details  on  equations  (4)-(6)  can  be  found  in  Ly  (1991). 


3.1  Boundary  Conditions 

3.1.1  At  the  air-sea  interface 


The  surface  boundary  conditions  for  the  turbulent  kinetic  energy  can  be  written  in  the 
traditional  form  for  the  atmospheric  boundary  layer  over  land  (Monin  and  Yaglom,  1971) 
and  for  the  air-sea  boundary  layer  problem  (Ly,  1986),  as  follows 


E.izi)  U. 


(7) 


where  u.^  is  the  friction  velocity  of  the  air  cind  water;  and  zoi  is  the  roughness  length  for 
both  air  (i  =  a)  and  water  {i  =  s). 

By  ignoring  the  buoytincy  terms,  Ly  (1991)  solved  the  e-equation  (5)  analytically  at  the 
air-sea  interface  and  obtained  boundary  conditions  for  e  at  surface  in  the  form: 


u: 


ea(^a)  u.=  7r^(9ia  +  (l  “  <Zio)  exp  (-(720)] 


for  the  atmosphere  (i  =  a)  where 


—  _  /  ^2a  ^  I  /2 

1  7  ?2o  (  too  ) 


<7la 


4a 


and 


(8) 

(9) 


eAz,)  Uo,=  +  (1  -  qi,)exp{-q2,)] 


(10) 
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tor  the  ocean  i  i  =  >  l  where 

J  1/2  3 

■^li  _  ,  '^25  ^  1/2  Mil 

3.1.2  Ocean  Breaking  Wave  Effect 

Ocean  wave  breaking  is  another  important  source  of  turbulent  energy  dissipation  besides  the 
dissipation  of  shear  production.  The  ocean  turbulent  energy  dissipation  of  shear  production 
at  the  air-sea  interface  is  shown  in  (10)  and  (11).  In  the  case  where  the  breaking  wave  effect 
is  taken  into  account,  the  turbulent  energy  dissipation  e  at  the  ocean  surface  has  the  form: 

U.=  +  -<7x*)  exp  (-<723)1+ £ou,  (12) 

where  and  are  the  same  as  in  (11),  and  £o^  is  the  turbulent  energy  dissipation  of  the 
wave  breaking  at  the  air-sea  interface.  Following  Benilov  (1973),  and  Monin  and  Ozmidov 
(1986),  the  energy  dissipation  £ow  can  be  expressed  in  term  of  ocean  surface  wave  elements: 


Vph  (h 


A  VA 


where  h,  and  A  are  the  wave  phase  velocity,  height,  and  length,  respectively,  and  7  is 
a  dimensionless  constant.  It  is  noted  that  in  the  case  of  no  surface  waves,  equation  (12) 
reduces  to  equation  ( 1 0 )  and  we  have  a  situation  described  by  Ly  (1991).  This  is  the  first 
attempt  ever  to  model  breaking  wave  effect. 


3.1.3  Away  from  the  interface 

■Away  from  the  air-sea  interface  the  turbulent  kinetic  energy  and  dissipation  (c)  tend  toward 
zero.  The  procedure  and  algorithm  for  solving  the  E-e  problem  are  similar  to  the  one  for 
the  E-f.  model  (see  Ly.  1989). 
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4  Simulation  and  Comparison 


The  basic  coupled  model  has  been  developed  to  study;  il)  the  barociinic  effects  in  the 
atmosphere  and  ocean  on  boundary  layer  structures.  (2)  the  surface  wave  effects  on  the 
upper  turbulent  layer  in  the  ocean.  (.3)  E-e  turbulence  scheme  in  a  coupled  model,  and  (4) 
the  effects  of  the  angle  between  wind  stress  and  wind  velocity  vectors  on  the  aerodynamic 
drag  coefficient.  Several  numerical  experiments  were  carried  out  to  study  the  ocean  surface 
etfect  and  E-e  turbulence  scheme  in  a  coupled  model.  The  numerical  results  are  compared 
with  available  data.  The  results  of  previous  studies  (Ly,  19S6.  1990.  1991)  are  summarized 
here.  .An  example  of  the  model  output  on  the  neutral  drag  coefficient,  Co,  in  comparison 
with  published  drag  coefficients  is  shown  in  Fig.  1.  Curve  10  is  the  model  Co  for  Ci  =4 
and  C2  =  0  (with  ocean  wave  effect).  The  study  shows  that  the  model  which  takes  the  wave 
effect  into  account  gives  better  Co  in  comparison  with  observed  data.  Model  results  are  in 
good  agreement  with  the  few  available  observations  of  TKE  budget,  friction  velocity,  drag 
coefficient,  roughness  length  and  the  accepted  understanding  of  the  influence  of  the  WL. 

The  geostrophic  drag  coefficient  {Cog)  and  wind  speed  at  10-m  height  [Uiq)  determined 
by  the  E-e  method  agree  very  well  with  observations.  The  E-e  method  gives  even  better 
neutral  Co  than  the  ^-method  in  comparison  with  observations  for  all  wind  speeds  greater 
than  5  m  The  u,  obtained  by  the  E-e  method  agrees  very  well  with  observations  for 
winds  less  than  15  m  Roughness  length  zq  also  agrees  we'l  with  observations  and  the 
mixing-length  method  results. 

5  Summary 

WL  strongly  influences  turbulent  and  air-sea  interaction  characteristics.  WL  is  an  additional 
source  of  TKE  for  the  lower/upper  parts  of  .AOBL.  As  expected,  the  atmospheric  and  oceanic 
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'urbuience  properties,  more  than  the  air-sea  interaction  characteristics,  are  sensitive  to  the 
'urbulence  scheme  used.  The  E-i  turbulence  scheme  with  appropriate  constants  can  give 
good  results  in  modeling  coupled  air-sea  boundary-layer  flows.  It  must  be  noted  that  the 
ertect  of  surface  waves  on  energy  transfer  at  the  air-sea  interface  is  extremely  complex  and 
requires  additional  study.  The  next  effort  is  the  modeling  of  the  ocean  breaking  wave  effect 
by  using  E-z  turbulence  scheme  in  a  coupled  model. 
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Figure  Captions 


Fig.  1  The  published  neutral  drag  coefficient.  C d-  as  a  function  of  wind  speed  at  10-m 
height,  f'lo  (after  Ly  et  ah.  1992).  Curve  10  shows  from  the  model  taking  surface 
waves  into  account  with  Ci  =  4  and  C>  =  0. 
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Recent  concentrated  research  in  the  Eastern  Mediterranean  involving  both  observa¬ 
tional  and  modeling  efforts  has  led  to  substantial  new  information  and  understanding  of 
the  general  circulation  and  its  variabilities.  A  new  conceptual  picture  is  emerging.  As  part 
of  phase  I  of  the  international  research  program  for  Physical  Oceanography  of  the  Eastern 
Mediterranean  (POEM)  (Malanotte-Rizzoli  and  Robinson,  1988),  a  comprehensive  pooled 
hydrographic  data  base  was  collected  during  1985-1987  by  scientists  from  several  insti¬ 
tutions  and  nations.  Multi  spatial  and  temporal  scales  are  associated  with  the  general 
circulation  and  its  variabilites  (POEM  group,  1991).  Three  dominant  and  interactive  spa¬ 
tial  scales  have  been  identified:  full  basin  scale,  sub-bzisin  scale  and  mesoscale.  It  has 
been  confirmed  that  a  single  slow  deep  vertical  thermohaline  cell  circulates  throughout  the 
Ionian  and  Levantine  basms  (Roether,  et  a/.,  1991).  (Quasi)-permanent  sub-basin  scales 
gyres  are  interconnected  by  jets  and  currents  to  form  the  basin  scale  horizontal  general 
circulation  of  the  main  thermocline  (Robinson,  et  ai,  1991).  An  energetic  mesoscale  eddy 
field  e.xists  and  interacts  with  the  sub-beisin  scale  features  with  variability  occurring  on 
many  time  scales  (Robinson,  et  ai,  1987).  In  this  paper  we  will  briefly  overview  novel  ex¬ 
perimental  results  on  the  three  dominant  scales  and  introduce  the  data  assimilative  local 
and  basinwide  dynamical  studies  underway  at  Harvard  for  the  Eastern  Mediterranean. 
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Figure  1  shows  the  data  coverage  of  the  POEM  coordinated  dataset  in  summer  of 
1987,  which  is  comprised  of  coordinated  CTD  surveys  with  nominal  1/2  °  latitudinal 
and  longitudinal  spacing,  each  carried  out  within  a  few  weeks.  In  conjunction  with  this 
quasisynoptic  general  circulation  survey,  two  basinwide  surveys  related  to  both  the  large, 
long  time  deep  thermohaline  circulation  and  to  the  synoptic/mesoscale  variability  were 
carried  out.  A  deep  transient  tracer  experiment  was  carried  out  by  the  F/S  Meteor. 
Consistent  with  the  small  Rossby  deformation  radius  of  about  12  km,  620  XBTs  were 
taken  every  10  km  along  the  ship  track  to  provide  a  uniformly  sampled  mesoscale  dataset 
throughout  the  basin. 

The  basinwide  thermohaline  circulation  involves  deep  water  formation  in  the  Adriatic, 
which  then  spills  through  the  Straits  of  Otranto.  An  E-W  section(section  A-B  in  figure  1) 
of  F-12  concentration  is  displayed  in  figure  2.  Deep  water  convective  renewal  rate  has  been 
estimated  to  be  0.29 ±.09  m  ^ /sec,  approximately  15  to  25  %  of  the  exchange  through  the 
straits  of  Sicily.  The  single  coherent  convective  cell  which  connects  the  two  basins  has  a 
turn  over  time  of  126  years  below  1200  meters  (Roether,  ct  ai,  1991). 

The  general  circulation  is  shown  on  figures  3a  and  3b,  respectively,  by  the  dynamic 
height  anomaly  of  the  upper  thermocline  layer  (30m  relative  to  450m)  and  the  associated 
relative  geostrophic  flow  pattern  (Robinson,  et  a/.,  1991).  In  general  the  flow  is  quite 
complex  and  variable.  Figure  3a  shows  a  jet  of  Atlantic  water  entering  through  the  Straits 
of  Sicily,  meandering  through  the  interior  of  the  Ionian  sea,  which  we  believe  feeds  the 
Mid-Mediterranean  Jet,  flowing  through  the  central  Levantine.  The  Mid-Mediterranean 
Jet  bifurcates,  one  branch  flows  towards  Cyprus  and  then  northward  to  feed  the  Asia  Minor 
Current,  and  a  second  branch  separates,  flows  eastward  and  then  southward.  Important 
sub-basin  features  are:  the  Rhodes  cyclonic  gyre;  the  Mersa  Matruh  anticyclonic  gyre;  and 
the  southeastern  Levantine  system  of  anticyclonic  eddies,  among  which  is  the  recurrent 
Shikmona  eddy  south  of  Cyprus.  Flow  in  the  upper  thermocline  is  of  the  order  10-20 
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cm/ sec  as  shown  in  figure  3b.  A  synthesis  of  the  features  and  their  variabilities  is  shown 
in  figure  4.  In  summary,  the  general  circulation  consists  of  sub-basin  scale  gyres  and 
eddies  interconnected  by  jets  and  filaments.  Important  variabilities  include:  1)  shape, 
position  and  strength  of  permanent  sub-basin  gyres  and  their  unstable  lobes,  multi  centers, 
mesoscale  meanders  and  swirls,  2)  meander  pattern,  bifurcation  structure  and  strength  of 
permanent  jets,  3)  occurrence  of  transient  eddies  and  aperiodic  eddies,  jets  and  filaments. 

It  was  found  by  Robinson  et.al.  (1987)  that  a  mesoscale  e.xisted  in  the  Levantine 
basin  consisting  of  energetic  eddies,  jets  and  filaments.  This  result  was  confirmed  by 
the  1987  mesoscale  experiment  and  extended  throughout  both  the  Levantine  and  Ionian 
basins.  Figure  5  is  a  sample  mesoscale  vertical  temperature  section  (section  CEDE  in 
figure  1).  This  figure  reveals  both  mesoscale  structure  in  sub-basin  scale  features  and  an 
energetic  eddy  field  dynamically  interacting  with  the  general  circulation.  Mesoscale  eddies 
occur  with  diameters  on  the  order  of  40-80  km  and  notable  features  include  energetic  sub¬ 
basin/  mesoscale  interaction  in  the  Levantine  basin  and  a  remarkable  thermostad  in  the 
Ionian. 

In  a  dedicated  modeling  effort  at  Harvard  University,  the  POEM  coordinated  datasets 
are  assimilated  into  dynamical  models.  The  Harvard  modeling  hierarchy  consists  of  an  eddy 
resolving  Quasi-Geostrophic  dynamic  model(Robinson  and  Walstad,  1987)  with  attached 
surface  boundary  layer  (Walstad  and  Robinson,  1991),  and  coastal  (Milliff,  1990;  Milliff 
and  Robinson,  1991)  and  island  (Ozsoy,  Lozano  and  Robinson,  1991)  constraints  and  a 
primitive  equation  dynamics  option  (Spall  and  Robinson,  1989). 

Assimilation  of  the  hydrographic  data  into  the  dynamical  models  is  carried  out  in 
order  to  perform  dynamical  adjustment,  dynamical  interpolation  of  fields  into  data  gaps, 
and  to  achieve  physical  process  and  realistic  simulations.  Detailed  parameter  and  sensi¬ 
tivity  studies  are  carried  out  and  thorough  energy  and  vorticity  balances  of  dynamically 
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filtered  data  reveal  dynamical  processes  which  are  occurring.  Basinwide  and  (embedded) 
subregional  assimilations  and  simulations  are  carried  out.  Figure  6  shows  the  complex  to¬ 
pography  of  the  Eastern  Mediterranean  and  the  domains  of  two  e.xamples,  to  be  discussed 
here,  namely  the  Mersa  Matruh  gyre  region  and  the  full  Levantine  basin. 

Studies  have  been  carried  out  in  the  region  of  the  Mersa  Matruh  anticyclonic  gyre 
in  order  to  determine  the  barotropic  mode,  to  fill  in  a  datagap,  and  to  study  the  local 
dynamics.  Figures  7a  and  7b  show  the  streamfunction  field  at  40  meters  and  the  corre¬ 
sponding  error  field  using  the  data  available  in  the  region  from  the  October/November 
1985  (ON85)  .  Figure  7c  is  the  simulated  field  of  stream  function  at  40  meters  after  11 
days.  Extended  parameter  and  sensitivity  studies  indicate  deep  barotropic  flow  of  2-3 
cm/sec,  and  successful  dynamical  interpolation  and  robustness  of  data. 

A  series  of  simulations  have  been  performed  to  investigate  robustness  of  features  in 
the  data,  to  determine  the  distribution  of  the  bznotropic  mode  amplitude  and  stability 
and  variability  of  the  fields  in  the  Levantine  basin.  Figures  8a  and  8b  are  simulated  fields 
of  streamfunction  at  30  meters  after  14  days  of  dynamical  interpolation  for  the  Octo¬ 
ber/November  1985  and  the  March/April  1986  datsets,  respectively.  Quasi-permanent 
features,  namely  the  Rhodes  Gyre,  Mersa  Matruh  gyre  and  the  southeastern  Levantine 
system  of  anticyclonic  eddies  are  notable  in  both  figures  8  and  9.  However,  variabilities  in 
the  quasi-permanent  gyres  and  jets  do  occur.  Our  studies  to  date  have  yielded  dynamically 
adjusted  absolute  flow  and  indicated  robustness  and  stability  of  the  sub-basin  scale  fea¬ 
tures.  This  provides  the  basis  for  internal  dynamical  process  studies  and  the  investigation 
for  their  competition  with  directly  forced  responses. 

In  summary,  circulation  and  variability  of  the  Eastern  Mediterranean  is  now  known  to 
be  comprised  of  three  dominant  interactive  scales.  A  single  deep  convective  cell  connects 
the  Levantine  and  Ionian  basins,  with  an  estimated  turnover  time  of  126  years  below  1200 
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meters.  The  thermocline  horizontal  general  circulation  is  comprised  of  sub-basin  scale 
features,  gyres  and  interconnecting  jets  with  variabilities,  and  fluctuations  consisting  of 
recurrent  and  transient  structures.  Mesoscale  eddies,  meanders  and  oscillations  exist  and 
interact  with  the  sub-beisin  scale,  more  energetically  in  the  Levantine  than  the  Ionian  basin. 
Assimilation  of  POEM  hydrographic  data  into  Harvard  models  and  a  series  of  numerical 
simulations  indicate  robust  and  stable  sub-basin  scale  features  and  provide  the  basis  for 
dynamical  process  studies  including  local  dynamical  processes  and  sub-basin/mesoscale 
interactions.  Internal  dynamical  processes,  cross-scale  interactions,  multiple  and  variable 
forcing  functions,  topographic,  coastal  and  air-sea  interactions  are  among  the  factors  to 
be  investigated  in  terms  of  these  assimilated  fields.  Such  assimilated  data  will  also  provide 
the  best  estimates  of  transport  fields  for  chemical,  biological  and  management  studies  and 
applications. 
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Figure  Captions 

Fig.  1  Station  locations,  August-September  1987(AS  87).  The  line  of  the  deep  stations  ob¬ 
tained  from  the  F/S  Meteor  starting  at  A  and  running  to  B  at  almost  constant 
latitude.  CBDE  is  the  shortest  path  along  the  connected  Meteor  track  through  the 
annotated  points. 

Fig.  2  \  Freon-12  vertical  cross  section,  in  the  Ionian  Sea. 

Fig.  3  a)  Dynamic  height  anomaly  (dm),  upper  thermocline  (30/450)  with  coastal  constraint, 
.\S  37.  b)  Horizontal  flow  vectors  at  30  m,  AS  87. 

Fig.  4  Schematic  general  circulation.  Dashed  features  are  recurrent  or  transient,  upper  thcr- 
mocline. 

Fig.  5  .4  vertical  temperature  cross  section  from  the  mesoscale  XBT  experiment. 

Fig.  6  The  bottom  topography  of  the  Eastern  Mediterranean  Sea  (contour  interval  is  500 
meters),  two  subregions  of  interest,  namely  the  Mersa  Matruh  gyre  region  and  Die 
Levantine  basin  are  shown. 

Fig.  7  a)  Initial  held  of  Stream  Function  at  40  meters  using  data  from  October-November 
1985  (ON  85).  b)  Associated  error  field,  c)  Simulated  stream  function  at  40  meters 
after  14  days. 

Fig.  8  Simulated  field  of  stream  function  at  30  meters  after  14  days,  a)  for  ON  85,  b)  for 
March-April  1986. 
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Figure  3  b 


u 


u 

<N 


Ui 

oo 

CN 


u 

« 

O 

<N 


UJ 

OO 


A-47 


Figure  4 
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1.  INTRODUCTION 

The  Mediterranean  Sea,  although  relatively  small,  is  of  sufficient  size  for  its  circulation 
to  be  governed  by  large-scale  oceam  dynamics.  It  is  among  the  most  interesting  of  our 
planet's  semi-encloseed  seas  because  of  the  great  range  of  processes  and  interactions  which 
occur  within  it  (Robinson  and  Golnaraghi,  1991,  hereafter  RG91).  Most  physical  processes 
which  characterize  the  global  general  ocean  circulation,  many  of  which  are  not  well  known 
or  understood,  occur  directly  or  emadogously  in  the  Mediterranean.  All  major  forcing 
mechamisms,  surface  wind  and  buoyancy  fluxes  and  lateral  mass  exchange,  are  present. 
Air-sea  interaction  is  vigorous  aind  both  deep  water  masses  and  intermediate  water  masses 
are  formed.  Increased  density  due  to  saltiness  caused  by  intense  surface  evaporation  is 
important.  The  analogy  to  the  global  ocean  process  is  that,  where  deep  water  masses  form 
in  polar  regions,  dense  brine  is  present  as  the  result  of  ice  formation.  Freezing  sea  water 
produces  fresh  ice  and  salt. 

The  Mediterranean,  connected  to  the  .Atlantic  by  the  shallow  Straits  of  Gibraltar,  is 
composed  of  two  similar  size  basins  separated  by  the  shallow  and  narrow  Straits  of  Sicily. 
The  Eastern  Mediterranean  is  thus  itself  an  isolated  basin  with  multiple  forcings  and  water 
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mass  formations  driving  the  processes  and  interactions  of  global  interest.  Logisticaily  the 
basin  is  convenient  for  operations.  Resource-wise,  it  is  both  observationally  and  computa¬ 
tionally  feasible  to  conduct  basin- wide  research  in  the  Etistern  Mediterranean  with  grids 
of  fine  enough  resolution  to  define  the  physical  processes  of  interest.  These  factors  provide 
a  unique  opportunity  for  oceanographic  research  of  general  interest.  Such  research  has 
been  underway  now  for  several  years,  coordinated  by  the  international  POEM  (Physical 
Oceanography  of  the  Eastern  Mediterranean)  program  (Malanotte-Rizzoli  and  Robinson, 
1988:  The  POEM  Group.  1991).  There  has  resulted  recent  rapid  progress  in  knowledge 
zmd  understanding.  .Although  certain  aspects  of  the  circulation  have  been  known  generally 
or  speculatively  for  some  time,  locations,  rates  and  dynamics  of  water  mass  formations 
and  transformations  are  beginning  to  emerge.  Furthermore,  a  new  picture  of  the  scales, 
structures  and  dynamics  of  the  horizontal  general  circulation  and  its  variabilities  is  be¬ 
coming  apparent.  Thus  the  basis  now  exists  for  research  in  understanding  climate  related 
dynamical  processes  of  global  relevance,  in  detail  and  in  a  basin-wide  context. 

Modeling,  of  course,  plays  a  significant  role  in  modern  ocean  scientific  research.  Mul¬ 
tiple  and  interacting  scsdes  and  processes  require  a  multi-faceted  modeling  approach.  Fur¬ 
thermore,  the  assimilation  of  data  into  models,  the  construction  of  melded  estimates  from 
models  and  observations,  and  data-based  simulations  now  provide  efficient  mechanisms 
for  exploiting  observations.  Data  assimilative  models  are  not  only  useful  research  tools. 
They  also  provide  the  best  mechanism  for  estimating  the  four-dimensional  fields  required 
for  transport  estimates  for  heat,  chemicals,  particles,  etc.  necessary  for  climate  related 
research  and  other  interdisciplinary  marine  science  and  technology.  realistic  set  of  high 
resolution  multiscale  basin-wide  management  models  can  be  based  upon  the  research  mod¬ 
els  developed  for  general  circulation  and  climate  change  research. 

The  following  sections  will  discuss  processes,  describe  the  multiscale  circulation,  and 
illustrate  data  assimilation  for  the  Eastern  Mediterranean  Sea. 
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2.  PROCESSES 


Processes  of  global  relevance  for  ocean  climate  dynamics  include  aspects  of  the  ther¬ 
mohaline  circulations,  water  mass  formations  zind  transformations,  dispersion  and  mixing. 
Processes  are  sketched  schematicaily  on  Figure  1.  Fresh  water  of  Atlantic  origin  (AW) 
enters  through  the  Straits  of  Sicily.  As  it  circulates,  driven  by  the  wind  (r)  and  other 
forces,  evaporation  (E)  much  exceeds  precipitation  at  the  surface  and  the  water  becomes 
more  sadine  and  dense.  Net  heating  (Q)  occurs  in  summer  and  cooling  in  winter.  In  the 
Levantine  basin,  predominantly  in  March,  Levantine  Intermediate  Water  (LIW)  is  formed 
during  storm  events,  the  newly  formed  water  mass  circulates  and  disperses  at  depths  down 
to  a  few  hundred  meters.  Some  LIW  recirculates  within  the  basin  and  some  exits  beneath 
the  .\W  through  the  Sicily  Straits.  The  ratio  of  recirculation  to  eflux  is  not  well  known 
but  the  net  inflow  and  outflow  transports  through  the  Sicily  Straits  are  estimated  between 
1-1.5  Sverdrup  (10®  m^/sec).  LIW  is  the  predominant  source  of  the  Mediterranean  water 
masses  in  the  Atlantic. 


Figure  1.  Processes  of  water  mass  formations,  dispersions  and  transformations. 
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Deep  water  formation  occurs  primaiily  in  the  Adriatic  Sea  which  produces  bottom 
water,  the  Eastern  Mediterranean  Deep  Water  (EMDW).  However.  Aegean  Deep  Water 
'  .AGDW)  also  forms  which  sinks  to  depths  somewhere  in  the  middle  of  the  water  column. 
The  site  of  EMDW  formation  in  winter  is  the  southern  Adriatic  which  is  several  hundred 
meters  deep.  Very  cold  and  dry  air  storms  apparently  cause  deep  convection  throughout 
the  water  column.  The  detailed  formation  process  is  not  well  known  but  it  involves  surface 
water  admixed  with  LIW  which  has  entered  from  the  Mediterranean.  The  newly  formed 
EMDW  then  exits  the  .Adriatic  through  the  Straits  of  Otranto,  plunges  to  the  bottom,  and 
moves  off  along  the  deep  western  boundary  of  the  Ionian  basin.  .A  single  vertical  thermoha- 
line  cell  exists  for  the  whole  Eastern  Mediterranean,  although  the  partitioning  of  the  rising 
motion  between  the  Ionian  and  Levantine  basins  is  not  known.  Other  essential  unknowns 
are  the  locations  and  rates  of  mixing  between  waters  of  LIW  and  EMDW  origins.  This  is 
symbolized  by  the  crossed  corkscrew  arrow's  on  Figure  1.  There  is  also  no  observationally 
based  quantitative  estimate  at  present  of  the  amount  of  (mixed)  EMDW  which  exits  at 
the  base  o'  the  LIW  flow  through  the  Straits  of  Sicily. 

The  thermohaline  circulations,  which  aire  included  on  Figure  1  are  redrawn  for  clarity 
on  Figure  2.  For  simplicity  the  schematic  mixing  of  EMDW  and  LIW  has  been  omitted. 
There  is  a”  internal  cell,  a  single  vertical  cell  encompassing  both  the  Ionian  and  Levantine 
basins,  and  an  external  cell.  The  external  thermohaline  cell  involves  an  exchange  of  water 
through  the  Straits  of  Sicily.  The  remote  e.xcursions  of  EMDW  and  the  return  paths  in 
the  Western  Mediterranean  and  .Atlantic  Ocean  are  not  known. 

There  is  special  interest  attached  to  the  internal  cell  in  that  it  is  a  closed  thermohaline 
cell  in  which  the  deep  water  formation  involves  a  preconditioning  process  of  density  increase 
due  'o  saltiness.  Thus  it  is  to  some  extent  an  analog  of  the  global  scale  "conveyor  belt" 
thermohaline  cell  whose  deep  water  forms  in  the  North  .Atlantic  (Gordon.  1QS6:  Broeker. 
19S7  i.  .A  detailed  understanding  of  the  Eastern  Mediterranean  ceil  should  be  achievable 
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Figure  2.  The  thermohaJine  cells. 

in  the  near  future  which  should  be  expected  to  contribute  insights  into  globally  relevant 
processes. 

3.  THE  MULTISCALED  CIRCULATION 

It  is  now  known  that  the  general  circulation  and  its  variabilities  consists  of  three  pre¬ 
dominant  spatial  scades  which  must  be  understood  together  with  their  nonlinear  interac¬ 
tions.  This  concept  is  depicted  schematically  in  Figure  3.  The  slow  (vertical)  thermohaline 
cell  is  basin  scale.  The  basin  scale  general!  circulation  of  the  main  thermocline  of  the  upper 
ocean  is  composed  of  subbasin  scale  gyres  linked  by  subbasin  scale  jets.  There  are  a  vari¬ 
ety  of  subbasin  scale  variabilities  which  will  be  described  below.  The  energetic  mesoscale 
which  is  present  is  illustrated  in  Figure  3  on  one  of  the  subbasin  scale  gyres  by  a  field  of 
interned  eddies  and  a  meandering  of  the  outer  swirl  flow,  and  as  the  meandering  of  a  jet 
segment. 
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Figure  3.  Scales  of  circulation  variabilities  and  interactions. 

The  POEM  program  conducted  five  quasi-synoptic  nominally  baisin-wide  hydrographic 
surveys  from  1985  through  1987  (Robinson  et  ai,  1991).  The  sampling  patterns  were 
grids  with  nominal  spacings  of  one  half  degrees  of  latitude  and  longitude,  which  provides 
fine  general  circulation  (subbasin)  scale  and  coarse  mesoscale  resolutions  (Robinson  et  a/.. 
1987).  The  best  coverage,  achieved  with  six  ships  during  August-September,  1987.  is  shown 
in  Figure  4.  A  beisin-wide  multiscale  experiment  was  conducted  simultaneously  from  the 
F/S  Meteor,  whose  track  is  superimposed  over  the  survey  grid  on  Figure  4.  Transient  tracer 
measurements  were  made  to  investigate  the  deep  thermohaline  circulation  and  XBTs  were 
taken  uniformly  every  ten  kilometers  along  the  ship  track  to  explore  the  characteristics  of 
mesoscale  phenomena  throughout  the  beisin.  CTD  stations  were  also  carried  out. 

3.1  The  Basin  Scale  Thermohaline  Cell 

The  transient  tracer  experiment  consisted  of  62  stations  where  chlorofluoromethane 
I  CFM)  and  tritium  were  measured  !  Schlitzer  et  al..  1991 ).  The  CFM  Freon  12  distribution 
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Figure  4.  Multiscaie  survey  data  sampling  map  (from  RG91  Figure  1). 

cdong  the  east-west  section  AB  shown  on  Figure  4  is  presented  on  Figure  5.  The  source  of 
Freon  12  is  the  atmosphere,  and  there  is  of  course  a  large  concentration  in  the  upper  ocean. 
But  there  is  cdso  a  maximum  in  the  Freon  distribution  in  the  deepest  water  in  the  western 
basin  (off  the  Malta  escarpment)  which  is  clear  evidence  for  the  presence  of  recently  formed 
deep  water  of  Adriatic  origin.  The  minimum  concentration  is  at  mid-depths  in  the  eastern 
Ionian  basin. 

In  order  to  obtain  quantitative  information  about  the  deep  water  renewal  rate  from  this 
data  set.  Roether  and  Schlitzer  ( 1991)  have  constructed  a  22  box  model  with  compartments 
for  the  Ionian  and  Levantine  deep  water  below  1000  m  depths.  They  find  that  the  water 
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Figure  5.  Freon  12  section  across  the  Ionian  basin  (pmd  kg  ’ )  (from  Rcether  and  Schlitzer. 
1991  Figure  2). 

formed  in  the  Adriatic  is  a  mixture  of  surface  water  (AW)  and  intermediate  water  (LIW) 
from  the  Mediterranean  which  therefore  has  also  previously  been  subject  to  surface  fluxes. 
Their  results  indicate  unambiguously  that  the  thermohaline  circulation  consists  of  a  single 
vertical  cell  throughout  the  entire  Eastern  Mediterranean,  flowing  coherently  throughout 
both  the  loniein  and  Levantine  basins.  The  turnover  time  for  the  deep  water  is  126  years. 

The  rate  of  deep  water  supply  is  0.29  i  09  Sverdrups  (10^  m^/s).  This  is  about 
one  third  of  the  inflow  or  outflow  rate  of  all  waters  through  the  Sicily  Straits.  There  is 
some  qualitative  evidence  observationally,  from  transient  tracer  measurements  (Roether 
and  Schlitzer.  1991)  that  some  EMDW  does  exit  through  the  Straits  of  Sicily.  The  only 
quantitative  information  comes  from  a  numerical  general  circulation  model  for  the  East¬ 
ern  Mediterranean  which  has  active  thermodynamics  (Malanotte-Rizzoli  and  Bergamasco, 
1991 !.  The  results  must  be  regarded  as  tentative;  the  stratification  is  imposed  and  the 
deep  convective  process  is  not  modeled.  The  interesting  result  however,  is  that  about  one 
third  of  the  water  efluxing  through  the  Straits  of  Sicily  is  deep  water  which  has  upwelled  in 
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the  southern  Ionian  basin.  At  this  time  it  would  seem  reasonable  to  assume  that  internal 
and  external  thermohaline  cells  are  of  comparable  magnitude. 

3.2  The  Subbasin  Scale  Components  of  the  General  Circulation 

The  currents  and  circulation  of  the  upper  ocean  thermocline  are  shown  on  the  ob¬ 
jectively  analyzed  dynaunic  height  (30  m  relative  to  450  m)  and  current  vector  maps  of 
Figures  6a  and  6b,  respectively.  The  complex  pattern  of  cyclonic  and  anticyclonic  gyres 
linked  by  jets  and  current  segments  is  real:  no  contours  are  plotted  where  the  e.xpected 
error  of  the  objective  analysis  is  less  than  fifty  percent.  A  synthetical  study  has  been  car¬ 
ried  out  for  the  five  independent  quasisynoptic  realizations  contained  m  the  pooled  data 
set  of  POEM  (Robinson  et  ai,  1991).  The  results  for  the  circulation  ,  ^  variabilities  are 

summarized  in  the  schematic  of  Figure  7  and  the  companion  Table  1. 

There  axe  a  number  of  permanent  albeit  variable,  features.  Other  features  etre  recur¬ 
rent,  i.e.  disappear  but  reappear  similarly.  Transient  gyres  and  jets  also  occur  randomly 
throughout  the  basin.  Veiriabilities  characterizing  the  permanent  and  recurrent  structures 
include:  the  center,  size  cind  strength  of  gyres  and  their  multilobe  mesoscale  instabilities; 
the  ajcis  location,  meander  pattern  and  bifurcation  structure  of  jets.  The  time  scales  as¬ 
sociated  with  the  subbasin  scaile  variabilities  and  mesoscale  interactions  are  not  yet  well 
defined  but  appear  to  range  from  yeeirs  (e.g.,  recurrence  time  of  the  Shikmona  gyre)  to  days 
(mesoscale  interactions).  General  circulation  modeling  studies  indicate  that  subbcisin  scale 
structures  are  the  form  of  the  Eastern  Mediterranean’s  response  to  the  large-scale  forcing 
functions  (Malaxiotte-Rizzoli  and  Bergamaisco,  1991)  due  to  factors  including  topographic 
effects  and  the  seasonal  power  in  the  wind-stress  forcing  (Pinardi  and  Navarro.  1991). 

3.3  Levantine  Intermediate  Water 

The  formation  of  Levantine  Intermediate  water  was  for  many  years  believed  to  occur 
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F’.sure  6.  ai  Dynamic  height  anomaly  in  the  upper  thermociine.  30  m  relative  to  450 
m  I  dynamic  meters)  August-September  1987.  b)  Associated  horizontal  geostrophic  flow 
vectors  i  scale  vector  cm  sec-l  i  'from  Robinson  et  al.  1991  Figure  7), 
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Figure  7.  General  circulation  schematic  (from  Robinson  et  al.  1991  Figure  2). 

TABLE  1  -  Upper  Thermocline  Features 
(from  Robinson  et  al.  1991  Table  1) 


Features 

Type 

ON  85 

MA  S6 

MA  87 

.AS  87 

lonizin  Atlantic  Stream 

P 

— 

- 

Y 

Y 

Mid  Med.  Current 

P 

Y 

Y 

- 

Y 

.Asia  Minor  Current 

P 

Y 

N 

Y 

N 

Cilician  Current 

R 

Y 

Y 

- 

Y 

Southeast  Levamtine  Jets 

T 

Y 

Y 

- 

Y 

Rhodes  C 

P 

Y 

Y 

Y 

West  Cyprus  C 

P 

Y 

Y 

- 

Y 

Mersa  Matruh  AC 

P 

Y 

Y 

- 

Y 

Cretan  Sea  C 

P 

Y 

9 

- 

Y 

Shikmona  .AC 

R 

Y 

Y 

- 

Y 

Latakia  C 

R 

Y 

N 

N 

Y 

.Antalya  AC 

R 

7 

Y 

- 

N 

Pelops  AC 

7 

- 

Y 

9 

Y 

Ionian  eddies  .AC 

T 

- 

- 

- 

Y 

SE  Levantine  eddies  .AC 

T 

Y 

Y 

- 

Y 

Cretan  Sea  eddies 

T 

Y 

Y 

- 

Y 

Definitions:  P=Permanent.  R=Recurrent,  T=Transient.  "’  =  Not  enough  information  to 
classify,  C=Cyclone.  AC  =  Anticyclone. 
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predommantly  in  a  limited  region  between  Rhodes.  Cypress  and  the  Turkish  coast  and 
to  involve  a  preconditioning  process  afforded  by  the  domed  pycnocline  characteristic  of 
the  Rhodes  gyre.  Recent  evidence,  however,  indicates  that  formation  occurs  over  a  much 
larger  domain  which  could  include  the  centers  of  anticyclonic  gyres  (Oguz,  1991)  via  a 
process  related  to  thermocline  deepening  in  warm  core  Gulf-Stream  rings  (Dewar,  19S6). 
The  horizontal  circulation  structures  and  pattern  observed  and  analyzed  by  Ozsoy  et  al. 
(1989)  for  the  upper  thermocline  of  the  Levantine  basin  in  October  of  1985  is  shown  in 
Figure  Sa.  composite  vertical  salinity  section  around  the  basin  (.\BCD  on  Figure  3a)  is 
shown  on  Figure  8b  with  LIW  (salinity  greater  than  38.95  ppm)  shaded.  LIW  is  pre'.'alent 
and  ubiquitous  throughout  the  basin  including  the  southern  sector. 


3.4  Mesoscale  Variability 

The  internal  Rossby  radius  of  deformation  is  the  horizontal  scale  which  characterizes 
mesoscale  instabilities,  meandering  and  eddying.  It  is  essentially  the  vertical  scale  (e.g. 
thermocline  depth)  of  the  phenomenon  amplified  by  the  ratio  of  the  Brunt-Vaisla  buoyancy 
frequency  to  the  Coriolis  parameter.  Since  this  length  scale  is  only  ~  10-12  km  in  the 
Eastern  Mediterranean,  four  times  smaller  than  typical  values  for  much  of  the  worlds  ocean, 
mesoscale  phenomena  were  not  defined  in  observations  until  recently  when  dedicated  high 
resolution  sampling  was  undertaken  (Robinson  et  ai,  1987). 

The  mesoscale  experiment  conducted  in  1987  by  .A.R.  Robinson,  W.G.  Leslie  and  N. 
Pinardi  took  620  XBT  samples  at  nominal  10  km  spacing  along  the  entire  FS  Meteor 
ship  track.  A  composite  section  across  the  basin  from  the  mouth  of  the  .Adriatic  Sea  to 
the  southeastern  Levantine  corner  (CEDE  on  Figure  4)  is  shown  on  Figure  9.  .4  number 
of  interesting  features  are  present  which  characterize  the  highly  resolved  structures  of 
the  subbasin  scale  gyres  and  jets  (e.g.  the  Rhodes  gyre  in  the  vicinity  of  D).  Mesoscale 
eddies  e.xist  in  the  open  sea  (e.g.  northwest  of  E)  and  interactively  with  the  subbasin  scale 
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100  km 

Figure  9.  Composite  XBT  derived  temperature  section  along  track  CBDE  of  Figure  4 
(from  RG91  Figure  5). 

features.  The  Levantine  basin  appears  more  energetic  in  the  mesoscale.  with  the  Ionian 
(C-B)  characterized  by  a  remarkable  thermostad  at  14.0-14.5°C  from  75  to  450m. 

4.  DATA  ASSIMILATION  INTO  DYNAMICAL  MODELS 
4.1  Methodology 

The  Harvard  dynamical  model  hierarchy  consists  of  a  set  of  compatible  and  inter¬ 
comparable  numerical  dynamical  ocean  models  constructed  under  physical  assumptions 
suitable  for  various  regions  and  phenomena  in  the  oceans.  Domains  and  geometries  are 
general,  rangi^’g  fmm.  an  a.rbitrary  piece  of  open  ocean  (water-water  boundaries)  to  an 
entire  basin  scale.  The  model  hierarchy  approach  was  adopted  in  order;  to  achieve  the 
closest  possible  relationships  of  modeling  to  experimental  data  bases;  to  allow  the  focus  of 
computing  resources  and  scientific  ertort  on  regions  of  special  interest;  and  to  facilitate  the 
physical  interpretation  of  results  via  utilization  of  the  simpie,st  physical  model  required  for 


a  simulation  of  specified  accuracv. 
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The  present  model  hierarchy  consists  of  a  quasigeostrophic  model  for  a  block  of  open 
ocean  i  Robinson  and  W'alstad.  19S7)  or  for  a  region  (partially)  bounded  by  an  arbitrary 
coastline  (Milliff.  1990)  which  contains  islands  (Ozsoy.  et  al.  1991).  A  surface  boundary 
layer  may  be  attached  i  Walstad  and  Robinson.  1991).  More  general  primitive  equation 
dynamics  may  be  utilized  with  arbitrary  open  boundary  conditions  (Spall  and  Robinson. 
1989). 

The  POEM  data  sets  have  been  assimilated  into  the  dynamical  models  and  studies 
are  ongoing.  Initialization  of  the  dynamical  model  with  a  real  ocean  data  set  serves  to 
dynamically  adjust  the  data  and  generally  improves  the  field  estimiates.  Data  gaps  can 
be  interpolated  dynamically.  Dynamics  together  with  hydrography  yields  absolute  flow 
estimates;  dynamics  provides  the  consistent  barotropic  mode.  Four  dimensional  fleld  es¬ 
timates  best  suitable  for  the  transport  flows  necessary  as  input  to  interdisciplinary  and 
applied  studies  are  provided.  Furthermore,  the  results  of  simulations  based  on  real  data 
sets  can  be  analyzed  in  detail  qualitatively  and  quantitatively  in  order  to  reveal  real  phys¬ 
ical  processes  operative  in  the  ocean.  Four  exeimples  of  the  Eastern  Mediterraneain  will  be 
briefly  described  for  the  regions  shown  on  Figure  10:  the  full  beisin,  the  eastern  Levantine 
basin,  and  the  Mersa  .Matruh  and  Rhodes  gyres. 

4.2  Results 

Figure  11  shows  the  stream-function  field  (proportional  to  upper  thermocline  dynamic 
high)  for  the  winter  1986  POEM  quasisynoptic  survey  after  fourteen  days  of  dynamical 
adjustment.  The  assimilation  is  remarkably  robust.  Overall,  the  features  resemble  the 
objective  analysis  initialization  (not  shown).  There  has  been  dynamical  smoothing  and 
natural  variabilities  (meandering  of  the  mid-Mediterranean  jet)  are  underway.  The  re¬ 
markable  robustness  of  the  subbasin  scale  features  to  dynamical  assimilation  in  this  data 
set  is  characteristic  of  all  the  data  sets. 
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Figure  10.  Modeling  domains  and  bottom  typography. 

Dynamical  interpolation  is  illustrated  for  the  Mersa  Matruh  gyre  in  the  region  indi¬ 
cated  in  Figure  11  (which  is  to  appear  in  the  doctoral  dissertation  of  Ms.  M.  Golnaraghi  of 
Harvard  University).  The  objective  analysis  initialization  of  Figure  12a  can  be  seen  from 
the  expected  error  maps  of  Figure  12b  to  suffer  from  an  absence  of  data  in  approximately 


Figure  11.  Streamfunction  held  at  40  m  depth  14  days  after  initialization.  The  stream 
tuncf.on  IS  proportional  to  dynamic  height  ( nondimensionai  units  i  i  from  RG91  Figure  S  i. 
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Figure  12.  a)  Messah  Matruh  stream  function  at  40  m  objective  analysis  map.  b  i  .Associ¬ 
ated  expected  error  map.  c)  .As  (a)  after  14  days  of  dynamical  run  (from  RG91  Figure  i  i. 


one  third  of  the  domain  of  interest  north  of  the  Egyptian  coastline.  .After  fourteen  days 
(Figure  12c)  dynamics  provides  a  consistent  complete  gyre  circulation. 

The  1987  mesoscale  experiment  accomplished  a  mapping  pattern  in  the  region  of  the 
Rhodes  gyre  as  evidenced  by  the  zig-zag  ship  track  in  that  region  on  Figure  4.  This 
served  as  the  basis  of  a  simulation  and  dynamical  process  study  of  mesosctde-subbasin 
scale  interactions  for  the  gyre  system  (Milliff  and  Robinson.  1991).  The  main  thermocline 
stream-function  for  three,  seven  and  ten  days  into  the  simulation  are  shown  in  Figure  13. 
There  is  a  four-lobed  oscillation  of  the  gyre  center  with  a  period  of  about  a  week.  This 
appecLTS  to  be  driven  by  intermittent  bursts  of  baroclinic  instability  associated  with  the 
tightening  of  the  jet  between  the  cyclonic  Rhodes  gyre  and  the  anticyclonic  gyre  to  its 
northeast,  which  can  be  seen  most  intensively  on  Figure  12c. 

.A  full  basin  assimilation  of  the  summer  19S7  data  shown  on  Figure  6  is  presented  in 
Figure  14.  Dynamical  interpolation  and  wind  forcing  have  filled  in  the  large  data  gap  in 
the  southern  Ionian  and  also  the  gaps  south  of  Crete  and  in  the  mid-Levantine  basin.  The 
flow  depicted  is  forty-five  days  into  a  dynamical  run  initialized  with  the  objective  analysis 
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Figure  13.  Mesoscale  simulation  of  Rhodes  gyre  system  after  a)  3  days,  hj  7  days,  c!  10 
days  after  MillifF  and  Robinson  1991  Figure  6). 


of  Figure  6a  and  also  driven  by  the  August-September  meam  winds.  The  Atlantic  Ionian 
Stream  which  projects  to  the  northeast  into  the  Ionian  threads  its  way  across  the  basin 
continuously  to  the  Mid-Mediterranean  Jet.  This  is,  to  our  knowledge,  the  first  time  that 
observations  and  dynamics  have  been  combined  to  produce  a  quasisynoptic  realization  of 
a  full  basin  circulation  with  subbasin  scale  resolution. 


5.  SUMMARY  AND  CONCLUSIONS 

The  Eastern  Mediterranean  Sea  is  an  open  but  isolated  basin  of  interest  as  a  laboratory 
basin  for  global  ocean  general  circulation  processes  including  oceanic  climate  dynamical 
processes.  The  circulation  and  variabilities  of  the  Eastern  Mediterranean  Sea  are  known 
to  be  constituted  from  multiple  and  interactive  space  and  time  scales.  .A.  single  full  basin 
scale  thermohaline  cell  exists  with  a  turnover  time  of  one  and  one  quarter  centuries.  The 
main  thermoclinic  general  circulation  is  built  up  from  subbasin  scale  features,  cyclonic 
and  anticyclonic  gyres  linked  by  interconnecting  jets.  A’ariabilities  include  fluctuations  in 
the  characteristics  of  the  permanent  features  as  well  as  the  occurrence  of  recurrent  and 
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Figure  14.  Full  basin  simulation  45  days  after  initialization.  Lpper  thermociine  function. 

transient  features.  Mesoscaie  eddies,  meanders  and  oscillations  occur  and  interact  with 
the  subbasin  scale  circulation,  most  energetically  in  the  Levantine  basin  (in  the  summer 
1987  data  set). 

Important  climate  processes  include  the  formation,  dispersion,  transformation,  mixing 
and  circulation  of  both  intermediate  and  deep  water  masses.  The  basin  wide  deep  ther¬ 
mohaline  circulation  has  both  an  internal  cell,  confined  to  the  Eastern  Mediterranean,  and 
am  external  cell  which  communicates  with  the  world  ocean  through  the  Western  Mediter¬ 
ranean.  The  internal  cell,  whose  deep  water  formation  involves  preconditioning  via  salt 
concentration,  may  present  interesting  analogies  to  the  global  conveyor  belt. 

The  assimilation  of  the  new  data  sets  into  models  provides  absolute  flow  from  hydro- 
graphic  data  plus  dynamics,  and  improves  the  data  via  processes  of  dynamical  adjustment 
and  dynamical  interpolation.  The  subbasin  scale  general  circulation  features  are  stable 
and  robust  under  assimilation.  Data  driven  extended  simulations  and  process  studies  for 
both  internal  dynamics,  general  circulation  forcings  and  multiscale  nonlinear  interactions 
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are  in  progress.  The  realistic  synoptic  time  series  of  fields  estimates,  provided  regionally 
or  basin  wide  by  such  models,  also  provide  best  estimate  transport  fields  for  biological, 
chemical  and  management  studies  and  applications. 
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abstract 


A  priminve  eqiianon  ocean  model  which  makes  use  of  a  ciin-iiinear  orthogonal  grid  and 
r  a  sigma-coordinate  system  was  used  to  study  the  general  circulation  oj  the  Meaiterranean 
Sea.  The  model  was  forced  with  monthly  values  of  wind  stress,  heat  and  salinity  flux 
aiues.  In  one  experiment,  the  fresh-water  river  discharge  of  the  Rhone.  Ebro  and  Po  rivers 
was  also  used. 

The  model  simulates  the  circulation  reasonably  well  and  the  results  indicate  that,  in  the 
western  Mediterranean,  the  wind  stress  is  very  important  in  establishing  the  summer 
northward  shift  of  the  .Atlantic  inflow.  The  lateral  boundary  forcing,  the  surface  heat  and 
laiinir.-  flux  are  necessary  to  the  maintenance  of  the  cyclonic  circulation  in  the  northern 
Balearic  Basin  and  enhance  the  summer  reversal  of  the  circulation  in  the  Tyrrhenian  Sea. 
particularly  interesting  result  is  the  existence  of  a  seasonal  variation  in  the  Levantine 
Intermediate  Water  path  . 

In  the  Ionian  Sea  ,  the  dynamic  consequence  of  the  wind-stress  and  the  surface 
buoyancy  flux  tend  to  offset  each  other  and  the  path  of  the  Ionian  Atlantic  stream  is 
seasonally  invariant.  In  the  Levantine  Sea  the  dominant  forcing  mechanism  is  probably  linked 
to  the  buoyancy  fluxes,  but  we  do  not  consider  our  results  conclusive  because  it  is  probable 
that  the  climatological  buoyancy  fluxes  are  underestimated. 

The  Levantine  Intermediate  water  formation  in  the  north-central  Levantine  Sea  (Rhodes 
gyre)  and  the  deep  water  formation  in  the  southern  Adriatic  Sea  are  reproduced  quite  well. 
However,  deep  water  formation  in  the  western  Mediterranean  is  not  produced  by  the  model. 
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/;  I.\TRODLCTIO.\ 


This  paper  deals  with  the  simulanon  of  the  general  circulanon  of  the  Medlterranear.  5ea 
ana  augments  previous  basin-scale  and  regional  modeling  efforts.  On  a  basin  scale, 
modeling  with  pnmitive  equation  models  were  carried  out  by  Menzin  and  Moskalenko 
:9S2].  who  studied  the  barocropic  wind  driven  circulation.  Stanev  et  al.  i  1989)  utilized  a 
coarse  resolution  model  to  investigate  the  response  of  the  circulation  to  annual  and  seasonal 
rorcing  while  the  recent  work  of  Pinardi  and  Navarra  1 1991)  focused  on  the  baroclinic  wind 
caven  circulation  in  a  basin  with  simplified  bathymetry'.  .An  inverse  model  was  used  by 
Tziperman  and  Malanone-Rizzoli  (1991)  to  calculate  the  seasonal  climatological  circulation 
rrom  hydrological  data. 

Regional  circulation  studies  were  carried  out  by  Heburn  (1987)  with  a  one-layer 
reduced  gravity  model  of  the  western  Mediterranean  and  by  Malanotte-Rizzoli  and 
Bergamasco  (1989;  1991)  who  performed  an  extensive  set  of  e.xperiments  with  a  primitive 
equation  model  to  study  the  barocropic  and  baroclinic  seasonal  circulanon  of  the  eastern 
Mediterranean.  Preller  (1986)  carried  out  a  regional  study  for  the  Alboran  Sea  and 
Malanotte-Rizzoli  and  Bergamasco  (1982)  for  the  Adriatic  Sea. 

In  this  study,  we  focused  on  the  seasonal  variability  of  the  circulation  and  on  water 
mass  formation  processes  to  identify  the  most  imponant  forcing  functions  and  the  areas  and 
mechanisms  involved  in  the  water  mass  transformation  processes. 

The  following  section  is  devoted  to  a  description  of  the  Mediterranean  Sea. 
hydrological  and  oceanographical  characteristics  in  the  light  of  the  most  recent  observational 
findings.  Sections  3,  4  and  5  describe  the  model,  the  bathymetry,  grid  and  the  surface 
forcings  respectively.  Section  6  and  7  describe  the  characteristics  of  the  numerical 
experiments  and  the  calculated  results.  Finally  in  section  8  we  offer  discussion  and 
conclusions. 
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2:  BACKGROUSD 


The  Ntedicerranean  Sea  (Fig.  1)  is  an  enclosed  Basin  connected  to  the  Atlantic  Ocean 
cy  the  narrow  (width  -  13  km)  and  shallow  (sill  depth  -  300  mi  Strait  ot  Gibraltar  and  the 
Black  Sea  is  connected  by  the  Dardanelles  -  Marmara  Sea  -  Bosphorus  system.  It  lies 
cetween  30°  and  46°  N  and  between  -4°  and  36°  E  and  is  essentially  formed  by  two  sub- 
basins.  the  western  (WIMED)  and  the  eastern  (EJvtED)  Mediterranean,  connected  by  the  Strait 
of  Sicily  whose  sill  depth  is  ~  300  m.  and  width  ~  35  km. 

The  W\IED  has  a  triangular  shape  with  the  Strait  of  Gibraltar,  the  island  of  Sicily  and 
•.he  Ligurian  Sea  at  its  ape.xes.  Excepting  the  Gulf  of  Lions,  the  continental  shelves  are 
f'arrow.  The  Tyrrhenian  Sea.  between  the  Italian  peninsula  and  the  islands  of  Sardinia  and 
Corsica,  is  the  easternmost  and  deepest  (depth  -  3500  m)  pan  of  the  Basin  and  joins  the  rest 
of  the  WMED  to  the  south  in  a  wide  opening  between  Sardinia  and  Sicily.  The  Balearic  Sea 
(depth  ~  2500  m)  occupies  the  central  pan  of  the  Basin  and  extends  nonheast  to  the  Ligurian 
Sea  and  west,  toward  Gibraltar,  to  the  Alboran  Sea. 

The  EMED  has  a  more  complicated  structure  than  the  WMED.  Four  sub-basins  can  be 
defined:  the  Ionian,  the  Levantine,  the  Adriatic  and  the  Aegean  Seas.  The  Ionian  Sea  lies 
between  southern  Italy  and  Greece  to  the  north  and  Lybia  and  Tunisia  to  the  South,  has  a 
depth  of  about  3500  -  4000  m,  with  a  maximum  of  5000  m  south  of  Greece.  The  Gulf  of 
Sine  is  the  southernmost  area  of  the  Mediterranean  Sea.  The  Levantine  Sea  has  a  roughly 
rectangular  shape  with  depths  of  about  2500  -  3000  m  in  the  center  of  the  Basin  and  a 
maximum  depth  of  about  4500  ra  in  a  depression  located  southeast  of  the  island  of  Rhodes.  It 
merges  with  the  Ionian  Sea  through  the  Cretan  Passage  (depth  >  1500  m),  between  Crete  and 
the  Lybian  coast,  fhe  Adriatic  Sea  is  connected  with  the  Ionian  by  the  Strait  of  Otranto 
(width  ~  75  km,  sill  depth  ~  800m);  its  topography  is  characterized  by  an  extended 
continental  shelf  to  the  north  and  by  a  relatively  deep  southern  section.  The  nonhem  end  of 
the  .Adriaric  (Gulf  of  Venice)  is  the  nonhemmost  area  of  the  Mediterranean  Sea.  The  Aegean 
Sea  joins  the  Levantine  Sea  through  several  passages  and  straits,  located  between  the  Greek 
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-ind  Turkish  coast  and  the  islands  of  Crete  and  Rhodes,  with  depths  and  widths  ranging  from 
.'00  to  “00  m  and  from  15  to  60  km  respectively.  It  has  very  irregular  coastlines  and 
topography  with  many  islands:  its  maximum  depth  is  about  1500  m. 

The  hydrology'  and  the  circulation  of  the  Mediterranean  Sea  waters  are  strongly  affected 
by  the  sunace  energy  and  mass  exchange  with  the  atmosphere;  the  fresh-water  budget  is 
negative;  i.e..  the  loss  by  evaporation  exceeds  the  input  by  precipitation  and  river  run-off. 
The  esrimated  fresh- water  deficit  amounts  to  about  2100  km^/year,  corresponding  to  a  sea 
level  drop  of  about  1  m/year  and  to  a  residence  time  of  about  75  -  100  years  {Hopkins.  197S: 
Beihonx.  1979:  Tchernia.  1980  V  The  fresh-water  deficit  (about  S.S'Td  of  the  total  inflowing 
volume)  is  mainly  replenished  by  the  .Atlantic  surface  waters  flowing  into  the  .Mediterranean 
having  a  salinity  of  36.25°/oo-  The  excess  of  salt  is  balanced  by  a  deep  outgoing  flow  of  high 
salinity.  38.40o/oo- 

The  e.xistence  of  the  two  superposed  flows  in  the  Strait  of  Gibraltar  requires  a 
transformation  of  the  inflowing  .Atlantic  waters  and  vertical  movement  which  provide  the  link 
between  the  two  flows  and  which  determines  the  formation  of  three  distinct  water  masses 
whose  vertical  distribution  is  schematized  in  Fig.  2.  The  surface  waters  flowing  into  the 
-Vlediterranean  are  subject  to  evaporation  and  mixing  with  the  underlying  waters,  causing  a 
progressive  increase  of  the  salinity;  the  surface  value  rises  from  36.25°/oo  in  the  Gibraltar 
area  to  37.25o/oo  in  the  Strait  of  Sicily  and  to  values  higher  than  38.50o/oo  in  the  Levantine 
Sea.  Such  water  forms  the  Modified  Atlantic  Water  (MAW).  Its  west  to  east  path  across  the 
Mediterranean  can  be  tracked  by  a  subsurface  salinity  minimum  {Lacombe  and  Tchernia. 
1 960)  progressively  deepening  from  20  to  50  m  depth,  representing  the  signature  of  their 
.Atlantic  origin. 

The  MAW  overlies  the  Levantine  Intermediate  Water  (LIW).  Their  depth  range  is  300  - 
"^00  m  in  the  WMED  and  200  -  400  m  in  the  EMED.  They  are  the  result  of  convection 
processes  in  the  EMED  occurring  in  the  Rhodes  -  Cyprus  area,  and  probably  in  other  sectors 
of  the  Levantine  Sea  in  winter  {Marcos.  1972).  .At  their  source  the  LIW  have  a  salinity  of 


about  39.IO0/00  ^nd  spreads  to  the  whole  Mediterranean  (S  =  38.70°/oo  at  the  Stratt  of 
Sictiv).  reaching  Gibraltar  where  they  form  the  Mediterranean  outflow  into  the  Atlantic 
Ocean. 

The  deep  waters  of  the  WMED  and  of  the  E\LED  are  the  result  of  intensive,  deep 
convection  processes,  but  their  sources  are  distinct  and  separate.  .VIoreover  their  interr'ace 
depth  with  the  LIVV  is  deeper  than  the  depth  of  the  sill  of  Sicily,  preventing  any  direct 
e.xchange  between  them.  The  source  of  the  western  Mediterranean  deep  water  is  located  in  the 
Gulf  of  Lions  area  and  the  water  mass  formation,  the  so-called  MEDOC  deep  water 
formarion.  takes  place  in  winter  (MEDOC  Group.  1969;  Sankey.  1973:  Leaman  and  Schot:. 
1991'].  There,  intensive  convective  movements  occur  under  the  influence  of  cold  and  dr\' 
winds,  causing  the  sinking  and  mixing  of  the  relatively  cold  and  salty  surface  waters  to 
depths  of  about  1200  -  1500  m.  The  resulting  water  mass  has  a  potential  temperature.  0.  of 
about  12.70°  C  and  a  salinity  of  38.40°/oo-  Such  water  mass  can  enter  the  Mediterranean 
outflow  at  Gibraltar  without  mixing  with  the  LIW  {Kinder  and  Parrilla,  1987).  The  source  of 
the  eastern  Mediterranean  deep  waters  is  the  Adriatic  Sea  {Poliak,  1951),  where  a  water  mass 
with  9  =  13.60°  C  and  S  =  38.70°/oo  is  formed  in  winter,  again  under  the  influence  of  wind, 
and  spreads  to  the  whole  EMED  deep  Basin.  The  traditional  interpretation  of  the  deep  water 
formation  in  the  Adriatic  Sea  assigns  a  key  role  to  the  cooling  and  evaporative  processes 
acring  on  the  shelf  waters  of  the  nonhem  Adriatic,  causing  their  sinking  along  the  continental 
slope.  However,  Ovchinnikov  et  al.  (1987)  observed  a  MEDOC-like  deep  water  formation 
in  the  southern  Adriatic  producing  dense  waters  flowing  into  the  Ionian  Sea  over  the  Otranto 
sill,  while  the  dense  water  formed  in  the  northern  Adriatic  were  not  involved  in  the  exchange. 

Recent  international  research  projects,  such  as  the  "Western  Mediterranean  Circulation 
experiment".  WMCE  {La  Violate,  1990),  "Mediprod  5"  {Millot,  1991)  and  the  "Physical 
Oceanography  of  the  Eastern  Mediterranean",  POEM  {Robinson  et  al.,  1991),  provide  new 
and  important  information  about  the  general  circulation  of  the  Mediterranean  Sea.  revealing 
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general  patterns  more  complex  than  those  depicted  by  earlier  studies  (Lacombe  and  Tchernia. 
I960:  Ovchinnikov  and  Fedoseyev.  1965:  Ovchinnikov.  1966'). 

The  sunace  and  intermediate  circulation  of  the  WlVtED  is  illustrated  in  Fig  3.  according 
:o  Milloi  1 1987a:  b:  1991).  At  the  surface  the  Atlantic  water  entering  the  Basin  forms,  in  the 
Alboran  Sea,  two  anticyclonic  gyres  (Heburn  and  La  Violetta,  1990).  The  eastern  side  of  the 
eastern  .Alboran  gyre  constitutes  the  .Almeria-Oran  front  (Tintore  et  al.,  1988:  .Arnone  e:  al.. 
1990).  a  strong  density  front  (representing  the  contact  between  the  inflowing  .Atlantic  water 
and  the  resident  water  of  the  Mediterranean)  which  deflects  the  MAW  flow  toward  the 
African  coast  and  originates  the  Algerian  current  flowing  along  the  .Algerian  coast  between 
0'^  and  2°  E.  The  current  becomes  progressively  unstable  and  generates  mesoscale  eddies 
both  cyclonic  and  andcyclonic.  Observations  (Millot,  1985)  showed  that  only  the  anticyclonic 
eddies  survive,  increase  in  size,  propagate  eastward  and  can  leave  the  coast  and  drift  into  the 
Balearic  Basin.  Thus  the  MAW  do  not  only  flow  eastward  toward  the  Strait  of  Sicily,  but 
can  be  advected  nonhward  by  the  eddies,  finally  coming  in  contact  with  more  modified 
MAW  in  the  northern  Balearic  Basin,  giving  rise  to  the  north  Balearic  front.  Such  a  pattern 
also  implies  that  the  MAW  flow  into  the  EMED  is  controlled  by  "eddy  shaped  sporadic 
events"  {Millot,  1987a). 

The  LIW  circulation  in  the  WMED  has  been  the  subject  of  many  studies  and 
speculations.  It  is  accepted  that  at  the  exit  of  the  Strait  of  Sicily  The  LIW  turns  around  the 
western  rip  of  Sicily  into  the  Tyrrhenian  Sea  where  it  cyclonically  flows  along  the  Italian 
coast.  Before  entering  the  Ligurian  Sea,  the  LIW  current  originates  a  branch  flowing  south 
and  subsequently  turning  north  around  Sardinia.  According  to  Wust  (1961)  and  Ovchinnikov 
( 1966),  the  LIW  exiting  the  Tyrrhenian  Sea  originates  a  further  branch  proceeding  directly 
westward  to  the  Strait  of  Gibraltar,  explaining  the  presence  of  a  LIW  layer  off  the  North- 
African  coasts.  This  view  has  been  questioned  by  Katz  (1972)  and  Millot  (1985)  on  the 
basis  of  field  studies  and  simple  theory.  The  LIW  presence  in  the  Algerian  Basin  is  now 
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assumed  :o  result  from  sporadic  erosion  of  the  LIW  current  by  the  mesoscale  eddies 
ongmating  from  instability  of  the  Algerian  current. 

The  c’j-culation  of  the  EMED,  reiadve  to  the  upper  thermocline.  is  shown  in  Fig.  4.  The 
general  circulation  appears  to  be  dominated  by  sub-basin  gyres  interconnected  by  jets.  The 
NLAW  enters  the  EMED  and  flows  through  the  Ionian  Sea  ( lonian-.Atlantic  Stream)  into  the 
Levantine  Sea  i\Iid-Mediterranean  Jet),  where  the  following  gyres  are  present  throughout  the 
year  {Oszoy  et.  ai..  1989:  1991):  a  cyclonic  gyre  between  Rhodes  and  Cyprus  (Rhodes 
gyre)  and  one  south  of  Crete  (Cretan  gyre),  an  anticyclonic  gyre  off  the  coast  of  Egypt 
(.Mersa  -Matruh  gyre)  and  a  wide  anticyclonic  area  in  the  north-eastern  Levantine  Basin 
Shdtmona  eddy).  The  Mid- Mediterranean  jet  flo%vs  between  the  Rhodes  and  the  Mersa- 
.Macruh  gyres  and  then  bifurcates  twice,  east  and  west  of  Cyprus.  The  nonhward  branch 
originating  west  of  Cyprus  feeds  the  Asia-Minor  current,  while  the  southward  branch  enters 
the  south-eastern  Levantine  anticyclonic  region.  Although  not  shown  in  Fig.  4.  another 
permanent  feature  in  the  EMED  circulation  is  the  southern  Adriatic  cyclonic  gyre  {Malanotie- 
Rizzoli  andHecht,  1988). 

3:  THE  OCEAN  MODEL  AND  THE  GRID. 

The  ocean  model  is  the  Princeton  ocean  model,  described  in  detail  by  Bliimberg  and 
Mellor  (1987).  The  description  of  the  model  computer  code  can  be  found  in  Mellor  (1991a). 
It  has  been  applied  to  the  South  Adantic  Bight  {Blumberg  and  Mellor,  1983),  Delaware  Bay 
{Galperin  and  Mellor  1990a,  b),  the  Hudson-Raritan  Estuary  {Oey  et  al.,  1985 ..  b,  c  ),  the 
Gulf  of  Mexico  (Blumberg  and  Mellor,  1985),  and  the  Gulf  Stream  {Ezer  and  Mellor,  1991) 
and  some  thirty  or  more  additional  oceanic  regions.  It  is  a  primitive  equation  model  with  a 
free  surface,  a  split  mode  time  step  and  solves  the  following  equations  for  the  ocean  velocity 
Ui  =  (LLV.W),  potential  temperature  9  and  salinity  S; 
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The  vertical  mixing  coefficients  K.v/  and  K//  in  (2),  (3)  and  (4)  are  calculated  using  the 
\fellor  and  Yamada  (1982)  turbulence  closure  scheme,  while  the  horizontal  diffusion  terms 
•F;;'.  Ft'),  F^  and  Fj  in  (2),  (3)  and  f4)  are  calculated  using  a  Smagorinsky  horizontal 
diffusion  formulation  {Mellor  and  Blwnberg,  1985).  The  term  R  in  (3)  is  the  ponion  of  the 
shon-wave  radiation  which  penetrates  the  sea  surface. 

The  hydrostatic  approximation  yelds: 

where  P  is  the  free  surface  elevation,  Po  the  reference  density  and  p  =  /?  (0,  S,  P)  is 
the  density  calculated  by  an  adaptation  of  the  UNESCO  equation  of  state  by  Mellor  (1991b) 
which  is  computationally  simpler  and  which  has  potential  temperature  rather  than  in  sim 
temperature  as  its  argument. 

The  model  uses  a  bottom- following  sigma  coordinate  system  (  cr  =  (z  -  77)/(H  +  p  ) 
where  Hfx,  y)  is  the  bottom  topography).  The  transformed  equations  can  be  found  in  the 
papers  cited  above.  In  the  present  study  the  model  has  16  vertical  layers,  corresponding  to 
the  following  sigma  coordinates:  0.000,  -0.003,  -0.006,  -0.012,  -0.025,  -0.050.  -0.100, 
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-0.:00.  -0,300.  -O.-OO.  -0.500.  -O.6OO.  -O.TOO.  -0.800.  -0.900.  -1.000. 

The  model  gnd  is  the  cundlinear  onhogonal  grid  shown  in  Fig.  5  and  contains  113  x 
39  points.  The  grid  resolution  ranges  from  a  minimum  of  8  to  15  km  in  the  .Alboran  Sea 
region  to  a  maximum  of  50  to  60  km  in  the  nonhemmost  pan  of  the  Basin. 

The  model  has  one  open  boundarx'  in  the  .Atlantic  Ocean  about  200  km  west  of  the  Strait 
of  Gibraltar.  The  normal  velocities  are  governed  by  a  Sommerfeld  radiation  condition. 
Temperature,  salinity  and  tangenrial  velocities  on  the  boundary  are  upwinded:  when  the  flow 
is  into  the  model  dom.ain  they  are  prescribed  from  the  observed  climatology  and  the  boundarx' 
values  are  advected  into  the  model  domain. 

4:  THE  MODEL  BATHYMETRY  AND  THE  INITIAL  CONDITIONS 

The  model  bathymetry  was  obtained  from  the  U.S.  Navy  bathymetric  data  base. 
DBDB5,  by  bilinear  interpolation  of  the  depth  data  into  the  model  grid.  The  resulting 
topography  is  illustrated  in  Fig.  6  and  shows  the  major  features  of  the  basin's  geometry, 
which  retains  the  island  of  Sardinia,  Corsica,  Sicily,  Mallorca,  Crete,  Cyprus  and  Rhodes, 
The  connections  with  The  Black  Sea  is  closed. 

The  model  has  been  initialized  with  the  temperature  and  salinity  annual  climatological 
averages  prepared  by  the  National  Institute  for  Oceanography  and  Fisheries  at  Alexandria 
(Egypt),  using  the  NODC  (National  (Dceanographic  Data  Center)  and  the  WDC  (World  Data 
Center)  data-bases.  The  annual  averages  are  mapped  on  a  1/4°  x  1/4°  grid  at  the  following  27 
depth  layers:  0,  10,  20,  30,  50,  75,  100,  125,  150,  200,  250,  300,  400,  500,  600,  700, 
800,  900,  1000,  1100,  1200,  1300,  1400,  1500,  2000,  2500,  3000  m  depth.  The  grid 
points  with  missing  data  were  filled  by  interpolation  of  the  surrounding  values.  The  surface 
temperature  and  salinity  annual  averages  are  shown  in  Fig.  7. 
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5;  THE  SURFACE  BOL'yDARY  CONDITIONS 

The  surface  boundarv'  condirions  are  surface  heat  flux,  surface  salinity  flux  and  wind- 
stress.  The  monthly  surface  heat  flux  and  wind-stress  averages  were  obtained  by  May  i  I9S2: 
1986).  These  climatologies  are  mapped  on  a  1®  x  1°  grid.  The  surface  heat  flux  is:  short 
wave  radiation.  Qj,  long-wave  back  radiation,  Q/,  sensible  heat  flux,  Q/,,  and  latent  heat  flux. 
Qt;:  the  total  heat  flux,  Q; ,  is  then  given  by: 

Q;  =  Qr  ~  Q/  +  Q/i  +  Q.e  (  6) 

In  the  model  Q/  +  Q/i  +  Qg  plus  a  portion  of  the  shon-wave  radiation  is  absorbed  at  the 
surface  where  the  boundary  condition  is: 

{Kh~)  =  (l-Tr)Q,+  Q/+  Qa+  a 

dz  (7) 

the  remainder  is  directly  added  to  the  heat  equation  (3)  as  9R/3z,  where: 

R  =QrTr  exp(Xz)  (8) 

The  transmission  coefficient  Tr  =  0.31  and  the  attenuation  coefficient  X  =  0.042  m*^ 
are  chosen  to  faJerlov's  (1976)  tabulated  attenuation  functions  for  "clear "  water  type. 

The  total  heat  flux  fields  for  the  months  of  February  and  August,  calculated  according 
to  (6),  are  shown  in  Fig.  8.  The  February  field  shows  a  large  cooling  area  over  the  Gulf  of 
Lions  (WMED),  while  the  EMED  is  characterized  by  a  zonal  pattern  with  the  largest  negative 
values  in  the  nonhem  pan  of  the  Basin  (around  Rhodes).  In  August,  the  heat  flux  is  positive 
everywhere  but  again  the  minimum  values  are  found  in  the  Nonhem  Balearic  Basin  and  in 
the  Rhodes  region. 
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The  '.vind-stress  fieids  for  the  same  months  are  liiustrated  in  Fig.  9  and  e.xempiify  the 
'.vind  regime  over  the  Mediterranean  throughout  the  year.  One  sees  evidence  of  the  gaps 
cetween  tne  mountain  ranges  bordenng  the  northern  Mediterranean  coasts.  The  most 
prominent  winter  feature  characterizing  the  W\IED  is  the  occurrence  of  a  strong  northerly 
wind,  the  .Mistral,  blowing  from  the  Rhone  river  valley  through  the  opening  between  the 
.Alps  and  the  Pyrenees.  In  the  EMED  the  wind-stress  is  mostly  zonal.  In  the  nonhern  .Adnanc 
there  occurs  an  easterly  wind,  the  Bora,  blowing  over  the  sea  from  the  opening  between  the 
.Alps  and  the  Dynaric  .Alps;  this  wind  can  be  very  strong  but  it  is  not  as  frequent  as  the  Mistral 
and.  climatologically.  it  results  in  only  a  slightly  larger  stress.  In  .August,  the  .Mistral  still 
affects  the  WA/IED  although  with  reduced  stress.  In  the  EMED.  the  dominant  w'inds  a.^e  the 
Etesians,  blowing  from  the  opening  between  the  Balkans  and  the  .Anatolian  mountains  over 
the  Levanrine  Sea  with  a  cyclonic  curvature.  The  wind-stress  direction  over  the  Adriatic  and 
the  Ionian  Sea  is  mainly  meridional. 

The  monthly  evaporation  rate  (E)  was  calculated  from  May's  latent  heat  flux  data 
according  to: 


E  =  ^  (9) 

Lv 

where  Ly  is  the  latent  heat  of  condensation.  The  monthly  precipitation  rate  (P)  was  obtained 
from  Jaeger  (1976),  whose  data  are  mapped  on  a  5°  x  2.5°  grid.  Finally  the  fresh- water 
surface  flux  (Ws)  was  calculated  according  to: 

W^  =  E-P  (10) 

and  the  surface  boundary  condition  for  salinity  is: 

(K^^)  =WsSo  (11) 

a:  :  =0 
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The  calculated  total  annual  fresh-water  loss  amounts  to  92.0  cm/year  (a  value  which 
agrees  reasonably  well  with  the  results  of  many  budget  calculations)  and  the  monthly 
averages  range  between  a  minimum  of  3.6  cm/month  (May)  and  a  ma.ximum  of  10.5 
cm/month  (  August).  The  fields  for  the  months  of  February  and  August  are  shown  in 
Fig.  10  and  are  very  similar  to  the  total  heat  flux  fields  of  Fig  S..  the  reason  for  which  is  the 
dominant  influence  of  the  latent  heat  flux  component  over  both  fields,  as  pointed  out  by 
Malanoiie-Rizzoli  and  Bergamasco  11991)  in  their  analysis  of  surface  flux  fields. 

.All  the  surface  boundary  condirions  were  applied  to  the  model  by  interpolaring  linearly 
between  adiacent  monthly  values. 

In  one  numerical  experiment  (see  below),  the  annually  averaged  fresh-water  nver 
discharge  for  the  Rhone.  Ebro  (WMED)  and  Po  (Adriatic  Sea)  rivers  were  added  to  the 
model.  The  averages  values  (Rhone  =  1712  m^/s.  Ebro  =  550  m^/s.  Po  =  1550  m^/s)  are 
those  reponed  in  UNEP  (1984). 


6:  THE  NUMERICAL  EXPERIMENTS 

We  give  here  a  brief  description  of  the  numerical  experiments  performed  using  different 
kinds  of  surface  forcings.  Their  characteristics  are  summarized  in  Table  1. 


Table  1:  Characierisacs  of  the  numerical  experiments. 


EXPERIMENT 

WIND-STRESS 

THERMAL  SURF. 
CONDITION 

SALINITY  SURF. 
CONDITION 

RIVER 

RUNOFF 

.VlED-1 

Monthly 

Surf,  temperature 

Surf,  salinity 

no 

equal  to  initial 

equal  to  initial 

conditions 

conditions 

.MED-2 

Monthly 

Monthly 

Constant 

no 

heat  flux 

salinity  flux 

MED-3 

Monthly 

Monthly 

Surf,  salinity 

no 

heat  flux 

equal  to  initial 

conditions 

MED-A 

Monthly 

Monthly 

Monthly 

Annual 

heat  flux 

salinity  flux 
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Expenment  NLED-i:  In  this  expenmeni  the  model  was  forced  witn  the  montni\  '.v  mu- 
stress  averages,  while  the  surface  temperature  and  salinity  were  prescnbed  to  be  equal  to  the 
annual  averages  unmai  conditions)  in  order  to  simulate  perpetual  annual  average  conditions. 

E.xpenment  MED-2;  The  forcing  were  the  monthly  averaged  wind  stress  ana  the 
monthly  sun'ace  heat  flux.  The  surface  salinity  flux  was  constant  in  time  and  space  and  was 
obtained  by  prescribing  a  surface  fresh-water  flux  of  1  m/'year.  corresponding  to  the 
estimated  annual  W5  average  for  the  whole  Mediterranean  Sea  iBeihoiLX.  1979). 

Experiment  MED-3;  The  wind-stress  and  surface  heat  forcing  were  analogous  to 
E.xpenment  MED-2.  but  the  surface  salinity  value  w-as  prescnbed  to  be  equal  to  the  annual 
averages. 

E.xpenment  .MED-4  (the  central  experiment):  the  monthly  climatological  averages  were 
used  for  all  the  surface  fluxes  and  the  fresh-water  inputs,  described  above,  were  introduced. 

The  integration  dme  was  5  years  for  experiments  MED-1,  MED-2,  .VIED-3  and  6  years 
for  experiment  ,MED-4,  but  we  show  here  the  results  relative  to  the  fifth  year  of  integration  in 
order  to  be  consistent  with  the  other  experiments;  the  results  relative  to  the  sixth  year  results 
do  no  show  significant  differences  with  those  of  the  fifth  year. 

7;  RESULTS 
Experiment  MED-1 

In  experiment  MED-1,  we  examine  the  effects  of  the  seasonal  wind-stress  forcing  by 
eliminanng  any  seasonality  in  the  surface  heat  and  salinity  fluxes,  fhe  surface  elevation  fields 
and  the  velocity  field  at  50  m  depth  for  the  months  of  February  and  August  are  shown  in 
Figs.  11  and  12  respectively.  The  surface  circulation  in  the  WMED  is  dominated  by  the 
MAW  inflow  along  the  Nonh  African  coast  and  by  the  cyclonic  circulation  in  the  northern 
Balearic  Basin  and  in  the  Tyrrhenian  Sea.  The  water  incoming  from  the  Atlantic  Ocean  flows 
m  the  .Alboran  Sea  along  the  coast  of  Spain  and  bend  southward  at  the  exit  of  the  Sea. 
forming  a  front  resembling  the  Almeria-Oran  front,  from  which,  in  winter,  originates  a 


A-88 


jurren:  moving  along  the  Afncan  coast.  In  summer  the  circulation  in  the  southern  Baieanc 
Basin  undergoes  a  major  change;  the  flow  along  the  African  coast  is  reversed  and  two 
anticyclonic  gyres  are  formed.  Throughout  the  year,  the  flow  from  the  .Atlantic  Ocean  feeds 
both  the  cyclonic  circulation  in  the  T>TTheman  Sea  and  the  flow  into  the  EMED.  In  winter, 
the  \LAW  lonian-.Atlantic  Scream  enters  the  Basin  via  the  Strait  of  Sicily,  bends  south  and 
subsequently  bifurcates  into  two  branches,  one  closely  following  the  .African  coast  and 
directly  entering  the  Levantine  Sea.  the  other  meandering  in  the  Ionian  Sea  a..d  becoming 
entrained  in  the  large  cyclonic  gyre  developed  in  the  nonhem  Ionian  Sea.  In  summer,  the 
lonian-.Atlandc  Scream  shifts  r'l'nhward  along  the  eastern  coast  of  Sicily  and  is  entrained  into 
a  broad  anticyclonic  gyre,  occupying  the  whole  central  part  of  the  Basin.  The  jet  is  bloc.ked 
from  tlowing  into  the  Levantine  Sea  by  a  cyclonic  gyre  comparable  vvith  the  Cretan  gyre 
described  by  the  POEM  data  (Pig.4).  The  coastal  following  branch  of  the  winter  jet  reverses 
its  direction  in  summer.  In  the  Levantine  Sea,  the  circulation  along  the  coast  is  cyclonic 
throughout  the  year  and  appears,  in  winter,  directly  connected  to  the  lonian-Atlantic  stream, 
while  in  summer  the  current  weakens  and  originates  from  the  Cretan  gyre,  in  the  interior  of 
the  Basin  a  meandering  jet  comparable  with  the  Mid- Mediterranean  jet  flows  eastward, 
bordering  a  wide  cyclonic  gyre  whose  location  and  shape  correspond  to  the  Rhodes  gyre.  In 
summer  other  anricyclonic  gyres  develop;  one  south-east  of  Crete,  one  south  of  the  Rhodes 
gyre  similar  to  the  Mersa-Matruh  gyre  and  one  south  of  Cyprus,  similar  to  the  Shikmona 
gyre.  Funhermore  two  interconnected  cyclonic  gyres  close  to  the  Egyptian  coast  that  appear 
as  in  the  POEM  maps  of  Ozsoy  et  al.  (1989)  and  in  the  geostrophic  calculations  by  El-Din 
and  Karam  ( 1 976 ) . 

The  circulation  at  intermediate  depth  in  February  and  August  is  shown  in  Fig.  13.  In 
the  \VMED  the  LIW  circulation  look  very  similar  to  the  pattern  proposed  by  Millot  (1985) 
and  described  in  section  2.  The  return  flow  from  the  EJvLED  in  winter  crosses  the  Strait  of 
Sicily  and  turns  around  Sicily  entering  the  cyclonic  circulation  of  the  Tyrrhenian  Sea.  from 
which  the  LIW  e.xits.  then  flows  around  Sardinia,  circulates  into  the  cyclonic  gyre  in  the 
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Baleanc  Basin  and  reaches  the  Alboran  Sea  moving  along  the  French  and  Spanish  coasts.  In 
summer,  the  cyclonic  circulation  in  the  Tyrrhenian  Sea  becomes  very'  weak  and  the  LIW 
b>"passes  this  Basin  and  flows  direcdy  into  the  Balearic  Basin  along  the  coast  or  Sardinia. 

In  the  EMED.  more  precisely  in  the  Ionian  Sea.  there  are  seasonal  differences  in  the 
ongin  of  the  LFW  flow  into  the  WMED.  In  winter,  two  currents,  coming  from  the  Southern 
and  the  Nonhern  Ionian,  join  south  of  Sicily  to  form  the  outflow  and  there  is  no  direct 
contribution  from  the  Levantine  Sea.  In  summer,  the  flow  along  the  eastern  coast  of  Sicily 
weakens  and  the  flow  toward  the  WMED  is  formed  mainly  by  the  current  from  the  southern 
Ionian,  which  is  also  partially  fed  by  a  flow  exiting  the  Levantine  Basin  and  recirculating 
into  the  cyclonic  gyre  developed  in  the  Gulf  of  Sirte.  E.xcept  for  this  current,  the  circulation 
of  the  Levantine  Sea  is  more  or  less  constant  throughout  the  year  and  is  mainly  characterized 
by  a  cyclonic  gyre  with  a  SW-.NE  elongation  in  the  center  of  the  Basin  centered  o\  er  the  area 
of  maximum  depth  and  by  the  same  anticyclonic  gyres  observed  at  the  surface. 

The  deep  circulation  (Fig.  14)  does  not  show  a  pronounced  seasonal  variation.  In  the 
WlvlED  the  flow  in  the  northern  Balearic  Basin  and  in  the  Tyrrhenian  Sea  is  reversed  with 
respect  to  the  surface  and  intermediate  layers,  while  in  the  southern  pan  of  the  Basin  the  flow 
is  still  directed  eastward. 

In  the  EMED  The  Ionian  Sea  is  fed  by  water  coming  from  the  Adriatic  Sea:  part  of  it 
moves  toward  the  Strait  of  Sicily  and  pan  recirculates  into  the  cyclone  in  the  southern  Ionian 
and  enters  the  Levantine  Sea,  where  the  circulation  is  dominated  by  a  large  cyclonic  gyre, 
occupying  the  greater  pan  of  the  Basin,  whose  shape  and  position  are  clearly  influenced  by 
the  bottom  topography  (Figs.  1  and  6). 

Experiments  MED-2  and  MED-3 

As  a  second  step  in  our  study  we  performed  an  experiment  forcing  the  mode!  with  the 
monthly  surface  heat  flux  and  with  a  constant  salinity  flux  obtained  by  prescribing  a  fresh¬ 
water  flux  of  1  m/year.  Despite  the  utilization  of  a  more  realistic  seasonal  forcing  the  results 


A-90 


of  me  expenmenr  MED-2  were  not  entirely  satisfactory.  The  major  problem  was  that,  after  5 
years  of  integration,  the  cyclonic  circulation  in  the  northern  Balearic  Basin  almost  disappears; 
i.e.  u  IS  very  weak  in  winter  and  practically  non  existent  in  summer  (Fig.  15i;  funhermore 
the  anticycionic  area  developing  in  summer  in  the  meridional  part  of  the  Basin  is  ill  defined. 
The  circulation  features  in  the  EMED  are  more  satisfactory  and  comparable  with  the  results  of 
e.xperiment  MED-1. 

In  order  to  understand  if  the  problems  in  the  northwestern  Basin  were  caused  by  the 
surface  heat  flux  or  by  the  surface  salinity  flux,  experiment  MED-3  (propeny  fields  not 
shown)  was  carried  out.  whereby  the  model  was  driven  with  monthly  averaged  wmd-stress 
and  sunace  heat  flux,  but  the  surface  salinity  was  mantained  equal  to  the  annual  average. 
This  experiment  gave  a  more  satisfactory  result,  indicating  that  the  choice  of  a  constant 
salinity  flux  is  probably  incorrect;  this  led  us  to  the  experiment  which  is  described  below. 

Experiment  MED-4 

Monthly  averaged  surface  salinity  fluxes,  according  to  equation  (11)  and  river  runoff  in 
the  nonh- western  Basin  were  introduced  and  improved  the  performance  of  the  model.  The 
Basin  averaged  temperature  and  salinity  time  history  is  shown  in  Fig.  16  and  indicates  that 
the  model  acquires  a  periodical  seasonal  cycle.  The  free  surface  elevation  and  velocity  surface 
fields  for  the  months  of  February  and  August  are  shown  in  Figs.  17  and  18  respectively.  The 
cyclonic  gyre  in  the  nonhem  Balearic  Basin,  although  weak  in  summer,  is  mantained 
throughout  the  year,  demonstrating  the  importance  of  the  salinity  surface  flux  and  lateral 
boundary  runoff  in  this  area.  With  respect  to  the  results  of  experiment  MED- 1  a  major 
difference  arises  in  the  Alboran  Sea  region  in  winter,  where  two  small  anticycionic  eddies  are 
formed  close  to  the  Moroccan  coast,  while  the  main  flow  of  Atlantic  water  is  in  the  nonh. 
Immediately  west  of  the  Alboran  Sea,  the  model  generates  a  large  and  intense  (particularly  in 
summer)  anticycionic  eddy.  In  the  scientific  literature  there  is  no  strong  evidence  of  the 
existence  of  such  a  gyre.  The  only  indication  comes  from  the  work  of  Green  and  Price 
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1991 ;  who  racked  buoys  deployed  in  the  Alboran  Sea  and  in  the  Southern  Balearic  Basin. 
Their  results  show  that  the  buoys  were  occasionally  entrained  into  an  anticyclonic  eddy 
which  they  associate  with  the  Aimena-Oran  front.  It  is  interesting  to  note  that,  as  in  the 
summer  of  experiment  MED-1.  the  surface  flow  in  the  Southern  Balearic  Basin  is 
characterized  by  the  development  of  anticyclonic  eddies,  which  in  this  experiment  appear  to 
onginate  from  the  described  eddy/frontal  feature  and  to  propagate  eastward.  Thus  it  seems 
that  the  model  is  attempting  to  reproduce  some  features  of  the  instability  processes  of  the 
.Algerian  current  leading  to  the  formation  of  anticyclonic  eddies  observed  in  the  southern 
Baleanc  Basin. 

.Another  important  difference  ber-veen  this  experiment  and  experiment  .VIED-i  occurs  in 
the  Tyrrhenian  Sea,  where,  in  expenment  \IED-4,  and  in  agreement  with  the  results  of  Tail 
1 1984)  and  Krivosheya  and  Ovchinnikov  (1973),  in  the  interior  of  the  Basin  the  circulation 
exhibits  a  strong  winter  to  summer  variability,  reversing  from  cyclonic  to  anticyclonic,  while 
a  nonhward  flow  along  the  Italian  coast  persists  throughout  the  year.  A  consequence  of  such 
reversal  is  the  reduced  summer  inflow  of  the  MAW  in  the  Tyrrhenian  Sea. 

In  the  EMED,  the  seasonal  variability  induced  by  wind  stress,  seen  in  experiment 
MED- 1 ,  now  seems  to  be  offset  by  the  surface  heat  and  buoyancy  fluxes.  In  this  experiment 
the  summer  formation  of  the  large  anticyclonic  gyre  in  the  central  Ionian  Sea  does  not  occur 
and  the  path  of  the  stream  simply  shifts  northward  forming  a  wide  anticyclonic  meander 
extending  from  Sicily  almost  to  the  eastern  coast  of  the  Gulf  of  Sine.  Subsequently  the 
stream  turns  nonh  forming  a  narrow  cyclonic  meander  and  enters  the  Levantine  Sea.  The 
nonhem  ponion  of  the  Ionian  Sea  is  occupied  by  a  cyclonic  gyre,  more  intense  in  winter  and 
weaker  in  summer,  during  which  it  exchanges  water  with  the  permanent  cyclonic  gyre 
located  in  the  southern  Adriatic  Sea  (Malanotte-Rizzoli  and  Hecht,  1988).  The  persistence  of 
the  cyclonic  gyre  in  this  area  of  the  Ionian  Sea  constitutes  one  of  the  major  differences 
between  the  results  of  our  model  and  the  seasonal  maps  of  the  Mediterranean  climatological 
ctrculation  obtained  by  Tziperman  and  Malanotte-Rizzoli  (1991)  utilizing  an  inverse  model; 
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'.heir  results  show  a  seasonal  reversal  of  the  circulation,  being  cyclonic  in  winter  and 
anncycionic  m  summer.  In  Tne  Levantine  Sea  the  most  prominent  permanent  feature  is.  as  in 
expenment  .\LED-1.  the  Rhodes  gyre.  In  the  Cretan  Passage  the  Mid-Mediterranean  jet  forms 
a  cyclonic  meander,  but  a  clear  development  of  an  anticyclonic  gyre,  comparable  to  the 
Mersa-.Matruh  gyre,  does  not  occur.  Another  absent  feature  that  is  the  Shikmona  gyre  which 
developed  (weakly)  in  experiment  MED-1;  in  this  experiment  only  a  small  anticyclonic 
meander  appears  in  front  of  the  Israeli  coast.  The  eastern  Mediterranean  model  developed  by 
Maianotte-Rizzoli  and  Ber^amasco  (1991)  exhibits  the  same  behaviour.  They  atmbuted  this 
model  deficiency  to  the  ^^ay  f  1 982 )  sunace  heat  flux  data,  whose  magnitude  they  thought  to 
be  underestimated.  This  may  be  the  case,  but  the  Shikmona  eddy  may  also  be  controlled  by 
bottom  topography  and  particularly  by  the  Eratosthenes  sea-mount,  rising  from  the  sea 
bottom  (-2000  m)  to  -600  m  depth,  (Brenner  et.  al.,  1991).  Our  model  barely  resolves  such 
topographical  feature. 

Despite  these  possible  discrepancies,  other  surface  features  are  represented  reasonably 
well,  in  particular  the  nonhward  branching  of  the  Mid-Mediterranean  Jet  west  and  east  of 
Cyprus,  both  branches  feeding  the  Asia- Minor  Current  flowing  westward  along  the  Turkish 
coast  and  becoming  the  nonhem  side  of  the  Rhodes  gyre. 

.At  intermediate  depth.  Fig.  19,  the  circularion  pattern  in  the  WMED  substantially 
reproduce  the  patterns  obtained  in  experiment  MED-1,  with  the  imponani,  exception  of  the 
Tyrrhenian  Sea.  In  particular,  it  is  noted  again  that  no  direct  LIW  flow  from  the  Strait  of 
Sicily  to  the  Alboran  Sea  develops,  the  path  being  blocked  by  an  anticyclonic  gyre.  The  LIW, 
once  again,  reaches  the  Strait  of  Gibraltar  by  flowing  along  the  northwestern  Mediterranean 
coast,  transiting  by  the  Tyrrhenian  Sea  in  winter  and  bypassing  it  in  summer.  The  entrance  of 
the  LIW  into  the  Tyrrhenian  Sea  appears  controlled  by  the  flow  reversal  occurring  between 
winter  and  summer  and  which  was  not  observed  in  experiment  MED-1. 

In  the  EMED  the  LIW  outflow  through  the  Strait  of  Sicily  is  formed,  similarly  to 
experiment  MED-1.  by  two  currents  coming  from  the  Northern  and  Southern  Ionian,  but  no 
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■veakening  of  the  northern  current  occurs  in  summer.  The  cenrral  Ionian  and  the  Cretan 
Passage  are  permanently  occupied  by  two  andcycionic  gyres.  In  winter  the  LIW  recirculates 
into  the  gyre  located  south  of  Crete;  in  summer  the  two  gyre  become  interconnected  and  the 
LIW  exit  the  Levantine  Sea  without  any  recirculation.  In  the  nonhem  Ionian  and  in  the 
southern  .Adriatic,  the  flow  is  cyclonic  over  the  entire  year.  In  the  Levannne  Sea  the  cyclonic 
gyre  observed  in  the  previous  experiment  occupies  the  center  of  the  Basin,  is  less  intense  and 
east  of  Rhodes  an  anticyclone,  similar  to  the  .Anaximander  gyre  of  Fig.  a.  is  developed. 

The  circulation  at  1200  m  depth.  Fig.  20.  is  more  stable  throughout  the  year,  and 
shows  marked  differences  with  the  results  of  experiment  .VIED- 1  in  the  WMED.  where  the 
circulation  in  the  northwestern  Basin  does  not  reverse  and  remains  cycionic.  In  the 
Tyrrhenian  Sea,  dominated  by  a  double  gyre.  In  the  EMED  the  circulation  is  still 
characterized  by  the  cyclonic  circulation  encircling  the  whole  Ionian  Sea  (fed  by  water  of 
.Adriatic  origin)  and  by  the  cyclonic  flow  in  the  Levantine  Sea. 

The  distribution  of  scalar  properties  and  the  transformation  of  water  masses. 

We  describe  now  the  temperature  and  salinity  horizontal  distribution,  obtained  in  the 
last  year  of  integration  of  the  experiment  MED-4.  Fig.  21  shows  the  surface  temperature 
distribution  for  the  months  of  February  and  August.  The  WMED  winter  surface  temperature 
ranges  from  minimum  values  slightly  lower  than  12°  C,  approximately  in  the  center  of  the 
nonhem  Balearic  cyclonic  cell,  to  maximum  values  slightly  higher  than  14°  C  (western 
■Alboran  Sea).  Such  values  appear  to  be  about  1°C  colder  with  respect  to  the  monthly 
climatologies  published  by  ENEA  (1990) .  The  tendency  of  the  model  to  produce  temperature 
values  lower  than  the  observed  data  in  the  WMED  is  more  evident  in  summer  where  the 
difference  between  the  model  results  and  the  abovementioned  climatologies  is  about  2-3°  C. 
However,  the  horizontal  surface  temperature  distribution  pattern  in  summer  looks 
qualitatively  similar  to  the  observations.  The  minimum  values  are  located  approximately  in 
the  center  of  the  cyclonic  area  of  the  northern  Balearic  Basin. 


A-94 


In  the  EMED  the  temperature  surface  values  are  closer  to  the  observed  data  {A nan. 
19~~.  Ovchinnikov.  198A).  In  winter  the  temperature  increases  west  to  east  from  14®  :o  17° 
C.  In  the  center  of  the  Rhodes  gyre  the  water  temperature  is  about  15°  C.  a  value 
corresponding  to  the  LIW  source  water  temperature.  In  summer,  in  the  southern  Ionian  Sea 
and  in  the  Cretan  Passage  (the  area  of  the  Ionian-.A.tlantic  Streaml  the  isotherms  are  zonally 
oriented  m.aridng  the  \LAW  flow;  the  temperamre  decreases  nonhward  from  27°  to  21°  C.  In 
the  Levantine  Sea.  on  the  contrary,  the  isotherms  tend  to  assume  a  meridional  orientation. 

The  surface  salinity  fields  for  the  same  months  are  shown  in  Fig.  22.  In  the  surface 
field  the  MAW  inflow  in  the  southern  Balearic  Basin  is  evident  throughout  the  year  as  a 
progressive  increase  of  the  salinity  value  from  west  to  east.  However,  from  winter  to 
summer  the  MAW  inflow  influences  more  the  interior  of  the  Balearic  Basin,  as  indicated  by 
the  nonhward  displacement  of  the  isohalines,  which  tend  to  group  in  the  area  between  the 
Balearic  Islands  and  Sardinia,  giving  hints  of  the  formation  of  a  structure  similar  to  the  nonh 
Balearic  front.  To  the  nonh  the  salinity  distribution  reflects  the  cyclonic  circulation  already 
described. 

As  in  the  case  of  temperature,  the  surface  salinity  values  of  the  WMED  generated  by  the 
model  do  not  match  the  observed  data  very  well;  the  model  values  are  0.15  -  0.20  °/oo 
lower,  although  the  distribution  is  qualitatively  similar  to  observed  fields.  Also  in  the  EMED 
the  surface  salinity  values  increase  from  west  to  east  reaching  the  maximum  values  (  >  39.20 
°/oo  in  winter  and  >  39.60  °/oo  in  summer)  in  the  Levantine  Sea.  The  isohalines  have  a 
predominantly  meridional  orientation.  The  transformation  of  the  MAW  salinity  value  under 
the  effect  of  the  surface  forcing  is,  obviously,  more  rapid  in  summer  than  in  winter,  for 
example  the  salinity  of  the  MAW  entering  the  Levantine  Sea  in  winter  is  about  38.50  o/oo, 
while  in  summer  they  have  a  value  of  about  39.00  %o- 

The  temperature  and  salinity  horizontal  distribution  at  300  m  depth  in  February  and 
.August  (Fig.  23  and  24  respectively)  reflects  the  circulation  patterns  described  above;  in 
panicular  the  temperature  field  shows  clearly  the  cyclonic  circulation  in  the  nonhem  Balearic 
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Sea  .'.nd  ihe  LIW  ouitlow  in  the  Cretan  Passage,  whereas  the  salinity  fields  rerlects  the 
-nac;.  clonic  circulation  m  the  Ionian  Sea.  Fig.  24  also  shows  that,  at  300  m  depth  the  salinity 
aiues  in  the  WMED  are  lower  than  reality,  in  parricular  in  the  southern  Baleanc  Sea. 

Fig. 25  shows  the  location  of  several  cross  sections  obtained  by  bivanate  interpolation 
:'rom  :ne  curv'liinear  grid  and  descnbed  below  in  order  to  illustrate  the  vertical  dismbution  of 
■vater  masses.  Fig.  26.  shows  the  February  and  .August  salinity  distribution  along  the  cross- 
>ecr.on  .A.  The  discrepancy  between  the  observed  salinity  data  and  those  created  by  the  model 
;s  panicularly  evident  at  depths  between  300  and  the  600  m.  the  LIW  layer,  wnc.^e.  in 
comparison  with  the  nearby  cross-sections  by  Katz  (1972)  and  ENEA  ‘  1990).  the  salinity  is 
fresher  by  about  0.20  -  0.40  A'oo.  The  salinity  decrease  at  interm.ediate  depths  occurs 
progressively  throughout  the  integration  time  of  every  numerical  e.xpenmient  performed  and  it 
does  not  seem  likely  that  such  decrease  is  caused  by  the  surface  flux  forcings.  It  is  possible, 
however,  that  this  freshening  depends  upon  the  anticycionic  circulation  developed  in  the 
southern  Balearic  Basin  at  these  depths  and  may  be  caused  by  excessive  vertical  mixing  with 
water  from  the  upper  layers. 

The  salinity  vertical  distribution  in  Fig.  26  shows  also  that,  despite  such  discrepancies, 
the  model  is  successful  in  reproducing  and  maintaining  one  of  the  most  imponant  and  typical 
features  of  the  Mediterranean  Sea  thermohaline  structure,  the  well  known  summer 
development  of  a  subsurface  salinity  minimum  marking  the  path  of  the  .VLAW  flowing  into 
the  Basin  which  is  modified  by  the  surface  warming  and  fresh-water  losses.  The  salinity 
minimum  is  clearly  defined  in  the  WMED  at  depths  of  about  30  -  50  m  and  it  extends 
throughout  the  whole  Basin,  reaching  the  Strait  of  Sicily  and  developing  also  in  the  EMED  at 
a  depth  of  about  50  m  (.cross-section  C,  Fig.27)  as  far  as  the  Cretan  Passage  and  vanishing 
into  the  Levantine  Sea.  In  both  basin  the  salinity  minimum  is  not  developed  in  winter  and 
appears  only  during  the  summer  season  throughout  the  entire  integration  time  of  the  model. 

In  Fig.  27  we  show  also  another  imponant  feature  that  the  model  reproduces,  the 
presence  of  a  LIW  core  with  a  salinity  value  greater  than  38.80  T'oo'  which  occupies  almost 
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ail  of  the  Levantine  Sea  and  extends  panially  into  the  Ionian  Sea  through  the  Cretan  Passage. 
The  development  of  a  LIW  core  led  us  to  investigate  the  location  of  its  source.  The 
northem/central  sector  of  the  Levantine  Sea  for  Rhodes-Cyprus  area)  is  a  well  investigated 
LT-V  source  area  (see.  for  instance  Ovchinnikov.  1984)  so  that  in  Fig.  2S  we  show  the  winter 
distribution  of  potential  temperature,  salinity  and  sigma-Q  along  the  meridional  cross-section 
C.  encompassing  the  whole  Levanrine  Sea  and  crossing  the  Rhodes-Cyprus  region.  The 
temperature  distnbution  shows  the  doming  of  the  isotherms  in  the  center  of  the  section 
determined  by  the  Rhodes  cyclone  circuladon  iAnaii.  1 984)  and,  overlying  the  dome,  a  layer 
of  water  with  a  thickness  of  about  300  m  having  homogeneous  characteristics:  potential 
temperature  between  15.25  and  15.50°  C.  salinity  greater  than  39.05  and  a  sigma-0 
value  comprised  between  29.00  and  29.05.  Such  values  are  In  close  agreement  with  the 
values  observed  during  LIW  formation  in  the  Rhodes  gyre  {Ovchinnikov ,  1984).  Tne  general 
features  of  the  thermohaline  structure  illustrated  in  the  section  indicate  that  two  factors 
concurrently  operate  to  break  the  surface  water  stratification  and  to  form  a  relatively  deep 
water  layer  with  homogeneous  properties:  one  is  the  divergence  bringing  dense  water  closer 
to  the  surface,  the  second  are  the  thermal  and  the  evaporative  losses,  which  determine  the 
winter  buoyancy  fluxes  resulting  in  the  increase  of  surface  water  density  and  convective 
mixing. 

Lasdy,  we  illustrate  the  model  results  relative  to  the  formation  of  the  deep  water  masses 
in  the  Adriatic  Sea  and  in  the  northern  Balearic  Basin  source  areas.  In  Figs.  29  and  30  we 
show  the  potential  temperature  and  sigma-9  winter  distribution  along  section  E,  crossing  the 
whole  Adriatic  Basin  from  the  Gulf  of  Venice  to  the  Strait  of  Otranto.  The  distribution  of 
such  properties  is  plotted  at  a  10  days  interval  for  the  period  February  1  -  February  30  of  the 
climatological  year.  The  potential  temperature  plots  show  the  shallow  nonhem  Adriatic 
occupied  by  a  cold  water  mass  whose  temperature  increases  to  the  south  from  -1 1°  to  13°  C. 
The  shapes  of  the  isotherms  indicate  the  flow  of  bottom  water  along  the  continental  slope. 
The  dme  evolution  of  the  potential  temperature  distnbution  is  clear  evidence  of  the  shelf  deep 
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•vaier  formation  process,  i.e.  the  flow  of  cold  waters  along  the  continental  shelf  iArregiani 
.ir.d  Saiusti.  1987)  sliding  into  the  southern  deeper  pan  of  the  Basin.  In  a  time  span  of  30 
days  the  greater  pan  of  the  southern  .Adriatic  water  mass  is  cooled  by  about  0.25°  C.  .Also 
the  salinity  distnbution  (not  shown  here)  reflects  the  same  patterns.  The  sigma-0  plots  for  the 
same  period  of  time,  conrirm  the  e.xistence  of  this  process,  but  also  give  indications  of  a 
second  deep-water  formation  process  acting  in  the  deeper  southern  .Adriatic  Sea.  In  fact  the 
piot  for  the  end  of  Februar>’  (lower  panel  of  Fig.  29)  shows  the  formation  of  a  surface  layer 
whose  sigma-0  value  (29.25)  equals  the  value  of  the  underlying  water  mass,  suggesting  the 
potential  for  intensive  convective  mixing,  developing  according  to  a  MEDOC  like 
mechanism,  as  observed  by  Ovchinnikov  ei.  al.  (1987)  in  the  southern  .Adriatic  Sea.  The 
newly  formed  dense  water  mass  flows  over  the  Otranto  sill  and  spreads  at  depth  to  the  whole 
E.\LED.  as  shown  in  fig  20. 

Both  processes  are  known  to  be  dependent  on  air  sea  interactions,  panicularly  wind 
related  winter  cooling  and  evaporation  which  determine  the  surface  density  increase  and 
initiate  the  ovenuming  and  decay  of  the  water  column  stratification  (Henderschott  and 
Rizzoli,  1976).  In  the  southern  Adriatic,  a  significant  preconditioning  role  is  also  played  by 
the  permanent  cyclonic  circulation  observed  there  (Fig.  18)  that  brings  dense  water  closer  to 
the  surface  as  described  previously  in  the  case  of  LIW  formation. 

On  the  other  hand,  we  were  not  able  to  find  any  significant  evidence  of  intensive  deep 
mixing  of  the  water  column  in  the  nonhern  Balearic  Sea.  Fig.  31  shows  the  salinity 
distribution  along  the  cross  section  B  in  .March.  It  can  be  noted  that  only  a  small  and  limited 
break  in  the  stratification,  extending  from  the  surface  to  about  250  -300  m  depth  occurs  in 
this  period,  during  which  deep  water  formation  processes  were  observed;  and  the  general 
distribution  of  salinity,  density  and  temperature  fields  (not  shown  here)  differ  from  those 
published  by  the  MEDOC  group  (1969),  Siommel  (1972)  and  Leaman  and  Schott  (1991). 
The  different  behaviour  of  the  model  in  reproducing  the  deep  water  formation  in  the  .Adriatic 
Sea  and  in  the  Gulf  of  Lions  is  difficult  to  explain.  The  adjustment  toward  fresher  salinity 
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'-aiues  ;n  :he  WMED.  and  the  resulting  surface  water  lighter  than  observed  are  probably 
pia>ing  an  inhibiting  role.  Moreover,  the  coarse  grid  resolution  is  probably  inadequate  to 
capture  tne  dynamic  of  the  process  and  the  use  of  climatological  flu.xes  to  force  the  model 
eliminates  the  strong  wind  (Mistral)  bursts  that  start  the  deep  water  formation  {Bunker. 
19~2).  but  these  two  factor  are  present  also  in  the  .Adriatic  Sea  where  the  deep  water  is 
formed.  We  leave  this  apparent  conaadiction  as  an  open  question  for  funher  studies. 


S:  DISCUSSION  AND  CONCLUSIONS 

The  differences  arising  from  the  comparison  of  the  results  obtained  from  experment 
NLED-1  and  MED-4  point  out  the  relative  importance  of  the  seasonal  wind  stress  and  surr'ace 
heat  and  salinity  fluxes  in  modeling  the  general  circulation  of  the  Mediterranean  Sea.  In  the 
WAtED,  both  experiments  provide  a  general  picture  of  the  surface  circulation  in  the  Baleanc 
Basin  characterized  by  a  cyclonic  circulation  in  the  north  and  by  the  N^LAW  inflow  in  the  south 
which,  in  summer,  generates  mesoscale  anticycionic  eddies  originated  by  instabilities  of  the 
.Almeria-Oran  ffont/Algerian  current  system,  and  extends  its  influence  farther  nonh  as  defined 
by  the  position  of  the  north  Balearic  front.  The  MAW  flow  in  the  WMED  (as  well  as  in  the 
EMED)  is  characterized  by  the  typical  summer  development  of  a  salinity  subsurface 
minimum.  The  substantial  similarity  between  the  two  experiment,  as  well  as  the  results  of 
other  modeling  effons  {Heburn,  1987;  Pinardi  and  Navarra,  1991),  indicates  that,  in  general, 
the  wind-stress  seasonality  is  a  very  imponant  forcing  factor.  However,  we  have  shown  that 
the  temporal  and  spatial  variability  in  the  surface  salinity  flux  and  the  effect  of  fresh-water 
runoff  arc  necessary  to  the  maintenance  of  a  realistic  horizontal  density  gradient  which  drives 
the  cyclonic  circulation  in  the  northern  Balearic  Basin.  Moreover,  the  seasonal  surface 
buoyancy  flux  seems  to  have  an  enhancing  effect  on  the  summer  generation  of  the 
anticycionic  mesoscale  eddies  in  the  south  and  strengthens  the  seasonal  reversal  of  the 
circulation  in  the  Tyrrhenian  Sea.  which  is  extended  in  depth  from  the  surface  to  the 
intermediate  layers.  This  imponant  seasonal  change  produces  a  variation  in  the  LIW 
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circuiadon  patrern.  as  the  current  reversal  in  the  Tyrrhenian  eliminates,  or  greatly  reduces, 
the  Lrvv  flow  into  the  Basin.  Thus  it  seems  that  throughout  the  year  the  LIW  circulation 
osculates  between  a  winter  configuration  similar  to  the  hypothesis  put  forward  by  Millot 
1 1987b),  involving  the  recirculation  in  the  Tvirhenian  Sea  and  a  summer  configuration  more 
in  agreement  with  the  interpretation  given  by  Katz  (1972)  involving  the  direct  flow  of  the 
LFW  into  the  Balearic  Basin  along  the  western  coast  of  Sardinia.  This  is  obviously  a  topic 
requinng  further  obser-'anonal  and  theoretical  studies. 

In  the  Ionian  sea.  the  wind-stress  and  the  surface  buoyancy  fluxes  have  a  balancing 
effect  on  the  seasonal  variability  of  the  loman-.Atlantic  stream,  as  the  introduction  of  the 
surt'ace  fluxes  in  expenment  MED-4  reduces  the  large  nonhward  shift  of  the  stream 
occurring  from  winter  to  summer  and  the  development  of  the  broad  anticyclonic  gyre  in  the 
southern  and  central  Ionian  Sea.  In  the  northern  Ionian  Sea.  on  the  contrary,  both 
experiments  gave  evidence  of  a  weakening,  in  the  transition  from  winter  to  summer 
conditions,  of  the  observed  cyclonic  circulation.  Thus  the  main  forcing  mechanism 
responsible  for  the  weakening  is  the  wind-stress. 

An  assessment  of  the  leading  forcing  in  the  Levantine  Sea  is  more  problematic.  The 
surface  circulation  is  constituted  of  sub-regional  gyres  of  both  sign  interconnected  by  the 
Mid-Mediterranean  jet  and  the  Asia  Minor  current.  Both  expenments  reproduced  the  Rhodes 
gyre,  the  bifurcation  of  the  Mid- Mediterranean  jet  east  and  west  of  Cyprus  and  the  Asia 
Minor  current.  However,  the  Cretan  cyclonic  gyre,  the  Mersa-Matruh  and  Shikmona 
ancicycionic  gyres  appear  bener  simulated  in  experiment  MED- 1  than  in  experiment  .MED-4. 
Such  features  appear  very  clearly  in  the  geostrophic  computations  derived  from  the  quasi- 
synoprical  data  collected  during  the  POEM  project  {Oszoy  et  aL,  1989).  The  Mersa-Matruh 
gyre,  in  parricular,  constitutes  the  stronger  signal  observed.  This  evident  discrepancy 
between  models  simulations  and  observations  indicates  a  possible  underestimation  of  the 
surface  heat  and  salinity  fluxes  in  the  Levantine  Sea  area,  whose  effects  on  the  circulation  are 
probably  counterbalanced  by  the  wind-stress  influence. 
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Ac  incermediate  depths  the  LIW  outflow  through  the  Strait  of  Sicily  is  reproduced  well 
by  the  model.  Their  path  across  the  EMED  undergoes  seasonal  changes,  panicuiarly  in  the 
connecting  flow  between  the  Levantine  and  the  Ionian  Seas,  across  the  Cretan  Passage, 
which  appears  reduced  in  winter.  The  LIW  current  in  the  Ionian  is  entrained  in  the 
ancicyclonic  cell  located  in  the  center  of  the  Basin  and  before  entering  the  Strait  of  Sicily 
rriLxes  funherly  with  water  coming  from  the  northern  Ionian  Sea. 

The  '.vater  mass  formation  processes  occurring  in  the  Mediterranean  Basin  are  strongly 
related  to  both  the  honzontal  circulation  and  the  surface  flu-xes.  We  have  shown  that  the  LIW 
and  the  deep  water  formation  areas  are  coincident  with  honzontal  cyclonic  flows  which  cau.se 
causing  doming  of  the  underlying  water  and  transpon  of  dense  water  to  the  suriace.  On  this 
structure,  the  surface  fluxes  create  thermal  and  evaporative  losses  which  lead  to  convective 
venical  mixing.  The  model  results  show  that  this  mechanism  is  active  in  the  Rhodes  gyre 
region.  .An  analogous,  but  more  intense,  situation  responsible  for  the  EMED  deep  water 
formation  occurs  in  the  southern  Adriatic  Sea.  Intensive  venical  mixing  over  the  deeper  pan 
of  the  Basin  and  the  shelf  deep  water  formation  are  reproduced  remarkably  well  by  the 
model.  However,  the  model  did  not  produce  realistically  deep  water  formation  in  the  WMED 
where  in  the  nonhem  Balearic  Sea  the  model  produces  only  a  small  and  spatially  limited 
homogeneous  layer  extending  down  to  200  m.  We  attribute  this  to  the  generalized  adjustment 
of  the  WMED  toward  lower  salinity  values  as  well  as,  probably,  to  an  insufficient  horizontal 
resolution,  but  a  precise  assessment  of  the  reason  for  the  lack  of  deep  water  formation  in  the 
Gulf  of  Lions  will  require  funher  studies. 

Finally,  we  note  that  the  use  of  climatological  surface  forcing  fields  eliminates 
completely  from  the  model  performance  any  influence  on  the  circulation  due  to  the 
interannual  variability  which  in  several  areas  of  the  Mediterranean  Sea,  such  as  the  Adriatic, 
the  Tyrrhenian  Sea,  the  LIW  and  the  deep  water  formation  areas,  can  play  a  very  imponant 
and  significant  role.  The  use  of  synoptic  surface  forcing  fields  together  with  a  higher 
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resolution  model  will  surely  provide  a  better  simulation  and  understanding  of  the  circulation 
processes  in  the  Mediterranean  Sea. 
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Tdoie  1 :  Characteristics  of  the  numerical  experiments. 
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FIGURE  CAPTIONS 


Fia.  1;  The  Mediterranean  Sea.  Upper  panel;  Geographic  features  (adapted  from  Hopkins. 
197S  )  lower  panei;  bathymetnc  map  (source:  U.S.  Navy  bathymetric  data-bae).  Contour 
interv'al  is  250  m. 

Fig.2;  Schematic  of  the  Mediterranean  Sea  water  masses  venical  distribution  (adapted  from 
Hopkins.  1985') 

Fig. 3;  Schematic  of  the  WMED  circulation.  Upper  panel:  MAW  circulation,  lower  panei: 
LI\V  circuladon  (adapted  from  Millot.  1987a). 

Fig. 4:  Schematic  of  the  ENIED  circulation  (adapted  from  Robinson  et  al..  1991). 

Fig. 5;  The  model  gnd. 

Fig. 6;  The  model  bathymetry'.  Contour  interval  is  250  m. 

Fig. 7:  .Annually  averaged  hydrological  climatologies  (model  initial  conditions).  Upper  panel: 
surface  temperature.  Contour  interval  is  1*^  C,  lower  panel;  surface  salinity.  Contour  interv'al 
is  O.lWoo- 

Fig.S:  Total  surface  heat  flux.  Upper  panel;  February,  lower  panel;  August.  Contour  interv-al 
is  10  W*m-2.  Not  all  the  vectors  are  plotted 

Fig.9:  Wind-Stress.  Upper  panel:  February,  lower  panel  August.  Units  are  dyne*cm’-. 

Fig.  10:  Surface  fresh-water  flux.  Upper  panei:  February,  lower  panel:  August.  Contour 
interval  is  10  mm. 

Fig.l  1;  Experiment  MED-1;  Surface  elevation.  Upper  panel:  February,  lower  panel:  .August. 
Contour  interval  is  2  cm. 

Fig.  12:  E.xperiment  MED-1;  Velocity  at  50  m  depth.  Upper  panehFebruary,  lower  panel: 
August.  Units  are  m'^s'^  In  this  picture,  as  well  as  the  other  velocity  field  picture,  the  vectors 
in  the  Gibraltar  area  are  not  plotted  for  clarity  sake. 

Fig.l3:E.xperiment  MED-1;  Velocity  at  300  m  depth.  Upper  panehFebruary,  lower  panel: 
August.  Units  are  m*s‘^ 

Fig.  14:  E.xperiment  MED-1;  Velocity  at  1200  m  depth.  Upper  panehFebruary,  lower  panel: 
August.  Units  are  m*s‘^ 

Fig.  15:  Experiment  MED-2:  Surface  elevation.  Upper  panel:  February,  lower  panel;  .August. 
Contour  interval  is  2  cm. 

Fig.  16.  Temperature  (left)  and  salinity  (right)  Basin  average  time  historv  for  experiment 
.VIED-4. 


Fig.  17:  E.xperiment  .VIED-4:  Surface  elevation.  Upper  panel:  February,  lower  panel;  .August. 
Units  are  m*s'F 
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Fig. IS;  E.xperiment  MED-4;  Velocity  a:  50  m  depth.  Upper  panel;  February-.  lower  panel; 
•August.  Units  are  m^s'F 

Fig.  19;  E.xperiment  .VLED-4;  Velocity  at  300  m  depth.  Upper  panel;  February',  lower  panel; 
.August.  Units  are  m*s‘l. 

Fig. 20;  E.xperiment  MED-4;  Velocitv'  at  1200  m  depth.  Upper  panel;  February,  lower  panel; 
•August.  Units  are  m^s'^. 

Flg.21;  Locatio'^  of  the  cross-section  described  in  the  te.xt. 

Fig. 22;  E.xperiment  MED-4;  Surface  temperature  distribution.  Upper  panel;  February 
(contour  interval  is  0.5°  C),  lower  panel;  .August  (contour  interval  is  1.0°  CA 

Fig. 23;  Experiment  MED-4;  Surface  salinity  distribution.  Upper  panel;  February  i contour 
interval  is  0.  l°/oo),  lower  panel;  .August  (contour  interval  is  0.2°/ooA 

Fig. 24:  E.xperiment  MED-4;  Temperature  distribution  at  300  m  depth.  Upper  panel; 
February,  lower  panel;  .August,  contour  interval  is  0.2°  C. 

Fig. 25:  E.xperiment  MED-4:  Salinity  distribution  at  300  m  depth.  Upper  panel;  February, 
lower  panel:  August.  Contour  interval  is  0.05  °/oo- 

Fig.26;  Experiment  MED-4:  Salinity  distribution  along  cross-section  A  in  February  (upper 
panel)  and  August  (lower  panel).  Contour  interval  is  0. 1  °/oo. 

Fig.27:  E.xperiment  MED-4;  Salinity  distribution  along  cross-section  D  in  February  (upper 
panel)  and  August  (lower  panel).  Contour  interval  is  0.1  °/oo. 

Fig.28:  E.xperiment  MED-4:  Cross-section  C.  Temperature  (upper),  salinity  (middle)  and 

sigma-0  distribution  in  February.  Contour  intervals  are  0.25°  C,  0.05  °/oo  and  0.05  sigma-0 
units. 

Fig. 29;  Experiment  MED-4:  Cross-section  E.  Temperature  distribution  on  February  10 
(upper),  20  (middle)  and  30  (lower)  of  the  climatological  year.  Contour  interval  is  0.25°C. 

Fig.30:  Experiment  MED-4:  Cross-section  E.  Sigma-0  distribution  on  February  10  (upper), 
20  (middle)  and  30  (lower)  of  the  climatological  year.  Contour  interval  is  0.05  sigma-0  units. 

Fig.31;  Experiment  MED-4:  Cross-section  B:  salinity  distribution  in  March.  Contour  interval 
is  0.05  °/oo. 
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THE  NRL  MEDITERRANEAN  SEA  MODEL 
RECENT  RESULTS 

by 

G.W.  Hebum 

Naval  Research  Laboratory 
NRL  Code  322 

Stennis  Space  Center,  MS  39529-5004 


INTRODUCTION 

The  circulation  in  the  Mediterranean  Sea  is  extremely  complex,  consisting  of  numerous  eddies  and 
current  meanders,  particularly  along  the  north  African  coast ,  in  the  Algerian  basin  and  throughout 
the  Ionian  and  Levantine  basins.  Satellite  IR  and  CZCS  imagery  as  well  as  recent  hydrographic 
surveys  (for  example  the  POEM  program)  reveal  a  picture  of  the  general  circulation  in  the  western 
and  eastern  Mediterranean  consisting  of  sub-basin-scale  gyres  and  interconnecting  jets.  There  are 
many  physical  factors  which  exert  an  influence  on  the  dynamics  of  the  circulation  of  the 
Mediterranean  Sea;  wind  stress,  hydraulic  controlled  inflow/outflow  through  the  Straits  of  Gibraltar 
and  Sicily,  inteimediate  and  deep  water  formation,  theimohaline  forced  circulation  and  topography. 

Various  versions  of  the  NRL  primitive  equations  ocean  circulation  model  are  used  to  systematically 
study  the  affects  of  these  various  forcing  mechanisms,  individually  and  in  combination,  on  the 
circulation  dynamics  of  the  Mediterranean  Sea.  While  the  general  results  from  various  case  studies 
involving  wind,  hydraulic,  theimohaline  and  topographic  forcing  suggest  that  no  one  mechanism 
dominates  the  complex  Mediterranean  circulation  and  that  non-linear  interactions  between  the 
directly  forced  motions  and  internal  flow  instabilities  play  an  important  role,  comparisons  of  model 
results  from  experiments  using  the  ECMWF  winds  for  the  period  1981  to  1989  with  observational 
data,  both  satellite  IR  and  hydrographic,  suggest  that  the  winds  play  a  major  role  in  the  evolution  of 
the  sub-basin-scale  gyres  seen  in  the  Mediterranean  Sea. 

MODEL 

The  NRL  multi-layer  primitive  equations  ocean  circulation  model  is  a  derivative  of  the  Hurlburt  and 
Thompson  ( 1980)  Gulf  of  Mexico  model.  For  a  full  description  of  the  basic  model  refer  to  Wallcraft 
( 1991 ).  The  version  used  for  most  of  the  recent  dynamical  studies  in  the  Mediterranean  Sea  is  a 
three  layer,  finite  depth,  hydrodynamic,  spherical  coordinate  version  of  the  basic  model.  The  model 
domain  is  presented  in  figure  1.  The  lateral  boundaries  for  the  model  are  taken  at  the  200m  isobath. 
The  initial  layer  thicknesses  are  take  to  be  200m  for  the  upper  layer  representing  the  Modified 
Adantic  Water  (MAW),  400m  for  the  second  layer  representing  the  Levantine  Inteimediate  Water 
(LIW)  and  the  remainder  in  the  third  layer  to  the  bottom  topography  representing  the  Mediterranean 
Deep  Water  (MDW).  The  inflow  boundary  condition  was  specified  as  a  steady  inflow  of  1.2sv 
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through  the  Strait  of  Gibraltar  in  the  upper  layer  with  conpensating  outflow  in  the  second  layer.  The 
formation  of  LIW  is  parameterized  as  a  steady  detrainment  for  the  upper  layer  to  the  second  layer  in 
the  eastern  Mediterranean  basin. 


Figure  1 :  Mediterranean  Sea  Model  Topography 

The  model  was  integrated  to  statical  equilibrium  using  the  Hellerman  and  Rosenstein  monthly  mean 
wind  stress  climatology  and  then  integrated  for  a  nine  year  period  with  a  hybred  wind  stress  data  set 
using  the  annual  mean  from  the  Hellerman  and  Rosenstein  climatology  and  the  monthly  mean 
annomaUies  from  the  ECMWF  winds  for  the  period  from  1981  to  1989. 

RESULTS 

The  results  from  the  simulations  using  the  hybred  wind  data  set  have  been  comparied  to  satellite  DR 
images  for  various  time  periods  during  the  1980s.  The  disappearance  and  reformation  of  the  eastern 
Alboran  gyre  during  July  to  August  1982  (Hebum  and  LaViolette,  1990)  is  well  represented  in  the 
model  simulation.  The  mean  annual  cycle  of  the  transport  of  LIW  through  the  Strait  of  Sicily  is 
replicated  in  the  model  simulations  using  the  climatological  wind  forcing  as  well  as  in  the  simlations 
using  the  hybred  wind  forcing  data  set.  The  simulations  using  the  hybred  data  set  also  show 
considerable  interannual  varibility. 

The  complex  series  of  sub-basin  scales  gyres  and  meandering  currents  observed  in  the  eastern 
Mediterranean  Sea  are  shown  in  the  model  simulations  as  well.  These  features  also  exibit 
considerable  seasonal  and  interannual  variability.  The  gyres  in  the  Ionian  basin  appear  to  be 
primarily  wind  forced  with  a  definite  annual  cycle. 
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The  modeling  effort  is  the  one  undertaken  by  Malanotte-Rizzoli  and  collabora¬ 

tors  in  which  three  approaches  have  been  tised  simxiltaneously.  In  the  first  study,  Tziper- 
man  and  Malanotte-Rizzoli  (1991)  use  a  simple  inverse  model  to  calculate  the  horizontal 
circulation  of  the  upper  SOO  m  of  the  entire  Mediterranean  Sea  for  a  seasonal  climatological 
hydrographic  dataset.  The  inverse  results  were  used  to  address  three  main  aspects  of  the 
general  Mediterranean  circulation,  namely  the  seasonality  of  the  surface  circulation;  the 
shape  of  the  general  circulation  of  the  Eastern  basin  and  the  shape  of  the  deep  circulation 
of  the  entire  Mediterranean.  Figure  1  shows  the  winter  circulation  (upper  panel)  and  the 
summer  one  (lower  panel)  in  the  sxirface  layer,  at  a  depth  of  10  m.  A  fairly  steady  yearly 
pattern  of  the  surface  circulation  in  the  Western  basin  emerges  from  an  examination  of 
Figure  11.  The  seasonal  signal  is  somewhat  stronger  in  the  Eastern  basin  and  especially 
in  the  Ionian  Sea  where  the  jet  of  Atlantic  water  entering  the  Sicily  Straits  shoots  south¬ 
eastward  in  winter  (Figure  11a)  while  forming  an  intense  meander  that  protrudes  into 
the  Northern  Ionian  in  stimmer  (Figure  lib).  Also,  Figure  1 .  clearly  shows  the  Eastern 
Mediterranean  circulation  to  be  constituted  by  a  succession  of  subbasin  scale  gyres. 

In  a  parallel  study,  Malanotte-Rizzoli  and  Bergamasco  (1991)  carried  out  an  intensive 
series  of  numerical  experiments  investigating  the  wind-driven  and  thermohaline  compo¬ 
nents  of  the  Eastern  Mediterranean  circulation  and  the  relative  importance  of  the  different 
forcing  functions  under  climatological  conditioxis  in  the  different  regions  of  the  Eastern 
basin.  They  used  a  linearized,  primitive  equation  level  model  with  coarse  resolution  but 
with  acti'.  e  thermodynamics,  and  realistic  topography  with  the  Sicily  Straits  open.  The 
model  is  forced  by  the  monthly  climatological  means  of:  wind  stress;  heat  fluxes  (incom¬ 
ing  solar  radiation,  backscatter,  latent  and  sensible  heat);  and  evaporation  minus  pre¬ 
cipitation  fluxes.  The  inflow  imposed  at  the  Sicily  Straits  is  derived  from  real  data  and 
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tba  hydrographic  temperature  and  salinity  dataset  is  the  same  used  by  Tziperman  and 
Malanotte-Rizzoli  (1991). 

The  numerical  results  indicate  that  the  inflow  at  Sicily  is  the  dominant  driving  mech¬ 
anism  in  the  Ionian  Sea,  with  a  strong  seasonal  modTiIation  imposed  by  the  change  in 
polarity  of  the  wind-stress  curl  from  winter  to  summer  over  the  Ionian.  However,  the  ther¬ 
mohaline  fluxes  emerge  to  be  the  dominant  driving  mechanism  for  the  circulation  in  the 
Eastern  Levantine,  where  the  seasonal  signal  is  greatly  attenuated.  Figure  2  shows  the 
total  transport  streamfunction  in  the  Central  Experiment,  the  one  both  wind  and  ther¬ 
mally  forced  in  which  all  the  fields  are  evaluated  prognostically,  respectively  for  January 
(upper  panel)  and  September  (lower  panel).  Examination  of  Figure  2  shows  the  circu¬ 
lation  to  be  quasi-steady  in  the  Eastern  Levantine.  In  the  Ionian  Sea,  the  jet  of  Atlantic 
water  entering  through  the  Sicily  Straits  is  deflected  southward  towards  the  African  coast 
in  winter  (upper  panel)  while  forming  a  big,  northward  protruding  meander  in  summer 
(lower  panel).  The  similarity  of  Figure  2  with  the  inverse  results  of  Figure  1  is-  striking 
and  also  shows  the  importance  of  the  barotropic,  forced  compouent  of  the  circulation  in 
the  Eastern  basin.  A  final  result  is  that  many  of  the  persistent  features  found  in  the  model 
circulation  patterns  can  be  related  to  the  strong  control  exerted  by  the  ambient  potential 
vorticity  (// H)  upon  the  circulation 


In  a  final  approach,  the  adjoint  method  of  control  theory  is  used  to  determine  the 
seasonal  circulations  of  the  Eastern  basin  (Bergamasco  ei  al.,  1992).  The  GFDL  model, 
adapted  to  the  Eastern  basin,  has  been  coupled  with  its  adjoint  developed  by  collaborators 
Long  and  Thacker.  The  prognostic  model  is  forced  with  the  climatological  wind  stress 
field  and  the  hydrological  data  set  of  temperature  and  salinity  is  used  to  constrain  the 
model  predicted  fields  as  part  of  the  cost  frmetion  of  the  adjoint.  The  cost  function  also 
includes  terms  penalizing  the  tendencies  of  the  prognostic  fields,  i.e.  horizontal  velocities, 
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-  jeffiperature  and  salinity.  The  GFDL  model  and  its  adjoint  are  used  with  1/4*  coarse 
resolution  but  with  17  vertical  levels,  thus  allowing  for  a  realistic  topography  in  the  Eastern 
basin.  The  Sicily  Straits  are  also  open  with  the  same  forcing  inflow  of  Malanotte-Rizzoli 
and  Bergamasco  (1991).  Thus,  a  seasonal  circulation  of  the  Eastern  Mediterranean  is 
reconstructed  that  is  consistent  both  with  the  model  dynamics  and  the  climatological 
hydrographic  dataset.  A  preliminary  result  in  Figure  3  shows  the  circulation  pattern 
in  the  surface  layer,  1  to  50  m  depth,  for  the  winter  season.  Again,  the  similarity  is 
striking  between  the  winter  circulation  and  the  corresponding  patterns  obtained  through 
the  independent  inverse  method  (Figure  la)  and  the  prognostic  calculations  (Figure  2a). 
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A  “CLOSE  TO  OPTIMAL”  SUB-OPTIMAL  FILTER  FOR  DATA 
ASSIMILATION  INTO  LINEAR  MODELS 

Ziv  Sirkes 


Uiveii  a  System  Model: 

■I'k  = 


where: 

.1';.  is  the  [A'  X  l|  system  state  vector  at  time  tk, 

Ak  is  the  [N  X  A"]  state  trausitiou  matrix  at  time  tk, 

Bk  is  the  [jV  x  L\  forcing  matrix  at  time  U-. 

Uk  is  tlie  [L  X  1]  forcing  vector  at  time 

Gk  is  the  [.V  x  A']  error  transition  matrix  at  time  tk, 

Wk  is  the  [A'  x  1]  vector  of  model  error  (assumed  to  be  Gaussian  with  0  mean  and  Q  covariance 
matrix). 

and  a  measurement  Model: 

Zk  =  Hk-tk  +  Vk  (2) 

where: 

Zk  is  the  [*\/  X  i]  measurement  vector  at  time  tk- 

Hk  is  the  [M  X  A’l  matrix  giving  the  interpolation  from  the  state  vector  to  the  measurement 
[joints. 

I'k  is  the  [.U  X  1]  vector  of  measurement  error  (assumed  to  Ije  Gaussian  with  0  mean  and  R 
covariance  matrix). 

One  can  predict  the  state  estimate  as: 

■i'k,k-i  —  .  +  Bk-ii'k-i  (■!' 
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whorp: 

i'’f  tlie  estimate  of  the  state  vector  at  time  tk.  using  measuremeiils  up  to  time  /,_| 
with  ail  error  covariance: 

Pic.k-i  =  ■  U-iA-i.A-i-'ft_i  + 

(  sing  measurements  l<lalal  at  lime  /;.  one  can  update  the  state  estimate  : 

n-.k-  =  h.K-i  -I-  l\\[  :k  -  )  ’  ■'  I 

’.vliere: 

Ai-  is  I  he  lilter  gain 

and  the  error  covariance  ; 

Pk.k  =  (/  —  I\.k^'h)Pk.k-i(I  —  hkHk)^  +  hkUk^k- 

If  one  wants  to  minimize  the  elements  along  the  diagonal  of  the  error  covariance 
giving  an  optimal  filter,  then  the  lilter  gain  becomes: 

/U,p.  =  Pk.k-iHkitikFk.k.iHl  +  Rk)-^.  (7) 

The  error  covariance  for  the  optimal  filter  can  be  simplified  to: 

Pk.k^^,  =  ( /  -  A\.,p,  Hk )  Pk.k-i  ■  ( ) 

The  two  main  problems  with  this  recursive  procedure  (Kalman  filter)  are: 

1.  The  covariance  matrix  is  very  big  ([.V  x  .V]). 

'2.  Ec|uation  4  requires  i,.'3  matrix  multiplications. 

Equations  8  and  4  can  be  combined  into  one  etjuation.  called  the  "Discrete  Ricatti  Eipia- 
tion": 

Pk^uk  =  Ak(Pk.k-i-Pk.k-iHl([IkPk.k-iHl+Rky' ffl Pk.k-i)  Al  (!)) 
+  Qk  C''[  . 
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l  or  a  lime  iiivariaul  case  i  where  A.  B.  G.  Q.  II.  R  rlo  uol  cliaiige  wiili  liiiiei  the  DHL 
rei-luces  to; 

=  -A  (P...-1  -  P;.A-i//^'  +  p)"‘  //"'P.-.A-i)  .A^’  +  '-'Or;''.  I  lU) 

f’silig  a  iloubling  algorithm  one  can  calculate  P^i+i from  Pja-.j^-i  by  [)erlorniiiig  1  l.'i 
matrix  multiplications.  This  is  more  expensive  than  the  l.b  mult  i[)lica  t  ions  in  erp  1.  liowe\er 
1  he  pr(jiires.sion  in  time  is  exponential. 

If  we  use  a  startegy  of  calculating  Pj^+i  _.a  with  the  iloubling  algorithm,  anil  use  a  i  on- 
'tatnt  A^.  until  the  next  update  of  the  error  covariance  matrix,  we  will  have  to  periorm 
less  matrix  multiplications  over  the  lomr  run  .  while  applying  a  subopt  imal  lilter.  Ibis 
idler  will  be  verv  close  to  optimal,  since  it  will  be  upilateil  ver\'  olten  at  the  betiinninii 
i  t  =  1. 2.  f.  8.  IG. ...  I  when  the  error  covariance  matri.x  is  changing  relatively  last,  and  will  l)e 
(  hangeil  less  and  less  freiiuenl  as  t  ime  pas.ses  (and  as  the  error  covariance  mat  rix  approai  lies 
a  steady  state). 

For  an  exi^eriment  with  2^*^  =  1024  assimilation  cycles  one  would  have  to  periorm  lib 
multiplications  using  this  method  (compared  to  l-bOC  multiplications  for  the  regular  Kalman 
filter). 
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PLENARY  GROUP  DISCUSSION 
David  Dietrich 


The  question  of  baroclinic  pressure  gradient  accuracy  was  declared  to  be  a 
non-issue,  based  on  recent  improvements  of  initialization  procedures  reported  by 
George  Mellor,  and  on  recent  vertical  resolution  sensitivity  results  reported  by 
Charles  Horton.  George  also  noted  that  the  energy  conservation  issue  (i.e.,  that, 
unlike  level  vertical  coordinates,  the  truncation  errors  lead  to  non-conservation  of 
available  potential  plus  kinetic  energy;  is  not  a  major  problem.  Perhaps 
geophysical  modelers  have  previously  put  too  much  emphasis  on  energy 
conservation,  when  overall  accuracy  is  the  ultimate  goal.  Also,  Tom  Rosmond  had 
noted  earlier  that  the  widely  publicized  baroclinic  pressure  gradient  issue  in  sigma 
coordinates  was  not  a  major  issue. 

In  particular,  George  noted  that,  in  the  case  of  exactly  stratified  initial  flow, 
the  erroneous  baroclinic  pressure  gradient  associated  with  truncation  errors  can  be 
compensated  by  a  slight  adjustment  of  the  initial  density  field.  In  particular,  the 
initial  density  field  can  be  adjusted  such  that  the  initial  baroclinic  pressure  gradient 
is  exactly  zero,  as  calculated  in  sigma  coordinates.  To  be  useful,  such  adjustment 
must  be  small  compared  to  the  local  density  fluctuations  that  normally  occur. 
Giarles  Horton  reported  results  which  suggest  that,  with  about  20  vertical  levels 
and  10  km  horizontal  resolution,  and  with  real  Mediterranean  bathymetry  and 
realistic  stratification,  this  adjustment  is  sufficiently  small  (equivalent  to  a 
geostrophic  baroclinic  velocity  change  of  less  than  10  cm/sec).  Presumably,  this 
kind  of  adjustment  could  also  be  made  when  realistic  nonzero  initial  horizontal 
density  gradients  occur. 

One  question  that  remains  is  whether  the  effects  of  the  initial  numerical 
density  redistribution  can  be  "unravelled"  after  a  calculation,  so  that  the  Navy  has  a 
sufficiently  accurate  density  field  for  its  intended  applications.  Because  the 
baroclinic  pressure  gradient  issue  relates  strictly  to  truncation  error  effects,  the 
rapporteur  believes  that  a  systematic  resolution  sensitivity  study  should  be  done  for 
a  small  basin  (i.e.,  with  dimensions  about  500  km)  with  bottom  slopes  and 
horizontal  and  vertical  density  gradients  comparable  to  those  in  the  Mediterranean. 
If  this  is  done,  no  serious  CFD  scientist  should  question  the  applicability  of  sigma 
coordinates  to  ocean  problems,  but.  until  this  is  done,  some  geophysical  modelers 
who  are  not  intimately  involved  with  the  details  of  initialization  using  sigma 
coordinates,  might  be  skeptical.  George  noted  that  it  seemed  reasonable  to  have  a 
z-Ievel  model  also  involved  in  the  Mediterranean  studies.  It  was  also  noted  that  the 
only  model  recognized  by  the  Naw  for  use  by  INO  in  such  studies  is  his  (the 
Princeton;  model.  Sigma  coordinates  are  a  good  choice  when  buoyancy  effects  do 
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not  strongly  dominate  over  boundary  iayer  effects,  such  as  in  shallow  continental 
shelf  regions.  However,  they  are  not  as  accurate  or  robust  as  corresponding  z-level 
models  in  deep  water  regions,  esoeciaily  near  the  edge  of  continental  shelves, 
where  buoyancy  effects  dominate,  in  spue  of  the  excellent  job  being  done  by  the 
Princeton  model  developers  in  improving  its  application  to  such  regions.  In  such 
deep  water  regions,  z-level  coordinates  are  closer  to  solution  adaptive  (density 
contour  following)  grids.  It  was  decided  that  it  is  very  difficult  to  represent  two 
very  small  scale  phenomena,  that  are  also  probably  significant,  using  present 
turbulence  closure  schemes.  These  are  coastal  scales  (tides,  diurnal  wind  effects, 
boundary  currents)  and  "chimneys  ’  of  rapidly  sinking  dense  water.  The  former 
can  probably  be  adequately  addressed  by  using  boundary-fitted  coastal  models, 
which  are  currently  under  development,  but  might  require  very  robust  numerics 
(because,  as  suggested  by  Dietrich's  recent  results,  realistic  modeling  of  the 
interaction  between  boundary  and  interior  flows  might  require  very  high  Reynolds 
number  calculations).  The  latter  require  process  studies  to  determine  how  they 
might  be  addressed  in  the  Mediterranean.  Both  of  these  phenomena  are  considered 
important  in  ocean  modeling.  It  was  decided  that  the  software  functionalities  that 
should  be  included  in  the  ECMOP  system  depend  on  the  expressed  needs  of 
individual  scientists.  New  needs  will  probably  evolve,  such  as  those  specific  to 
boundary-fitted  coastal  models.  Free  surface  effects  are  considered  important 
enough  ^at  they  should  be  addressed.  Although  they  are  not  strongly  coupled  to 
ocean  general  circulations  on  short  time  scales,  George  observed  that  free  surface 
models  have  the  advantage  that  they  are  massively  paraUelizable,  because  the  rigid- 
lid  approximation  leads  to  some  pressure  solver  operations  that  are  paraUelizable 
o,:ly  in  one  dimension,  unless  relatively  slow  SOR-like  solvers  are  used.  (The 
degree  of  SOR  convergence  required  again  raises  the  perhaps  over-emphasized 
energy  conservation  issue.)  The  rapporteur  notes  that  such  SOR  solvers  still 
require  much  fewer  operations  than  the  barotropic  free-surface  mode  requires,  and 
that  a  question  remains;  does  the  interaction  between  free-surface  and  internal 
waves  somehow  limit  the  "long '  time  step  used  by  sigma  coordinate  models  for  the 
internal  modes?  Different  kinds  of  data  that  will  be  useful  for  assimilation  for 
modeling  the  Mediterranean  were  discussed.  Wind  and  heat  flux  data  are  available 
from  atmospheric  models,  but  probably  will  have  to  be  modified  to  higher 
resolution  for  Mediterranean  modeling.  Basic  research  is  needed  to  get  better 
relations  among  winds,  wind-wave  interactions,  wind  stresses)  and  surface  heat 
fluxes.  SateUite  data  sets  are  available  for  sea-surface  heights  and  temperatures, 
but  the  temperatures  probably  require  correction  for  air  humidity  because  of  the 
low  amplitude  temperature  signal  in  the  Mediterranean.  Model-predicted  sea 
surface  temperature  can  be  used  to  distribute  atmospheric  model  and  satellite  heat 
exchange  data  appropriately  to  small  .Vlediterranean  scales.  The  new  high- 
resolution  NORAPS  data  will  be  most  welcome.  Meantime,  existing  multi-year 
data  sets  from  ECMWF  can  be  used.  Drifters  are  more  useful  for  model  validation 
than  for  assimilation. 
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Professor  Joe  Thompson 
Engineering  Research  Center 
Mississippi  State  University 


The  use  of  numerically-generated  curvilinear  coordinate  systems  in  partial 
differential  equation  (PDE)  solutions  allows  finite  difference  methods  to  treat 
irregular  boundaries  without  interpolation  or  stair-stepping,  since  some 
coordinate  surface  will  be  coincident  with  each  boundary  segment.  This  allows 
difference  representations  and  boundary  conditions  to  be  implemented  easily, 
regardless  of  the  shape  of  the  boundaries  of  the  region.  The  entire  problem  is,  in 
fact,  converted  to  one  that  is  logically  rectangular,  i.e.,  a  solution  on  an 
equally-spaced  Cartesian  grid  on  a  region  composed  of  rectangular  blocks. 

Thus,  solution  algorithms  can  operate  on  an  equally-spaced  rectangular 
grid,  regardless  of  the  shape  of  the  region.  It  is  also  possible  to  vary  the 
distribution  of  the  grid  lines  in  order  to  reduce  the  spacing  in  regions  of  large 
solution  variation,  either  statically  beforehand  or  through  dynamic  coupling  with 
the  developing  solution. 

The  price  to  be  paid  in  the  use  of  such  block-structured  grids  is  that  the 
PDE’s  acquire  additional  terms,  and  all  directional  components  in  the  equations 
are  coupled.  For  instance,  the  first  derivative  becomes,  in  two  dimensions, 

f  =  Vn  ~  ye 
^  Xe  yti  ■  Xti  Ye 

Stencils  that  are  5-point  in  2-D  on  orthogonal  grids  become  9-point  on 
general  nonorthogonal  systems.  Similarly,  7-point  stencils  in  3-D  become 
27-point. 

Also,  any  advantage  taken  of  the  x-momentum  equation  being  dependent 
only  on  the  x-component  of  the  pressure  gradient,  as  in  some  staggered-grid 
approaches,  is  lost  in  the  general  curvilinear  system,  since  in  2-D  the 
X-momentum  equation  will  contain  both  pe  and  prj.  Such  staggered-grid  codes 
are  not  directly  convertible  to  use  the  general  curvilinear  systems.  However, 
other  codes  can  be  converted  without  unreasonable  effort,  since  all  that  is  really 
involved  is  the  addition  of  terms  to  the  expressions  and  some  data  handling. 


A-155 


The  principal  advantages  of  general  curvilinear  systems  come  into  play 
when  (1)  boundary  shape  is  important  to  the  solution,  or  (2)  there  is  the  need  to 
resolve  strong  gradients.  A  third  consideration  can  be  storage,  since  all  points 
are  utilized  and  none  are  outside  the  boundary. 

The  Engineering  Research  Center  at  Mississippi  State  University  has 
experience  in  converting  and  developing  codes  for  general  curvilinear  systems 
and  can  assist  code  developers  in  the  effort. 


A- 156 


APPENDIX  B 


Descriptive  Oceanography 

This  appendix  includes  a  list  of  some  of  the  relevant  dynamical 
considerations  in  the  Mediterranean. 

1.  North  African  current  baroclinically  unstable.  -1.2  Sv  inflow. 

2.  Ligurian  Sea  flow:  x  -  driven. 

3.  Straits  of  Sicily  divides  basins,  modified  Atlantic  water  penetrates  to  the 
east. 

(I) 

4.  Topographic  control  \^l  important  everywhere. 

5.  (V  X  x)  z  probably  important  in  Ionian  Basin  and  west. 

6.  ~  1.0  Sv  of  thermohaline  circulation  exits  Gibraltar  at  intermediate  depths; 
0.2  Sv  evaporates. 

7.  Bottom  water  formation  occurs  in  South  Adriatic.  Does  it  also  occur  in 
Ligurian  or  Tyrrhenian?  How  to  capture  in  models? 

8.  Interactions  between  basins,  sub-basins,  and  mesoscale  flows. 

9.  Surface  exchanges  (waves,  etc.). 
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